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Abstract 


Defense  suppression  of  enemy  ground  forces  is  basic 
to  successful  counterair  objectives,  ^The  F-4G  Wild  Weasel 
(WW)  weapon  system  provides  the  teeth  in  getting  the  defense 
suppression  job  done — identifying,  locating,  and  killing 
enemy  ground  based  threat  radars.  The  objective  of  this 
thesis  was  to  develop  a  methodology  that  could  examine  and 
evaluate  the  WW  defense  suppression  mission.  The  problem 
was  developed  for  a  NATO/Warsaw  Pact  encounter  in  Central 
Europe . 

A  model  of  the  threat  environment  was  built  using 
the  SLAM  computer  simulation  language.  Threats  in  the 
defense  sector  can  be  moved  as  desired.  Friendly  aircraft 
can  enter  the  system  at  a  variety  of  intervals,  altitudes, 
and  airspeeds.  WWs  hunt  for  threats  to  attack  by  search¬ 
ing,  identifying,  locating,  and  then  launching  their 
weapons  at  the  threat.  WW  tactics  can  be  changed  as  the 
requirements  of  the  mission  dictate  or  at  the  desire  of 
the  WW  crew.  Self -protect ion  jamming  can  be  selected  by 
either  WW  or  attack  aircraft.  Enemy  threats  will  fire  at 
an  aircraft  when  the  aircraft  comes  within  the  threat's 
range  as  long  as  the  threat  is  not  engaged  with  another 
aircraft.  Early  warning  radars  account  for  threat  radar 


^ command  and  control  functions;  their  control  over  the 
associated  radars  can  be  changed  as  desired.. s 

Changing  the  WW's  altitude  from  60  meters  to  200 
meters  did  not  effect  friendly  attack  aircraft  survivabil¬ 
ity.  Leading  the  attack  force  into  the  threat  area  by 
30  seconds  as  opposed  to  accompanying  the  attack  force  did 
not  influence  attack  force  survivability.  Further  develop 
ment  of  the  model  to  include  turn-mode  capability  for  the 
WW  weapon  and  a  tactic  for  pre-emptive  weapons  launch  in 
anticipation  of  threat  radar  radiation  is  recommended. 
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A  WILD  WEASEL  PENETRATION  MODEL 


I .  Introduction 


Background 

The  requirement  for  a  weapon  system  that  would  deal 
exclusively  with  the  surface-to-air  missile  (SAM)  threat 
developed  from  air  battles  in  the  Vietnam  War.  The  North 
Vietnamese  used  Soviet  SA-2  SAMs  in  concert  with  AAA  and 
MIG  aircraft  to  counter  a  numerically  superior  air  force. 

To  counter  these  tactics,  the  USAF  developed  the  F-100/ 

F-105  Wild  Weasel  (WW) ,  a  weapon  system  dedicated  to  the 
anti-SAM  mission. 

The  current  strategy  of  the  Soviet  Union  and  Warsaw 
Pact  countries  emphasizes  air  doctrine  through  their  ground 
forces  to  gain  superiority  on  the  battlefield  (Ref  6:69). 
Their  extensive,  sophisticated  mobile  air  defenses,  consist¬ 
ing  of  mixes  of  guns  and  missiles,  provides  overlapping 
coverage.  The  net  effect  is  a  wide-ranging,  protective 
umbrella  for  their  ground  forces. 

NATO's  task  in  a  confrontation  against  the  Warsaw 
Pact  will  be  the  ultimate  test  of  its  airpower  capabili¬ 
ties.  Counterair  operations  must  gain  air  superiority  over 
the  battlefield  if  ground  forces  are  to  be  effective. 


USAF  Basic  Doctrine/  AFM  1-1,  recognizes  this  important 
requirement  and  now  includes  defense  suppression,  along 
with  offensive  and  defensive  counterair  operations,  as  a 
primary  task  within  the  counterair  mission  (Ref  1:2-16). 
Defense  suppression  is  a  fundamental  objective  for  an 
effective  air-ground  force. 

The  WW  weapon  system  represents  a  key  element  in 
the  defense  suppression  mission.  The  F-4G  aircraft,  the 
current  WW  used  by  the  USAF,  can  accomplish  defense  suppres 
sion  objectives  against  threat  radar  systems  by  either 
physically  destroying  the  radars  with  anti-radiation 
missiles  (ARMs)  or  bombs,  or  by  causing  the  radars  to 
cease  operation  as  a  precaution  against  an  attack  by  the  WW 

If  the  F-4G  WW  is  the  ideal  airpower  instrument  for 
suppressing  enemy  air  defenses  then  the  question  remains 
as  how  best  to  use  this  tool  to  maximum  effectiveness.  In 
light  of  the  fact  that  both  the  F-4G  and  its  specialized 
ARM  weapons  are  very  limited  resources,  WW  operations  must 
be  effective  preventing  threat  radars  from  attacking 
friendly  aircraft  while  at  the  same  time  surviving  attack 
from  these  very  same  radars.  The  highly  fluid  arena  of  the 
air-ground  battle  coupled  with  the  expected  heavy  concen¬ 
tration  of  enemy  radar  threats  and  the  diversity  of  their 
employment  make  WW  operations  a  complex  and  difficult  task. 

In  planning  for  a  defense  suppression  mission  not 
only  must  many  different  factors  be  considered  but  also  the 


underlying  relationships  between  them  clearly  understood. 
Aspects  such  as  ASMs  employment/  force  sizing,  ingress/ 
egress  routing,  threat  hierarchy,  electronic  countermea¬ 
sures,  and  command  and  control  must  be  carefully  thought 
out  and  planned  for.  In  addition,  defense  suppression 
operations  must  be  amenable  to  last-second  changes  as  a 
result  of  the  evolving  air  battle. 

Heedless  to  say,  choosing  the  best  tactic  or  best 
set  of  tactic  options  for  WW  operations  is  not  easy.  None¬ 
theless,  if  counterair  operations  are  to  be  successful 
then  neutralizing  enemy  air  defenses  must  be  accomplished 
by  WWs  using  optimum  defense  suppression  tactics. 

An  analysis  for  the  most  advantageous  weapons 
allocation  and  tactics  that  explicitly  evaluates  the  WW  as 
a  complete  system  has  yet  to  be  accomplished.  Although 
no  single  weapon  system  can  be  expected  to  successfully 
counter  every  aspect  of  a  highly  sophisticated  threat 
environment,  the  F-4G's  WW  defense  suppression  effective¬ 
ness  is  a  major  factor  upon  which  the  USAF's  counterair 
mission  rests. 

Problem  Statement 

The  ability  of  WWs  to  perform  defense  suppression 
operations  in  a  radar-rich  threat  environment  is  basic  to 
successful  counterair  operations.  The  objective  of  this 
thesis  is  to  develop  a  methodology,  through  a  simulation 
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model,  for  evaluating  the  WW  defense  suppression  mission. 

In  particular,  the  methodology  should  be  capable  of  ana¬ 
lyzing  force  sizing,  ARM  conf igurations ,  and  WW  tactics. 

These  parameters  were  selected  because  they  are  considered 
significant  areas  that  impact  on  the  WW  mission  (Refs  10;  11). 

Assumptions  and  Limitations 

Conclusions  from  this  thesis  can  only  be  applied 
within  the  context  of  the  developed  system.  Thus,  because 
of  the  complexity  of  the  WW  defense  suppression  operations, 
not  all  components  and  variables  of  the  system  are  included. 

It  should  be  noted  that  only  WW  defense  suppression  opera¬ 
tions,  and  not  any  other  WW  mission,  are  analyzed  in  this 
research  effort. 

The  following  major  assumptions  and  limitations 
apply  to  the  system. 

1.  WWs  carry  only  ARMs.  For  other  types  of  mis¬ 
sions  WWs  may  operate  in  a  mixed-mode  configuration  (ARMS 
and  hard  bombs)  or  with  only  hard  bombs.  For  the  escort- 
type  scenario  in  this  thesis  only  ARMs  are  used. 

2.  Self-protection  jamming  is  used  only  by  attack 
aircraft.  Jamming  pods  for  WWs  ar.e  available  but  are  not 
used. 

3.  Once  a  WW  begins  an  attack,  the  attack  is  con¬ 
tinued  until  the  WW  launches  an  ARM  at  the  threat  or  the 
WW  is  killed. 
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4.  wws  have  perfect  inter-aircraft  radio  communica- 


tions.  Although  this  assumption  is  open  to  debate,  secure 
voice  radio  equipment  may  insure  it. 

Results  of  the  thesis  are  only  useful  for  making 
relative  comparisons  between  alternatives  evaluated.  When 
five  out  of  ten  aircraft  are  predicted  to  survive  for  a  par 
ticular  simulation  run,  given  one  set  of  conditions,  and 
only  one  survives  for  another  set  of  conditions,  the  impor¬ 
tant  result  is  the  comparison  between  the  alternatives  not 
the  number  of  surviving  aircraft. 

Threat  Scenario 

A  WW-threat  environment  was  developed  based  on 
confrontation  between  NATO  and  Warsaw  Pact  in  Central 
Europe.  The  threat  consists  of  a  typical  Soviet  ground 
army  deployed  along  the  forward  edge  of  the  battle  area 
(FEBA) .  A  force  of  WWs  are  assigned  the  task  of  defense 
suppression  in  support  of  a  low-level  fighter  attack  force 
that  will  penetrate  the  FEBA  and  fly  through  the  enemy 
defenses  to  strike  a  target  behind  enemy  lines.  For  the 
scenario  the  WWs  can  realize  their  defense  suppression 
objectives  by  destroying  ground-based  threat  radars  or 
associated  early  warning  (EW)  radars. 

The  air  defense  elements  of  a  typical  Soviet  Army, 
located  within  an  area  50  kilometers  wide  by  100  kilometers 
long,  consists  of  approximately  1000  SAMs  and  AAA  units. 
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Concentrated  close  to  the  FEBA  will  be  hundreds  of  small 
arms.  However,  weapons  that  require  optical  information 
to  acquire  and  track  targets  are  marginally  effective 
against  high  speed,  maneuvering  aircraft.  In  addition, 
Soviet  interceptors  are  hypothesized  not  to  be  a  factor  for 
the  low  altitude  attack  force  as  their  doctrine  prescribes 
operation  along  the  FEBA  at  or  above  10,000  feet,  leaving 
the  ground  forces  responsibility  for  the  air  superiority 
mission  up  to  this  altitude  (Ref  6:70). 

For  these  reasons  the  following  defense  elements 
were  selected  as  potential  threats  to  the  attack  force. 

1.  AAA 

2 .  SAM-A 

3.  SAM-B 

4.  SAM-C 

5.  SAM-D 

Command  and  control  of  threat  radars  is  represented  by  EW 
radars.  For  this  particular  threat  environment  there  are 
eight  of  these  radars  positioned  in  the  defense  sector. 

Structural  Model 

Figure  1  is  the  structural  model  of  the  air 
defense  developed  within  the  constraints  of  the  scenario. 
Threats  are  located  in  belts  that  approximate  their 
expected  position.  Their  actual  position,  however,  will 
depend  on  battlefield  tactics,  terrain  features,  employment 
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SAM  A:  2  launchers/battery 
SAM  B,  C,  D:  3  launchers/battery 
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doctrine,  etc.  Note  that  the  number  of  each  threat  is 
specified  in  Figure  1  but  not  exact  location. 

The  enemy's  area  of  operation  behind  the  FEBA  is 
modeled  by  a  single  line  of  communication  (LOC)  road  net¬ 
work  extending  from  the  target  area  to  the  FEBA.  The 
LOC 1 s  layout  depends  on  the  area ' s  terrain  features  and 
the  requirements  to  carry  equipment  and  supplies  to  the 
front  line.  The  attack  force  will  avoid  any  major  LOC  due 
to  the  probable  concentration  of  weapons  along  it.  Accord¬ 
ingly,  the  LOC  for  the  scenario  is  hypothesized  to  lie  at 
the  upper  portion  of  the  defense  area.  Figure  2  depicts 
a  typical  LOC. 

Methodology 

WW  defense  suppression  operations  are  considerably 
complex  involving  many  dynamic  component  interactions.  Sets 
of  threat  radars  search  for,  acquire,  track,  and  attack 
aircraft  that  fly  through  their  defense  sector.  WWs  search 
for,  identify,  locate,  and  launch  their  weapons  at  threat 
radars  while  maneuvering  in  the  battle  area.  Attack  air¬ 
craft  penetrate  then  fly  through  the  threat  scenario  enroute 
to  a  target  far  behind  enemy  lines. 

Because  of  the  complexity  of  the  WW  system  and  the 
need  to  study  the  intricate  interplays  of  the  system  com¬ 
ponents,  a  computer  simulation  model  of  the  WW  defense 
suppression  mission  was  developed.  Simulation  was  chosen 
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Structural  Model  and  Lines  of  Communication 


because  the  problem  lends  itself  to  this  methodology. 
Shannon  lists  several  conditions  when  simulation  should  be 
considered  for  use;  two  of  these  cast  a  strong  vote  for  the 
defense  suppression  problem. 

1.  A  complete  mathematical  formulation  of  the 
problem  does  not  exist  or  analytical  methods  of  solving 
the  mathematical  model  have  not  yet  been  developed. 

2.  Simulation  may  be  the  only  possibility  because 
of  the  difficulty  in  conducting  experiments  and  observing 
phenomena  in  their  actual  environment  (Ref  14:11). 

An  added  and  perhaps  equally  important  advantage 
is  that  with  simulation  the  complex  system  under  investiga¬ 
tion  can  be  expressed  in  elementary  events  (Ref  12:3). 

By  formulating  the  simulation  and  thus  the  problem  in  terms 
of  elementary  events  complex  interactions  such  as  those  in 
the  WW  defense  suppression  system  can  be  synthesized. 

SLAM  was  chosen  as  the  simulation  language  because 
of  the  power  and  flexibility  it  offered.  The  dynamic  por¬ 
trayal  of  the  WW  attack  as  well  as  the  threat  radar- 
aircraft  interactions  were  accurately  represented  through 
the  use  of  SLAM's  language  and  concept  of  process  orienta¬ 
tion.  These  features  will  be  discussed  in  detail  in 
Chapter  III . 
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Overview 


The  thesis  is  explained  in  detail  in  the  following 
chapters.  The  system  structure  of  the  model  and  the  basic 
components  that  make  up  the  system  are  discussed  in 
Chapter  II.  The  simulation  model  is  explained  in  a  logical, 
sequential  manner  in  Chapter  III.  In  Chapter  IV,  data 
collection,  experimental  design,  the  results  of  the  experi¬ 
ment,  and  the  validation  of  the  experiment  are  all  discussed. 
Finally,  the  results  of  the  thesis  are  covered  in  Chapter  V 
and  recommendations  for  follow-on  areas  are  listed  in 
Chapter  VI. 
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II.  Systems  Structure 


Introduction 

The  purpose  of  this  thesis  is  to  develop  a  model 
for  investigating  the  interaction  between  a  friendly  air¬ 
craft  force  penetrating  the  forward  edge  of  the  battle  area 
•  (FEBA)  to  strike  targets  in  a  rear  area  and  the  enemy  air 

t 

defense  network  responsible  for  preventing  this  penetration. 
Two  elements  comprise  the  friendly  force:  a  strike  force  of 
fighter  aircraft  and  WW  aircraft  employed  in  a  defense  sup¬ 
pression  mission.  The  air  defense  network  is  composed  of 
SAM,  AAA,  and  early  warning  (EW)  systems.  Possible  inter¬ 
action  between  the  elements  include  the  following: 

1.  The  WW  attack  and  attempt  to  neutralize  any  of 
the  air  defense  systems.  Neutralize  includes  either 

,  physically  destroying  the  site  or  stopping  them  from  radi¬ 

ating. 

2.  The  air  defense  systems  attempt  to  destroy 
both  strike  and  WW  aircraft  penetrating  the  FEBA. 

The  model  will  be  used  to  study  specific  strike 
force  and  WW  tactics  such  as  ingress  altitude,  airspeed, 
and  spacing  between  aircraft  with  the  ultimate  intent  of 
developing  procedures  designed  to  increase  the  overall 
probability  of  the  strike  force  successfully  penetrating  to 
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the  target  area.  The  model  must  be  flexible  enough  to 
allow  both  the  size  and  mix  of  both  friendly  and  enemy 
forces  to  be  varied. 

Pertinent  to  the  model's  development  are  the  major 
system  elements  required  to  describe  the  system.  These 
elements  are  the  following: 

1.  The  WW  aircraft,  its  associated  characteristics 
and  tactics, 

2.  The  strike  force,  and 

3.  The  enemy  defensive  network. 

This  chapter  discusses  each  element  separately  and  in 
detail.  It  explains  the  analytical  methodology  required 
by  each  element  to  perform  its  mission. 

Wild  Weasel  Mission  Scenario 

Today's  WW  platform  is  a  modified  F4-G  aircraft 
containing  sensitive  radar  homing  and  warning  (RHAW)  equip¬ 
ment.  The  WW  crew's  mission  is  to  search  the  enemy's  area 
of  operation,  detect  and  identify  enemy  radar  signals 
associated  with  SAMs,  AAA,  or  EW  systems,  and  force  these 
radars  to  cease  operation  by  either  physically  destroying 
the  site  with  bombs  or  antiradiation  missiles  (ARMs)  or  by 
causing  the- sites  to  stop  radiating  to  preclude  being 
attacked  by  a  WW. 

Leek  and  Schmidt  investigated  the  FEBA  penetration 
and  the  associated  enemy  defense  network  by  a  flight  of 


strike  aircraft  and  evaluated  the  loss  rates  incurred  in 
the  penetration.  This  thesis  extends  their  analysis  by 
including  WW  aircraft  employed  in  a  defense  suppression 
mission  and  again  evaluates  the  loss  rates  of  both  strike 
and  WW  aircraft.  In  this  scenario  specific  threat  loca¬ 
tions  will  not  be  known  before  the  mission  by  the  WW  crews 
although  estimates  of  the  number  and  type  of  enemy  radars 
will  be  provided  by  intelligence  briefings.  The  WW  knows 
the  corridor  used  by  the  strike  force  during  ingress  to  the 
target.  This,  along  with  projected  capabilities  of  the 
enemy  radars,  allows  the  WW  to  estimate  an  approximate  dis¬ 
tance  either  side  of  the  attack  corridor  a  specific  SAM  or 
AAA  can  be  located  and  constitute  a  threat  to  the  strike 
force . 

Figure  3  depicts  the  four  phases  of  a  typical  WW 
mission  profile:  search,  detection,  ranging,  and  attack. 

The  figure  illustrates  the  dynamic  nature  of  the  scenario. 
For  the  mission,  the  WW  aircraft  assemble  50  km  prior  to 
the  FEBA  and  over  friendly  airspace.  In  the  search  phase 
the  WW  proceeds  toward  the  enemy  defenses,  its  RHAW  equip¬ 
ment  sweeping  through  radar  frequencies  attempting  to  detect 
and  identify  enemy  radar  signals.  At  detection,  the  RHAW 
system  alerts  the  crew  to  the  type  of  threat,  the  bearing 
to  the  site  from  the  aircraft,  and  an  approximate  distance 
to  the  threat  radar.  The  WW's  onboard  processor  will  refine 
the  site's  location  to  within  a  few  meters  during  the 
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ranging  phase.  Here,  the  system  calculates  numerous  lines 
of  position  (LOP)  between  the  aircraft  and  the  radar  site 
as  the  aircraft  flies  through  approximately  30  degrees 
of  bearing  change.  The  site's  location  is  estimated  to  be 
the  intersection  of  the  LOPs.  The  more  LOPs,  the  more  pre¬ 
cise  the  site's  location.  After  completing  the  ranging 
routine,  the  aircraft  enters  the  attack  phase  by  turning 
toward  the  site,  and  for  this  profile,  launching  an  ABM. 
After  launch  the  aircraft  enters  the  next  search  sequence 
and  the  process  repeats  until  the  WW  depletes  all  its  ARMs. 

Implicit  within  the  profile  description  are  several 
areas  that  must  be  determined  during  the  mission.  Did  the 
search  phase  start  while  the  threats  were  beyond  the  field 
of  view  (FOV)  of  the  aircraft  or  when  the  sites  were  inside 
the  FOV?  Did  the  WW  maneuver  in  a  direction  that  allowed 
it  to  minimize  the  amount  of  time  required  to  complete  each 
profile?  What  happens  if  the  attack  phase  occurs  inside 
the  minimum  launch  range  of  an  ARM?  These  questions  and 
the  logic  required  to  evaluate  them  will  be  addressed  in  the 
following  paragraphs. 

Threat  Detection  Criteria 

The  WW  search  phase  begins  well  before  the  FEB A. 


The  minimum  distance  the  WW  can  detect  specific  ground- 
based  emitters  is  limited  by  either  the  RHAW  receiver's 
sensitivity  or  the  radar  horizon  (or  FOV)  of  the  aircraft. 


Receiver  sensitivity  indicates  the  minimum  power  signal 
at  the  aircraft's  RHAW  receiver  that  can  be  processed  by 
the  equipment.  In  certain  instances  either  sensitivity  or 
limited  FOV  dominates  the  other  in  determining  the  maximum 
detection  range.  Which  one  dominates  in  the  WW  case  will 
now  be  determined. 

An  aircraft's  radar  horizon  depends  on  its  alti¬ 
tude.  The  earth's  curvature  limits  the  radar  horizon  of 
either  a  transmitter  or  receiver.  (See  Figure  4.)  The 
horizon  can  be  approximated  by  the  following  equation: 


where  h  *  aircraft's  altitude,  feet;  and 
R  »  radar  horizon,  NM  (Ref  5:36). 

Converting  both  R  and  h  to  meters,  the  equation 

becomes: 


R  «  4117.3  (2) 

where  h  *  aircraft’s  altitude,  m;  and 
R  ■  radar  horizon,  m. 

Thus,  an  emitter  at  a  height  of  60  meters  could 
detect  objects  out  to  32,000  meters  (disregarding  clutter 
and  terrain  blockage) .  Objects  at  some  height  above  the 
earth's  surface  beyond  this  32,000  meter  range  could  extend 
this  detection  distance.  For  example,  the  same  emitter  at 
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Fig.  4.  Radar  Horizon 


60  meters  could  detect  an  object  at  a  height  of  30  meters 
out  to  54,000  meters  (again  disregarding  clutter  and  ter¬ 
rain  blockage) .  For  ground  based  emitters  with  antenna 
extending  less  than  10  meters  above  the  surface,  the  radar 
horizon  calculations  of  equation  (2)  appear  sufficient  for 
this  model. 

Power  received  at  the  aircraft,  Pr,  is  propor¬ 
tional  to  the  radar's  effective  radiated  power  reaching 
the  WW  aircraft  (ptGt)  and  the  WW  effective  antenna  area, 

A  ,  where  P  is  the  transmitted  power  in  watts  and  G.  is 
the  transmitting  antenna  signal  gain  in  the  direction  of 
the  WW  aircraft.  A  ,  in  square  meters,  can  be  approxi- 
mated  as  follows: 


where  c  *  speed  of  light,  meters  per  second  (mps) ; 


f  =  radar  transmit  frequency,  in  Hertz  (Hz);  and 
Gr  *  WW  receiver  antenna  gain  (Ref  5:18). 

For  the  case  of  the  WW  receiver  antenna,  the  gain 
will  be  assumed  to  be  one,  or  the  antenna  is  omni¬ 
directional.  Power  received  at  the  aircraft  is  inversely 
proportional  to  the  receiver's  internal  losses,  L,  and  the 
square  of  the  distance  between  the  site  and  the  aircraft, 
R.  Expressed  mathematically,  Pr  becomes: 


P 


r 


P  G 

c-s-fe  (f ) 

4TTR^  1 


(4) 


where  all  terms  are  as  previously  defined  (Ref  5:18) . 

A  RHAW  receiver  is  designed  to  a  sensitivity  speci¬ 
fication  determined  by  the  type  of  signals  radiated  in  the 
receiver's  working  environment.  Once  designed,  the  minimum 
sensitivity  cannot  be  changed.  If  the  receiver  is  too 
sensitive  the  processor  will  be  saturated  by  too  many 
signals;  not  sensitive  enough  and  many  threats  will  not  be 
detected  by  the  WW  until  the  aircraft  is  too  close  to  the 
site.  Receiver  sensitivity  is  expressed  in  decibels  (dB) 
referenced  to  a  specific  power  level.  For  this  thesis, 
the  reference  will  be  1  watt  and  the  sensitivity  expressed 


as  dB,  watts  or  dBw.  Converting  a  parameter  such  as 
radiated  power  to  dBw  is  accomplished  as  follows: 


P*  *  10  log1Q(Pt  ) 


'dBw 


(5) 


"watts 


Converting  all  of  the  terms  of  the  Pr  equation  yields: 


rdB  tdB 


+  G*.  +  2(0)^  -  2  (4 it) 


"dB 


dB 


-  2<R>dB  ‘  LdB  -  2<f>dB 


(6) 


For  a  given  aircraft/radar  encounter  all  terms  in 
the  Pr  calculations  can  be  assumed  to  be  constant  except 
Pr,  R,  and  Gfc.  Table  I  depicts  typical  values  for  Pr,  f , 
and  L  in  the  FEBA.  Rearranging  equation  (6)  and  solving 
for  Gfc  in  terms  of  Pfc  and  R  yields: 


*  2  (R)  Jt»  -  P*  +  2(4tt)._  +  P^ 
fcdB  63  fcdB  63  rdB 


+  LdB  +  2(f)dB  -  2{c)dB 


(7) 


where  all  terms  are  as  previously  defined. 


TABLE  I 

AVERAGE  SYSTEM  PARAMETERS 


Parameter 

Normal  Units 

Decibels  (dBw) 

Transmit  Power,  Pr 

50-100  kw 

47-50 

Radar  Frequency,  f 

10-15  GHz 

100-102 

System  Losses,  L 

- 

10 

20 


Setting  Pr  to  the  receiver  sensitivity,  S,  at  the  air¬ 
craft,  equation  (7)  becomes  an  equation  for  in  terms 
of  R  alone.  A  receiver  sensitivity  of  -100  dBw  is  both 
technologically  feasible  and  operationally  effective. 
Setting  Pr  to  -100  dBw  equation  (7)  becomes: 

G.  »  2 (R) -  80.6  dB  (8) 

fcdB  63 

As  a  limit  of  the  Gt  required  to  process  a  signal  assume 
R  is  the  radar  horizon  of  an  aircraft  flying  at  an  altitude 
of  100  m  (R  *  41,173  m  or  46.15  dB) .  Solving  equation 
(8)  for  Gfc  yields: 

G.  -  92.3  -  80.6  *  11.7  dB  (9) 

rdB 

Typical  values  for  the  main  beam  Gfc  near  the  FEBA  are 
40  dB  with  an  average  sidelobe  gain  20  dB  down  from  this 
maximum,  or  the  average  sidelobe  Gfc  of  20  dB.  Since  the 
receiver  requires  only  11.7  dB  {or  28.3  dB  below  maximum) 
at  the  radar  horizon,  the  system  could  identify  radars 
when  receiving  sidelobe  level  energy  signals  or,  stated 
another  way,  the  receiver  could  process  radar  signals  from 
an  emitter  where  the  antenna  main  beam  is  randomly  oriented 
with  respect  to  the  WW.  For  this  reason,  the  aircraft's 
FOV  and  not  the  receiver  sensitivity  limits  the  WW  in 
detecting  threat  radars  in  this  model. 


Wild  Weasel  Heading  System 

The  continuous  time  model  allows  the  WW  aircraft 
to  maneuver  in  response  to  a  threat.  A  heading  reference 
system  (HRS)  similar  to  an  aircraft's  horizontal  situation 
indicator  or  compass  was  developed  to  follow  an  aircraft 
throughout  its  mission.  By  using  this  HRS  the  key  elements 
of  the  WW's  low  altitude  mission  could  be  determined.  Four 
quantities  defined  by  the  HRS  are  listed  below  and  described 
in  the  following  paragraphs: 

1.  The  WW's  heading,  H; 

2.  The  magnetic  bearing  from  the  aircraft  to  the 

site,  B; 

3.  The  absolute  bearing,  AB?  and 

4.  The  relative  bearing. 

The  WW's  heading,  H,  is  the  direction  of  the  air¬ 
craft's  flight  vector  referenced  to  North  or  000  degrees 
and  is  calculated  as  follows.  (Note:  the  model  assumes 
a  no-wind  condition,  thus  the  aircraft's  flight  vector, 
velocity  vector,  and  heading  are  aligned.)  An  initial 
heading,  H^,  and  a  rate  of  heading  change  with  time,  AH, 
are  initially  specified.  A  left  hand  turn  results  in 
negative  AH  and  a  right  hand  turn,  a  positive  AH.  The  new 

heading,  H  ,  At  time  units  later  is  determined  as  follows: 
n 

H  -  H.  +  (AH) (At)  (10) 

n  l 
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A  360  degree  correction  factor  limits  the  value  of 


H  to  0-360  degrees.  For  example,  if  the  WW's  heading,  , 
was  010  degrees  and  a  left  hand  turn  of  four  degrees  per 
second  for  five  seconds  is  required,  equation  (10)  calcu¬ 
lates  a  new  heading,  Hr,  of  -10  degrees.  The  correction 
factor  of  +  360  degrees  is  applied  to  negative  headings, 
and  H  equals  350  degrees. 

The  magnetic  bearing,  B,  locates  the  threat  radar 
site  with  respect  to  the  aircraft  and  as  was  true  with  H, 
is  referenced  to  North  or  000  degrees.  B  depends  on  the 
relative  position  of  the  aircraft  with  respect  to  the  radar 
Figure  5  depicts  the  airspace  around  an  emitter,  S,  divided 
into  four  quadrants.  For  each  quadrant  the  angle  0  is 
calculated  as  follows: 


0 


(11) 


where  Yg  *»  y  coordinate  of  the  threat  radar, 

Xg  *  x  coordinate  of  the  threat  radar, 

yA  »  y  coordinate  of  the  aircraft,  and 

xA  *  x  coordinate  of  the  aircraft. 

Knowing  the  magnitude  of  the  Tan  0  and  the  magni¬ 
tude  of  the  quantity  (Yg  -  yft)  the  aircaft.  A,  is  located 
in  a  specific  quadrant  around  the  threat  radar.  For 
example,  consider  an  aircraft  located  in  quadrant  IX. 

The  quantity  (Y_  -  Y.)  is  positive,  and  the  quantity 


Fig.  5.  Magnetic  Bearing  Determination 


(Xg  -  xA)  is  negative.  Tan  9  will  be  negative.  A  nega¬ 
tive  tan  9  and  a  positive  (Yg  -  yft)  occurs  only  in  quadrant 
II.  B  is  then  referenced  to  000  degrees  by  either  adding 
0  to  or  subtracting  0  from  one  of  two  cardinal  headings — 
090  or  270  degrees.  The  B  calculation  for  each  quadrant 
around  site  S  is  also  depicted  in  Figure  5. 

The  absolute  bearing,  AB,  is  the  abaolute  value  of 
the  difference  between  the  heading,  H,  and  the  magnetic 
bearing,  B,  or: 

AB  *  Absolute  value  (H-B)  (12) 

All  of  the  WW's  low  altitude  ranging  routines  are 
simulated  by  calculating  the  value  of  AB.  For  instance, 
after  detecting  a  threat  radar  and  with  both  H  and  B 
known,  the  WW  will  start  turning  until  AB  reaches  75 
degrees,  the  hypothesized  value  of  AB  required  to  start 
the  WW's  LOP  calculations.  The  WW  will  stop  the  turn  and 
allow  AB  to  increase  to  105  degrees  to  simulate  the  entire 
ranging  routine. 

As  will  be  shown,  AB  is  used  in  the  relative  bear¬ 
ing  calculations  and  is  limited  to  a  maximum  of  180  degrees. 
For  computed  values  greater  than  180  degrees,  a  new  value 
is  calculated  as  follows: 

AB  »  360  -  AB  (13) 
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Knowing  H,  B,  and  hence,  AB,  the  WW  can  accomplish  its  low 
altitude  attack  profile  in  both  a  minimum  time  and  air¬ 
space.  Also,  by  knowing  these  terms,  a  fourth  quantity  of 
the  BBS  is  specified — the  relative  bearing,  RB.  RB  is  the 
bearing  of  the  site  from  the  aircraft  referenced  to  the 
aircraft's  current  heading.  If  (H  >  B) ,  the  site  is 
located  AB  degrees  to  the  left  of  the  aircraft.  If  (B  >  H) , 
the  site  is  located  AB  degrees  to  the  right  of  the  aircraft. 
Subsequent  turns  away  from  and  toward  the  site  will  be  made 
in  reference  to  the  current  RB.  Figure  6  depicts  an  air¬ 
craft  heading  090  degrees  with  a  RB  of  045  degrees  to  the 
site,  S.  Since  (H  >  B)  ,  the  site  is  located  (90  -  45)  or 
45  degrees  to  the  left  of  the  aircraft. 

Slant  Range  to  Site  Calculations 

After  developing  a  HRS,  the  analytical  model  con¬ 
tinued  by  calculating  the  slant  range,  SR,  between  the 
radar  site  and  the  WW  aircraft  for  each  time  increment 
during  the  simulation  run.  Solving  for  SR  by  triangula¬ 
tion: 

SR  -  V (Ys‘yA) *  +  (Xs"xA) 2  +  ALt2  (14) 

where  ALT  -  the  aircraft's  altitude,  meters  (m) ,  and 
other  terms  as  previously  defined. 

As  the  aircraft  proceeds  toward  the  FEBA,  its  x  and 
y  position  change  continuously,  and  the  SR  to  a  specific 
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H  090  degrees 


— - — -B  045  degrees 

Fig.  6.  Aircraft/Site  Relationship 

site  also  changes.  When  the  SR  to  a  site  decreases  to  less 
than  the  WW's  FOV  the  model  assumes  the  site  has  been 
detected.  In  the  model  detection  includes  both  determining 
the  type  of  radar  and  its  RB  from  the  aircraft's  current 
position.  After  detecting  the  threat,  the  WW  maneuvers 
to  a  position  from  which  the  ranging  phase  of  the  profile 
can  be  accomplished.  For  this  thesis,  the  aircraft  will 
turn  in  the  shortest  direction  to  achieve  a  75  degree  AB. 
For  an  operational  WW  system  the  entire  ranging  routine 
requires  continuous  LOP  calculations  for  a  30  degree 
change  in  AB.  Thus,  the  ranging  requires  AB  to  increase 
from  75  to  105  degrees.  Figure  7  depicts  the  reason  75 


Wild  Weasel  Ranging  Problem 


degrees  AB  was  chosen  as  the  initial  rollout  bearing.  An 
aircraft  would  need  to  fly  the  distance  EF  to  achieve  a 
30  degree  AB  change  if  the  initial  AB  was  30  degrees,  but 
a  shorter  length  leg  (EC)  for  an  initial  AB  of  90  degrees. 
By  making  the  required  AB  change  symmetrical  about  the  90 
degree  AB,  a  minimum  ranging  time  duration  can  be  achieved. 
This  technique  is  similar  to  an  aircraft  flying  an  arc 
about  a  TACAN  station.  Thus,  the  initial  and  final  AB 
were  chosen  to  be  75  and  105  degrees  (90  ±  15) . 

Wild  Weasel  Turn  Direction  Algorithm 

The  first  pilot  input  to  the  WW  flight  scenario 
is  to  determine  the  initial  turn  direction  required  to 
achieve  the  desired  AB  of  75  degrees.  As  was  mentioned 
earlier,  time  is  critical  during  the  ranging  routine.  The 
WW's  self -protect ion  and  evasive  tactics  are  limited 
during  this  time  period,  so,  the  turn  to  75  degrees  should 
be  in  the  direction  which  takes  the  shortest  time.  In  the 
aircraft,  the  decision  is  simple.  The  pilot  looks  at  his 
RHAW  display  and  determines  the  KB  to  the  site.  He  turns 
either  left  or  right  depending  on  the  value  of  KB.  If  the 
site  is  at  the  pilot’s  10  O'clock  position  (as  in  Figure  6), 
the  turn  will  be  to  the  right. 

Converting  from  the  pilot  "seeing"  the  RB  to  the 
model  computing  the  direction  can  be  accomplished  once  the 
B  and  AB  have  been  calculated.  The  initial  decision  is  to 
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compare  the  site's  x  coordinate,  Xg,  with  the  aircraft's 
x  coordinate,  xft.  If  (Xg  >  xft)  the  bearing  to  the  site 
(in  degrees)  will  fall  in  the  range  (000<  B<  180).  Next, 
the  model  analyzes  which  side  of  the  site  the  aircraft 
will  pass  or  has  passed,  if  the  aircraft  maintains  a  con¬ 
stant  H.  From  Figure  6,  the  aircraft  will  pass  to  the 
right  of  the  site  if  it  continues  on  its  present  heading. 
Converting  to  the  model's  logic,  if  [B  <  H  <(B  +  180)], 
the  aircraft  will  pass  to  the  right  of  the  site.  Again, 
from  Figure  6,  the  aircraft's  H  is  090  degrees,  the  B  to 
the  site  is  45  degrees.  The  heading  falls  in  the  range  of 
(45  <  90  <  225)  and  satisfies  the  requirement  for  passing  to 
the  right  of  the  site. 

After  determining  the  site's  position  in  relation 
to  the  aircraft  and  whether  the  aircraft  will  pass  or  has 
passed  to  either  the  right  or  left  of  the  site,  the  final 
comparison  checks  the  magnitude  of  AB.  If  (AB<  75)  degrees 
the  turn  will  be  made  to  increase  AB  to  75  degrees.  If 
(AB  >75  degrees)  the  turn  will  be  in  a  direction  to 
decrease  AB  to  75  degrees.  This  completes  the  logic 
required  to  turn  the  aircraft  in  the  shortest  direction  to 
initiate  the  ranging  phase. 

At  the  completion  of  ranging,  the  aircraft  turns 
again,  this  time  toward  the  site  and  prepares  for  the 
attack  phase.  This  second  turn  direction  is  dependent 
only  on  the  site's  position  in  relation  to  the  aircraft 


prior  to  the  initial  ranging  turn.  Table  II  summarizes 
all  positions  and  relative  bearing  combinations  prior  to 
the  initial  turn  with  the  decision  blocks  indicating  the 
direction  of  the  initial  turn  to  start  the  ranging  phase 
and  the  second  turn  to  initiate  the  attack  phase.  For 
example,  for  the  aircraft  in  Figure  6,  it  is  to  the  right 
of  the  site  and  AB<  75  degrees.  To  start  the  ranging,  the 
aircraft  would  turn  right  and  to  start  the  attack  phase 
the  second  turn  would  be  to  the  left. 

TABLE  II 


WILD 

WEASEL  TURN 

LOGIC 

Side  of 
Closest 
Approach 
to  Site 

<  75 

>  75 

Ranging 

Turn 

Attack 

Turn 

Ranging 

Turn 

Attack 

Turn 

Left 

Left 

Right 

Right 

Right 

Right 

Right 

Left 

Left 

Left 

Addressing  the  requirement  for  turning  in  the 
shortest  direction,  a  question  arises  as  to  how  much  time 
is  saved  in  a  WW  mission  scenario  by  turning  in  the 
shortest  direction.  Figure  8  depicts  to  scale  the  situa¬ 
tion  where  a  WW  detects  a  threat  at  a  point,  turns  in  the 
shortest  direction  to  initiate  ranging,  then  turns  in  the 
shortest  direction  to  attack  the  site.  Figure  9  shows 
the  aircraft  under  the  same  initial  AB  and  H  conditions. 
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but  which  maneuvers  in  the  incorrect  direction  for  both 
turns.  The  increased  time  required  by  the  "wrong  turn" 
aircraft  in  completing  the  profile  can  be  translated  into 
fewer  sites  attacked  during  a  mission  or  increased  exposure 
to  the  enemy  threats.  This  increased  exposure  translates 
into  an  increased  probability  of  either  being  detected  by 
the  enemy  sites  or  by  being  destroyed  after  being  detected. 
Thus  maneuvering  in  the  shortest  direction  to  achieve  the 
required  AB  increases  the  number  of  sites  attacked  per 
mission  and  may  decrease  and  probability  of  the  WW  being 
destroyed. 

Attack  Phase  of  WW  Profile 

The  final  and  most  important  phase  of  the  WW’s 
mission  profile  is  the  actual  attack  of  the  site.  After 
completing  the  ranging  phase  and  with  an  accurate  estimate 
of  the  site's  location  the  WW  turns  back  toward  the  threat. 
In  this  thesis  the  aircraft  stops  turning  when  its  flight 
vector  aligns  with  the  site's  location  and  the  WW  is  bore- 
sighted  on  the  threat  (AB  *  0  degrees) .  Although  not 
required  by  all  types  of  ARMs ,  boresighting  on  the  threat 
increases  the  ARM's  delivery  accuracy.  Any  launch  azimuth 
off  of  boresight  results  in  a  decreased  delivery  accuracy. 

The  WW  carries  two  types  of  ARMs:  a  type  B  desig¬ 
nated  for  a  long  range  launch  and  an  older,  shorter  range 
missile  designated  type  A.  For  a  given  launch  speed  and 


33 


altitude  each  ARM  has  an  associated  envelope  that  denotes 
the  minimum  and  maximum  firing  distances,  representative 


flyout  velocities,  and  an  associated  P^.  For  the  WW  alti¬ 
tude  used  in  this  model,  these  hypothesized  parameters  are 
shown  in  Table  III.  Given  the  increased  for  ARM  type  B 
the  decision  logic  for  selecting  a  weapon  will  be  as  fol¬ 
lows.  If  any  type  Bs  are  left,  the  WW  crew  selects  a 
type  B  weapon. 


TABLE  III 

ANTI -RADIATION  MISSILE  (ARM)  CHARACTERISTICS 


Average  Velocity  (m/sec) 
Probability  of  Kill  (P^) 
Range  (m)  (Min/Max) 


8000/15,000 


8000/25,000 


This  thesis  assumes  weapons  delivery  accuracy  is  indepen 
dent  of  launch  range  when  launch  occurs  in  the  launch 
window.  After  boresighting  on  the  threat  the  WW's  next 
action  depends  on  the  slant  range  to  the  site.  Three 
possible  conditions  exist: 

1.  The  SR  falls  in  the  launch  envelope, 

2.  The  SR  is  greater  than  the  maximum  launch 
range  (Rj^)  ,  or 
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3.  The  SR  is  less  than  the  minimum  launch 
range  <Rmin) . 

If  the  SR  satisfies  the  first  case,  the  launch  point  will 

occur  somewhere  between  the  current  SR  and  the  R  .  . 

ram 

The  model  assumes  that  any  point  between  SR  and  R  .  is 

nun 

equally  likely  and  determines  a  launch  range,  Rlau»  based 

on  a  uniform  probability  distribution  between  SR  and  R  . 

mm 

For  case  2,  the  aircraft  must  delay  launch  until 

it  enters  the  launch  envelope.  Here,  the  launch  point  is 

assumed  to  be  uniform  between  R  _  and  R  .  .  For  both 

max  mm 

cases  1  and  2  the  WW  proceeds  toward  the  site  until  its 
SR  reaches  launches  the  ARM,  allows  the  weapon  to 

clear  the  aircraft  (assumed  to  be  launch  time  plus  five 
seconds) ,  and  begins  searching  for  the  next  threat. 

Case  3  requires  the  WW  to  maneuver  away  from  the 
site.  Using  the  turn  algorithm,  the  WW  aircraft  turns 
away  from  the  site  until  AB  *  180  degrees,  rolls  out  of 
the  turn,  and  proceeds  outbound  until  the  SR  to  the  site 
increases  to  18,000  meters,  reverses  its  turn  and  bore- 
sights  on  the  site.  This  insures  the  second  attack  will 
fall  in  case  2  and  a  launch  will  culminate  the  attack. 

The  case  3  profile  typifies  action  required  by  the  WW  in 
a  threat-rich  environment  given  the  launch  envelope  of 
the  ARMs. 


WW  Defense  Suppression 
Mission  Summary 

In  summary,  three  key  elements  comprise  the  WW's 
defense  suppression  system  structure: 

1.  The  heading  reference  system  (HRS)  and  its 
associated  logic, 

2.  The  aircraft  maneuvering  logic,  and 

3.  The  Defense  Suppression  Mission  Phases. 

The  HRS  can  be  compared  to  an  aircraft's  hori¬ 
zontal  situation  indicator.  It  allows  the  model  to 
describe  a  threat  radar's  position  in  relation  to  an  air¬ 
craft's  current  position  and  heading.  The  first  two  out¬ 
puts  calculated  by  the  HRS  are  the  aircraft's  heading  and 
the  magnetic  bearing  to  a  threat  radar  site.  These  two 
elements,  in  turn,  define  a  third  quantity,  the  relative 
bearing  of  the  threat  radar  to  the  aircraft's  current  head¬ 
ing.  Knowing  these  three  elements  allows  the  model  to 
solve  the  maneuvering  logic  to  accomplish  the  mission 
scenario. 

Four  key  phases  of  the  defense  suppression  mission 
scenario  are  the  following: 

1 .  Search , 

2 .  Detection , 

3 .  Ranging ,  and 


4.  Attack. 


The  search  phase  occurs  when  the  WW  has  weapons 
remaining  and  is  attempting  to  locate  a  target  to  destroy. 
At  detection,  the  aircraft's  RHAW  system  alerts  the  crew 
to  a  potential  target.  During  ranging,  the  WW's  onboard 
processor  refines  the  site's  location  estimate  to  within 
a  tolerance  required  by  the  aircraft's  weapons.  The  attack 
phase  culminates  the  mission  scenario.  The  aircraft 
selects  an  ARM  and  launches  it  against  the  site.  The  mis¬ 
sion  scenario  continues  until  the  WW  exhausts  its  ARMS. 

Although  important  to  the  overall  success  of  the 
mission,  the  WW  egress  from  the  target  area  was  not 
addressed  in  this  chapter.  The  model  will  handle  egress 
by  assigning  the  x  and  y  coordinates  of  (0,  25000  m)  as 
the  site  coordinates  of  the  next  target  when  all  ARMs 
are  depleted.  The  WW  turns  toward  this  point,  boresights 
on  it  and  departs  the  enemy  portion  of  the  FEBA. 

Strike  Aircraft  System  Structure 

The  strike  aircraft  flight  forms  with  the  WW  over 
friendly  territory  well  before  the  FEBA.  The  planned 
ingress  corridor  keeps  the  strike  force  as  far  away  as 
possible  from  the  enemy's  main  lines  of  communication  (LOC) 
(Ref  9:83).  The  enemy  uses  these  LOCs  to  transfer  troops 
and  equipment  forward  to  and  away  from  the  FEBA,  and  thus, 
concentrates  his  defenses  near  them.  The  mission's  target 
is  located  in  an  area  95  km  behind  the  FEBA.  The  aircraft 


penetrate  the  FEBA  near  the  target's  latitude  at  an  alti¬ 
tude  of  200  m.  The  speed  to  the  target  is  a  constant 
480  kts  (247  m/sec) .  No  aircraft  turning  in  response  to 
threats  is  modeled  and  the  ingress  heading  is  kept  constant 
at  090  degrees. 

Defensive  System  Structure 

The  final  system  element  is  the  defensive  array 
the  WW  and  strike  aircraft  will  attempt  to  penetrate.  The 
major  areas  addressed  in  this  section  are  the  following: 

1.  Radar  system  characteristics  and  limitations, 

2.  SAM  probability  of  kill  calculations, 

3.  Intercept  geometry  calculations, 

4.  Specific  differences  with  previous  modeling 

efforts, 

5.  AAA  probability  of  kill  calculations,  and 

6.  Command  and  control. 

Radar  System  Characteristics 
and  Limitations 

The  radar  systems  composing  the  defensive  network 
can  be  characterized  by  certain  system  parameters  listed 
in  Table  IV  (Ref  9:18) .  The  Pfc,  Gfc,  and  RF  are  the  same 
parameters  defined  in  the  receiver  sensitivity  section  of 
this  thesis.  The  maximum  range  and  minimum  altitude 
represent  restrictions  on  either  the  radar  itself  or  the 
missile  associated  with  the  radar.  The  acquisition  and 
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TABLE  IV 


tracking  time  is  the  time  required  by  the  system  operators 
after  detecting  an  aircraft  to  evaluate  the  flight  path 
and  to  determine  if  the  probability  of  destroying  the  tar¬ 
get  warrants  further  tracking  by  the  site.  At  its  minimum 
limit,  it  represents  the  minimum  time  after  detecting  a 
threat  for  a  site  to  schedule  a  missile  launch.  The  missile 
velocity  is  the  average  velocity  of  the  missile  from  launch 
to  detonation.  In  this  thesis  missile  acceleration  will 
not  be  modeled.  The  confounding  delay  is  the  minimum 
time  required  by  the  site  to  prepare  the  system  for  the 
next  engagement.  Typical  examples  of  confounding  delays 
could  be  the  time  after  missile  detonation  required  by 
the  site  before  searching  for  a  new  target  or  the  time  for 
a  site  to  access  a  target  disappearing  from  its  radar 
scope  (an  aircraft  going  out  of  range  or  being  destroyed 
by  another  site) . 

For  this  thesis'  model,  the  defensive  radar  site 
must  be  capable  of  analyzing  aircraft  engagements  from  any 
attack  azimuth  and  heading.  Additionally,  the  aircraft 
can  employ  either  jamming  (wet)  or  nonjamming  (dry)  tactics 
(WW  aircraft  do  not  operate  ECM  equipment  during  a  defense 
suppression  mission  profile) ^ 

The  methodology  required  for  determining  the  proba¬ 
bility  of  destroying  the  attacking  aircraft  will  be  devel¬ 
oped  in  the  following  paragraphs. 


Radar  system  parameters  such  as  the  pulse  repeti¬ 
tion  interval  (PRI)  determine  the  maximum  range  most  radars 
can  detect  and  track  a  target.  The  accuracy  achieved  by 
these  radars  and  their  associated  missiles  has  led  the 
attacking  aircraft  force  to  employ  low  altitude  penetra¬ 
tion  as  one  means  of  increasing  its  probability  of  reach¬ 
ing  the  target  area.  The  low  altitude  penetration  compli¬ 
cates  the  radar  tracking  problem  because  it  introduces  the 
multipath  signal  problem.  Below  a  certain  altitude,  the 
reflected  radar  return  from  a  target  interacts  with  a 
mirror  image  signal  reflected  off  the  earth's  surface  at 
the  radar  processor  (Ref  15:172) .  Because  the  target 
return  interferes  with  and  cannot  be  distinguished  from 
the  mirror  image,  accurate  tracking  cannot  be  accomplished. 
The  signal  can  induce  wild  fluctuations  in  the  radar's 
range  gate  and  causes  the  return  to  "break  lock."  At 
some  elevation  angle,  referred  to  as  the  multipath  angle, 
the  difference  in  the  arrival  time  between  the  signal  and 
its  mirror  image  is  great  enough  that  the  two  signals 
can  be  distinguished  and  accurate  tracking  of  the  actual 
return  becomes  possible.  For  the  radars  in  this  thesis, 
an  average  multipath  angle  of  .25  degrees  will  be  used. 

This  represents  an  average  of  the  values  used  by  Leek  and 
Schmidt  (Ref  9:15).  Thus,  in  addition  to  being  inside 
the  maximum  detection  range  of  a  radar,  the  aircraft 
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must  be  at  an  elevation  angle,  a,  above  the  horizon  such 
that  a  >.25  degrees,  where  a  is  defined  as  follows: 


*  1  , ALT <  c \ 

a  *  sin  (“§5“)  (15) 

where  ALT  *  altitude  of  the  aircraft,  and 

SR  *  slant  range  from  the  radar  to  the  aircraft. 

After  detecting  the  aircraft,  the  radar  site 
attempts  to  track  the  target  and  launch  its  missile  in 
such  a  manner  as  to  achieve  the  highest  probability  of 
kill  (P^) .  In  a  similar  manner  to  Leek  and  Schmidt ,  this 
thesis  evaluates  the  P^  in  terms  of  the  missile's  lethal 
radius  (LR)  and  the  circular  error  probable  (CEP)  (Ref 
9:19-22).  LR  is  the  distance  from  the  SAM's  detonation 
point  where  as  many  aircraft  survive  as  are  killed  by  the 
detonation  beyond  it.  CEP  is  the  error  associated  with 
the  distance  measured  between  the  desired  and  actual  point 
of  impact  (Ref  2:44).  More  precisely,  CEP  or  spherical 
error  probable,  is  a  sphere  around  the  target  aircraft 
within  which  50  percent  of  the  missiles  fired  under  a 
given  set  of  conditions  will  detonate  (Ref  9:19) .  Thus, 

Pk,  evaluated  in  terms  of  CEP  and  LR  reduces  to  the  follow¬ 
ing  equation: 

Pk  -  l-(.5)  (LR/CEP)2  (Ref  9:22)  (16) 
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For  a  specific  SAM  system,  LR  is  a  constant.  By 
decreasing  CEP,  the  missile's  P^  increases;  for  the  high¬ 
est  P^,  the  missile  must  deliver  its  missiles  with  the 
smallest  CEP. 

CEP  calculations  depend  on  whether  the  target  air¬ 
craft  is  either  wet  or  dry — either  jamming  or  not  jamming. 
Leek  and  Schmidt  evaluated  the  wet  case  where  CEP  is  deter¬ 
mined  as  follows: 

CEP  «  yA(J/S)Ri2+B(J/S)+C  (17) 

where  A,  B,  and  C  *  constants  dependent  on  the  type  of  SAM, 

R.  *  range  from  the  launch  point  to  the 
1  target,  and 

J/S  *  the  jamming  to  signal  ratio. 

The  range  to  intercept,  K^,  will  be  evaluated  later 
in  this  section  for  all  intercept  cases.  (J/S)  is  depen¬ 
dent  on  the  range  from  the  radar  to  the  target  aircraft, 
the  effective  radiated  power  (ERP)  of  the  aircraft's 
jammer,  the  technical  parameters  of  the  radar,  and  aRCS,  the 
radar  cross-section  (RCS)  of  the  target  aircraft.  RCS  is 
dependent  on  the  intercept  geometry  between  the  aircraft 
and  radar  and  represents  the  amount  of  reflected  signal 
energy  reradiatad  back  toward  the  radar.  It  is  usually 
expressed  in  decibels  referenced  to  square  meters  of 
reflected  signal  strength.  Table  V  lists  the  RCS  in 
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TABLE  V 


RADAR  CROSS-SECTION  OF  FIGHTER 
AGAINST  DEFENSIVE  SYSTEMS 


Aspect 

(Degrees) 

Radar  Cross-Section 

SAM-A 

(in  dBm^) 

SAM-B, 

SAM-C 

SAM-D, 

AAA 

0. 

10.75 

3.60 

8.20 

5. 

4.75 

.40 

6.70 

10. 

4.65 

2.55 

1.90 

15. 

-2.33 

3.03 

3.48 

20. 

4.80 

1.70 

.85 

25. 

2.75 

2.50 

3.95 

30. 

1.20 

10.65 

5.20 

35. 

-1.65  ) 

7.95 

-1.20 

40. 

2.15 

3.45 

1.70 

45. 

4.70 

3.70 

6.35 

50. 

1.75 

-.95 

3.10 

55. 

1.10 

-.65 

1.00 

60. 

-4.00 

-.95 

1.90 

65. 

.43 

.05 

.25 

70. 

9.23 

8.45 

8.19 

75. 

13.73 

14.55 

13.43 

80. 

13.98 

16.35 

16.70 

85. 

16.38 

16.00 

16.08 

90. 

25.85 

24.98 

24.38 

95. 

20.88 

19.95 

19.23 

100. 

17.00 

15.75 

16.58 

105. 

7.03 

9.20 

5.83 

110. 

9.80 

8.73 

9.50 

115. 

1.75 

3.65 

6.20 

120. 

-.75 

2.33 

7.85 

125. 

-.23 

3.13 

4.28 

130. 

-2.33 

3.03 

3.48 

135. 

.08 

-1.08 

4.35 

140. 

-1.28 

-2.60 

3.90 

145. 

-1.90 

.50 

6.00 

150. 

2.78 

.55 

5.23 

155. 

6.55 

.43 

6.93 

160. 

.50 

.58 

2.95 

165. 

.55 

6.38 

4.73 

170. 

4.15 

6.53 

9.93 

175. 

4.63 

8.85 

13.50 

180. 

.98 

9.48 

15.05 

2 

dB  (meters  )  for  the  aircraft  in  this  thesis  and  listed 
as  an  intercept  angle. 

J/S  can  be  determined  by  the  following  formula: 


am  -  f >  ’  ^ 

Pt°t  JRCS 


where  all  terms  are  as  previously  defined  (Ref  9:101-102) 
Converting  all  terms  in  equation  (18)  to  dB  equivalents 
yields: 

( J/S) ds  =  2(R)dfi  +  ( 4tt ) dB  +  ERPdB  -  P 

CLB 

'  ‘  0ECSdB  (19> 

For  a  particular  aircraft/SAM  engagement,  the  terms  4it  , 
ERP,  Pfc,  and  Gt  are  constant  and  the  J/S  ratio  becomes: 

(J/S)^  -  2(R)dB  -  o  +K  (20) 


where  K  *  4ir  +  ERP  ”  ”  Tt' 


The  J/S  required  for  equation  (17)  is  obtained  by  convert¬ 
ing  the  dB  equivalent  back  to  the  actual  value,  or: 


(J/S)  -  10 


(J/SJ^/IO 


From  equation  (17) ,  a  smaller  value  of  (J/S) 
results  in  a  smaller  CEP.  The  radar  attempts  to  control 
(J/S)  by  insuring  the  missile  impact  occurs  at  the  minimum 


possible  range  from  the  site  and  the  largest  value  of  RCS. 
Both  values  can  be  achieved  simultaneously  when  the  missile 
impacts  near  the  aircraft  as  the  aircraft  nears  its  closest 


approach  point,  R  ,  to  the  site  or  the  aircraft  comes  abeam 

c 

the  site.  Here,  the  intercept  angle  is  90  degrees  and 
from  the  values  of  RCS  in  Table  V,  RCS  is  a  maximum.  The 
launch  rule  for  the  radar  sites  in  this  thesis  will  be  to 
attempt  a  launch  when  the  intercept  with  the  target  occurs 
at  the  closest  approach  point  to  the  site.  This  firing 
methodology  is  similar  to  the  constant  bearing  intercept 
discussed  by  Breuer  (Ref  3:166-173) .  A  constant  bearing 
intercept  provides  the  shortest  duration  trajectory  for  a 
constant  velocity  target. 

For  a  dry  aircraft,  CEP  can  be  evaluated  as 

follows: 


where  D,  E,  and  F  are  constants  dependent  on  the  type 
of  SAM  and  all  other  terms  as  previously  defined. 

The  aRcg  required  in  equation  (22)  is  the  one 
obtained  by  converting  the  dB  equivalent  to  the  actual 
values ,  or : 


a  -  10Odfl/1° 


(23) 


46 


Again,  the  lowest  CEP  and  the  highest  result  when  the 
intercept  range  is  a  minimum  and  the  RCS  is  a  maximum. 

For  both  wet  and  dry  aircraft,  the  same  launch  decision 
rule  will  be  applied:  attempt  to  launch  so  that  intercept 
occurs  at  90  degrees  angle  between  the  aircraft's  flight 
vector  and  the  site's  bearing  to  the  aircraft  (point  C  in 
Figure  11) . 

Intercept  Geometry 

At  the  completion  of  acquisition  and  tracking  the 
model  evaluates  the  specific  encounter  conditions  between 
the  threat  radar  and  an  aircraft.  During  acquisition  and 
tracking  the  site  determines  the  aircraft's  heading,  H, 
velocity,  V  ,  and  both  its  slant  range,  SR,  and  bearing 

a 

from  the  site.  For  each  encounter  two  major  areas  dictate 
the  intercept  computations.  Figures  10  and  12  depict 
these  conditions.  Based  on  its  position,  velocity,  and 
heading  the  aircraft  is  located  either  prior  to  or  past 
the  closest  approach  point  to  the  site,  C.  Within  the 
first  case  (prior  to  the  closest  approach  point)  two  pos¬ 
sible  conditions  exist:  the  missile  if  fired  immediately  at 
the  end  of  acquisition  and  tracking  will  intercept  the  air¬ 
craft  at  or  prior  to  the  closest  approach  point,  or  the 
missile  will  intercept  the  aircraft  past  the  closest 
approach  point.  Referring  to  Figure  10,  an  aircraft  at 
point  A  is  headed  in  a  direction  such  that  its  closest 
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Fig.  10.  Aircraft  Prior  to  Closest 
Approach  Point  (3-Dimension) 


Fig.  11.  Aircraft  Prior  to  Closest 
Approach  Point  (2-Dimension) 


Fig.  12.  Aircraft  Past  Closest  Approach 
Point  (3-Dimension) 


approach  to  the  site  S  occurs  at  point  C.  The  slant  range 
from  the  site  to  the  aircraft,  SR,  is  known  (see  equation 
(14)).  The  SR  is  projected  down  to  the  x-y  plane  as 
follows  (see  Figure  11) : 


SR  cos  a 


where  SR  *  the  projection  of  SR  in  the  x-y  plane,  and 
xy 

all  other  terms  as  previously  defined. 


Using  the  bearing  from  the  site  to  the  aircraft  and  the 
aircraft's  H,  the  angle  6  is  calculated.  The  time  for  the 
aircraft  to  reach  point  C,  tAC,  is  calculated  as  follows: 


SR _ cos  9 


where  AC  ■  the  length  of  the  leg  from  A  to  C  or  SR 
cos  6  , 

V  *  the  aircraft's  velocity,  and 

cl 

all  other  terms  as  previously  defined. 


The  time  required  for  the  missile  to  get  from  S  to  C,  tMC, 
is  also  known: 


Sic  “  v 
m 


.J(SR  sin0 )  *  +  ALT' 


where  SC  »  distance  from  point  S  to  point  C, 

Vm  *  average  missile  flyout  velocity,  and 
all  other  terms  as  previously  defined. 


Note:  For  a  strike  aircraft  with  a  constant  H,  the  R^ 
depends  only  on  the  difference  between  the  site's  and  the 
aircraft's  y  coordinates,  Ay,  and  the  aircraft's  altitude, 
ALT: 


SC 


+  (Ay)2 


(27) 


where  all  terms  are  as  previously  defined. 

The  model  now  compares  tAC  and  tMC-  If  (tAC >  tM(J 
the  missile,  if  fired  now,  would  arrive  at  point  C  before 
the  aircraft.  Since  the  intercept  must  occur  at  C  to 
achieve  the  lowest  CEP,  and  highest  P^,  a  delay  of 
(tj\c  “  SlC^  t*me  “nits  is  programmed  into  the  missile's 
firing  schedule,  resulting  in  both  missile  and  aircraft 
arriving  at  point  C  simultaneously. 

If  (tjjj,  >  tAC)  the  missile,  if  fired  now,  would 
arrive  at  C  after  the  aircraft.  In  this  case,  the  inter¬ 
cept  will  occur  past  the  optimum  90  degree  intercept  point. 
The  new  intercept  point,  0,  can  be  calculated  as  follows: 


®  -  <tMC-tAC>Va 


(28) 


where  CD  «  the  distance  past  point  C  the  intercept  will 
occur,  and 

all  other  terms  as  previously  defined. 

The  actual  range  to  intercept,  can  be  calculated 


as  follows: 


R.  «  WCD  +  (SR  sin  0)  +  ALT 

1  1  Xjf 

Ri 

Sjd  “  v" 

m 

where  t^  *»  time  of  the  missile  to  flyout  to  point  D,  and 
all  the  terms  are  as  previously  defined. 

Instead  of  the  RCS  angle  being  an  optimum  90  degrees 
it  will  be  a  new  angle ,  if) ,  where  $  can  be  calculated  as 
follows  (see  Figure  10  and  11) : 

5  *  sin”1  (31) 

Ri 


(29) 

(30) 


Y 


cos 


-1 


CD 

cos 


(32) 


$  *  180  -  y  (33) 

where  all  terms  as  previously  defined. 

Thus*  the  can  be  calculated  with  CEP  determined  by  the 
RCS  for  an  aspect  angle  $  and  the  intercept  range*  R^. 

At  this  time  the  site,  knowing  the  expected  P^* 
would  evaluate  its  decision  to  fire.  If  the  P^  is  above 
some  threshold  value*  the  decision  would  be  made  to  fire; 
if  below,  the  site  would  be  released  from  this  threat  and 
resume  searching  for  another  target  aircraft.  The 
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threshold  value  depends  on  other  criteria  such  as  the  num¬ 
ber  of  missiles  available  for  firing  and  the  priority  placed 
on  destroying  the  strike  force/WW  aircraft.  For  FEBA 
engagements  early  in  a  conflict  with  high  stockpile  levels 
of  missiles ,  this  threshold  value  will  be  low.  For  this 
thesis,  values  above  .05  will  be  considered  adequate 
to  warrant  the  site  firing  at  the  aircraft.  Leek  and 
Schmidt  used  an  even  lower  threshold  value  of  .02  in  the 
firing  logic  of  their  thesis  (Ref  9:50).  If  the  is 
above  the  threshold  value,  the  missile  is  launched  immedi¬ 
ately  with  intercept  occurring  t^  time  units  later. 

The  second  major  category  of  target  is  one  already 
past  the  closest  approach  point  at  the  completion  of  track¬ 
ing  and  acquisition.  Again,  the  site  has  determined  the 
SR  and  bearing  from  the  site  plus  the  aircraft's  H  and  V  . 

a 

(See  Figure  12  and  13,  page  49.)  The  range  to  missile 
intercept  can  be  calculated  as  follows.  The  intercept  will 
occur  at  a  point,  B.  The  time  for  the  aircraft  to  get 
from  its  present  position,  A,  to  the  intercept  point, 
tAB'  wil1  be  the  same  as  the  time  required  for  the  missile 
to  go  from  S  to  B,  t^,  or: 


AB 


tMB 


(34) 

(35) 
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and 


(36) 


V 

AS  -  -  VrR. 

in 

where  Vr  *  the  velocity  ratio,  Va/Vm,  and  all  other 
terms  are  as  previously  defined. 


For  small  angles,  S  ,  * 

follows  (see  Figure  12) : 


R.  ,  where  S 
xxy 


is  defined  as 


,  -1  .ALT. 
sin  (*r7_) 


L  *  R.  cos  6 

v  1 


where  all  terms  are  as  previously  defined. 


(37) 

(38) 


The  distance,  R^  cam  be  determined  using  the  Law  of 
Cosines: 


R^2  ■  (SR  cos  a)2  +  (AB)  2 
-  2 (SR  cos  a)  (AB)  cos  0 


(39) 


where  all  terms  are  as  previously  defined. 

The  angle  0  is  determined  because  the  aircraft's  heading 
and  its  bearing  from  the  site  are  both  known.  Substituting 
from  equation  (36)  for  AB,  equation  (39)  becomes: 

R,2  ■  (SR  cos  a)2  +  (vrRi^2 

-  2  (SR  cos  a)  (VyR^  cos  0  (40) 

where  all  terms  are  as  previously  defined. 
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Rearranging  terms,  the  equation  becomes  a  quadratic  in 
2 

R.  .  Equation  (40)  is  solved  by  the  quadratic  formula 
with: 


b  =  2  (SR  cos  a)  (Vr)cos  0 
2 

c  *  -  (SR  cos  a ) 

-  _  -b  ±\T b2  -  4ac 

Ri  *  25  (41> 

The  negative  of  the  term  under  the  radical  is 
neglected  since  R^  is  a  positive  number.  The  RCS  angle 
at  intercept  can  be  determined  since  the  bearing  to  the 
target ,  0 ,  and  the  target ' s  heading  is  known  at  intercept . 


Rc  - 

SRxy  cos  e 

(42) 

8  - 

-1  Ri 
cos  a(5±) 

(43) 

c 

Y  +6  - 

90° 

(44) 

$  m 

RCS 

180  -  y 

(45) 

(b  m 

yRCS 

90  +  0 

(46) 

Again,  knowing  R^  and  the  RCS  aspect  angle  $ ,  the  CEP  and 
P^  can  be  evaluated  at  intercept.  Similar  to  the  cases 
prior  to  the  closest  approach  point  P^  is  above  the  .05 
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threshold  value,  the  missile  is  fired.  Because  the  inter¬ 
cept  is  beyond  the  closest  approach  point,  the  missile  is 
fired  immediately  and  intercept  occurs  t^  time  units 
later. 

An  additional  consideration  necessary  before 
leaving  the  SAM  attack  geometry  calculations  is  the  neces¬ 
sity  of  the  calculations  to  handle  a  maneuvering  aircraft. 
The  initial  evaluation  at  the  end  of  acquisition  and 
tracking  was  made  considering  the  aircraft's  heading  to  be 
constant.  Since  the  WW  aircraft  can  maneuver,  this  may  not 
be  true.  To  account  for  a  turning  possibility  this  thesis 
incorporates  the  following  modification.  At  the  scheduled 
missile  firing  time,  a  new  is  calculated  (see  equation 

(29)  or  (41)).  If  the  R.  is  the  same  as  the  initial 

■^new 

R ,  calculated  at  the  end  of  acquisition  and  tracking, 
the  aircraft  has  not  turned  and  the  intercept  will  be  the 
same.  To  preclude  small  variations  in  the  Ri  calculations, 
a  .1  tolerance  is  allowed,  or,  intercept  will  be  considered 
the  same  if: 


(.9  R .  <  R .  <  1.1  R.)  (47) 

1  new 


If  (R.  <  R.) ,  the  aircraft  must  have  turned  toward  the 

1new  1 

site.  A  new  $  and  are  calculated  and  a  new  launch 

time  computed  based  on  the  R.  .  A  new  launch  time  is 

xnew 

computed  using  the  above  calculations  (equations  (25) -(26)) 
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and  a  new  launch  time  is  schedule.  If  (R.  >R.),  the 

xnew  1 

new  P^  is  evaluated.  If  the  is  above  the  .05  threshold 
value,  the  site  immediately  fires  the  missile  and  computes 
a  scheduled  impact  time. 

At  the  scheduled  impact  time  the  R^  is  again  com¬ 
pared  to  the  value  computed  at  the  missile  launch  time. 

If  both  the  R^  values  are  the  same,  the  P^  can  be  evaluated 
If  the  new  R^  is  less  than  the  launch  R^,  the  aircraft  has 
turned  towards  the  site.  The  and  <j>  are  evaluated  at 
the  new  R^.  If  the  new  R^  is  greater  than  the  launch  R^ 
a  new  impact  time  is  scheduled.  This  new  impact  time  is 


R.  R. 

,  xnow  v  ,  1  launch  . 

'D  *  V  y  ' 


1  launch 


(48) 


m 


or 


new 


(49) 


m 


where  all  terms  are  as  previously  defined. 

At  the  new  impact  time,  the  same  calculations  are 
made  and  the  reevaluated.  In  this  manner,  the  flight 
time  of  the  missile  can  be  approximated  without  requiring 
a  separate  missile  flyout  time  routine  as  part  of  the 
thesis. 
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A 


Differences  with  Previous 


The  SAM  engagement  portion  of  this  thesis  differs 
from  the  approach  used  by  Leek  and  Schmidt.  In  their  model 
Leek  and  Schmidt  assumed  that  if  the  attack  aircraft  main¬ 
tained  a  20  dB  J/S  ratio  at  the  radar  the  site  was  denied 
the  ability  to  track  in  both  range  and  azimuth.  Using  this 
criteria,  kill  zones  in  which  the  J/S  dropped  below  20  dB 
for  the  particular  SAM  were  established.  If  the  aircraft 
passed  through  this  kill  zone  at  the  scheduled  impact  time 
the  was  calculated.  If  the  aircraft  intercept  occurred 
outside  the  zone,  the  P^  was  set  equal  to  zero.  To  this 
end  the  Leek  and  Schmidt  model  attempted  to  control  the 
SAM  launches  such  that,  if  possible,  launches  were  delayed 
so  that  intercept  occurred  when  the  aircraft  entered  a 
kill  zone.  The  kill  zone  criteria  resulted  in  a  narrow 
azimuth  window  near  the  90  degree  relative  bearing  to  the 
site  (Ref  9:24-31). 

The  main  criticism  is  that  this  approach  does  not 
consider  even  the  most  rudimentary  electronic  counter 
countermeasures  (ECCM)  techniques.  Most  modem  SAMs  employ 
some  type  of  ECCM  such  as  .frequency  discrimination,  side- 
lobe  blanking,  or  polarization  mismatch  to  defeat  ECM  tech¬ 
niques.  Although  not  explicitly  stated,  attack  aircraft 
in  the  Leek  and  Schmidt  model  used  noise  jamming.  This 
type  of  jamming  normally  denies  range  information  to  a  site 
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A 


although  angle  information  can  be  obtained  since  the  out¬ 
put  from  the  noise  jammer  results  in  a  strobe  on  the  radar 
scope  of  the  site  eminating  from  the  azimuth  of  the  target 
aircraft.  In  this  case,  a  network  of  sites  such  as  the 
model  employs  could  triangulate  and  to  determine  the  air¬ 
craft's  approximate  location  {Ref  15:547-550).  For  this 
reason  kill  zones  are  not  defined.  Instead,  a  necessary 
assumption  in  the  analytical  portion  of  the  model  is  that 
the  sites  use  some  type  of  ECCM  techniques  and  the  air¬ 
craft's  PR  is  based  on  the  calculated  value  of  CEP,  R^, 
and  oRCS  determined  in  the  above  equations,  not  whether 
the  aircraft  at  missile  impact  falls  in  a  kill  zone.  For 
large  values  of  Ri  and  high  CEPs  an  appropriate  decrease 
in  PR  will  be  calculated  as  per  equations  (17)  and  (22) . 


AAA  Probability  of  Kill 

The  AAA  PR  calculations  follow  closely  the  work  of 
Leek  and  Schmidt  (Ref  9:35-40).  The  short  maximum 
engagement  range  of  the  AAA  (2990  m,  s*e  Table  IV)  means 
the  aircraft  must  pass  near  the  site  before  being  attacked. 
The  Pfc  calculations  are  based  on  the  vulnerable  surface 
area  of  the  aircraft  as  a  percentage  of  total  surface  area 
as  viewed  from  a  specific  aspect  angle,  the  dispersion  of 
the  AAA  rounds  around  the  actual  aim  point,  the  time  of 
flight  (TOF)  of  each  round  from  the  gun  to  target,  and  the 
number  of  rounds  fired  in  each  engagement. 
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The  projectile's  TOF  is  a  function  of  the  initial 
and  final  velocity  of  the  projectile.  The  velocity  of 
the  round  at  impact,  V^,  can  be  determined  as  follows; 

Vf  *  Vi  exp[-pCdARi/2m]  (50) 

where  V.  *  the  initial  muzzle  velocity  of  the  round, 

1  m/sec  (930)  ; 

p  »  density  of  the  air,  kg/m3  (1.2247); 

*  dimensionless  coefficient  of  drag  (.38); 

2 

A  *  cross-sectional  area  of  the  round,  m 
(4.16  x  10-4) . 

R  =  intercept  range,  km;  and 
m  =  mass  of  the  round,  kg  (.195)  (Ref  3:48). 

Substituting  the  average  values  used  in  the  calculations 
(in  parenthesis,  above)  the  equation  reduces  to  the  follow 
ing  in  terms  of  R  kilometers  to  intercept: 

Vf  «  930  exp  [-.4965  Rj  (51) 

From  the  velocity  calculation,  the  TOF,  in  seconds 
can  be  determined: 


TOF  a 


2m 

PcdA 


(52) 


where  the  terms  are  as  previously  defined  (Ref  3:48). 

Substituting  the  average  values  from  the  above  equation 
this  becomes: 


2014.46 


-  2.166 


Leek  and  Schmidt  use  an  average  vulnerable  area 

(Av)  of  55.65  ft2  (5.17  m2)  (Ref  9:39).  This  was  based 

on  an  average  viewing  aspect  of  a  total  projected  area 
2 

of  265  ft  and  a  21  percent  vulnerable  area. 

The  dispersion  of  the  rounds  about  the  aim  point 
in  a  combat  situation  was  assumed  to  be  20  mils  (Ref  9:39). 
This  angular  dispersion  represented  a  one  sigma  standard 
deviation.  Figure  14  depicts  the  angular  deviation  a  in 
terms  of  R.  Converted  to  range  at  intercept,  the  a  becomes 


<r  «  crR 

m 


where  R  *  intercept  range,  km; 

a  *  angular  deviation,  in  mils; 


a  *  20R 
in 


The  single  shot  probability  of  kill  (Pk  )  becomes: 


^ss  2irc2 


«p{-  $r»- «2«  tof2i21?  <«> 

+A„  1  *  L  2 ire  +  A„  -* 


where  g  *  the  number  of  "g's"  the  aircraft  is  undergoing 
for  its  current  position  in  the  battle  area 
(Ref  13:98) . 

Sustituting  the  constant  terms  the  equation  becomes: 


5*17  aim  - 

L  i 

'[9.8 

g  top2]2 

2  axP 

2iro  +5.17 

\  I 

-  2ira2 

+  5.17  . 
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Fig.  14.  AAA  Dispersion  Pattern 

The  overall  Pfc  is  dependent  on  the  and  the 

ss 

number  of  rounds,  n,  fired  by  the  AAA.  Fifty  rounds  would 
be  typical  for  the  combat  arena:  adequate  to  cover  the 
target,  yet  not  enough  to  overheat  and  damage  the  gun 
barrels.  For  the  AAA  weapons  of  the  model,  this  represents 
a  2/3  second  burst.  Finally,  the  Pk  calculation  becomes: 

Pk  -  1.0  -  (1-Pk  )50  (58) 

SS 

or  the  probability  of  AAA  kill,  P.  ,  is  1  minus  the  proba- 
bility  of  the  aircraft  surviving  50  single  round  shots 


Command  and  Control 


The  Soviets  employ  numerous  EW  radars  to  augment 

the  command  and  control  of  the  air  defense  behind  the 

FEBA.  This  thesis  portrays  command  and  control  of  threat 

radars  by  associating  EW  radars  with  randomly  selected 

weapon  system  radars.  The  EW  sectors  are  defined  in  terms 

of  threat  belts  and  the  LOC .  The  EW  radars  and  their 

assigned  coverage  sectors  are  shown  in  Figure  15.  EW 

Nn 

radars  depicted  in  the  figure  as  an  [x]  are  used  for  area 

control  of  radars  located  in  sector  "n"  and  control  all 

systems  types  within  that  sector.  An  EW  radar  depicted 
w 

as  (xj  is  associated  with  a  specific  weapon  system  battery. 

Area  1  extends  from  the  FEBA  to  belt  5  (five 
kilometers  behind  the  FEBA)  and  north  of  the  LOC  (see 
Figure  1) .  Area  2  covers  the  same  east-west  distance 
south  of  the  LOC.  Area  3  extends  from  belt  5  to  belt  9 
(35  kilometers  behind  the  FEBA)  and  north  of  the  LOC. 

Area  4  covers  south  of  the  LOC  and  below  area  3 .  Area  5 
extends  east  from  belt  9  to  the  target  area,  covering  both 
north  and  south  of  the  LOC. 

Each  of  the  area  specific  EWs  provides  early  warn¬ 
ing  information  for  threat  batteries  located  within  its 
area.  By  attacking  the  EW  radars  the  WW  can  disrupt  the 
normal  handoff  communication  between  the  EW  and  its  associ¬ 
ated  threat  radars,  thus  increasing  the  time  required  by 
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Fig.  15.  Early  Warning  Network 


individual  sites  to  acquire,  track  and  fire  on  an  air¬ 


craft. 


The  threat  environment  of  this  model  associates 


50  percent  of  each  threat  type  radars  with  a  specific  area 
EW  radar.  (Note:  The  Soviets  have  many  types  of  EW  sys¬ 
tems.  This  thesis  models  a  general  system  associated  with 

only  a  percent  of  the  total  threat  radars.)  For  example, 

N1 

half  of  the  threat  radars  in  area  1  are  assigned  to  {&]  . 

If  the  WW  kills  this  EW  radar  or  forces  it  to  stop  radiating 
those  radars  associated  with  it  must  operate  autonomously 
and  the  corresponding  target  acquisition  and  tracking  time 
will  be  at  a  maximum  value.  The  remaining  threat  radars 


that  do  not  have  an  associated  EW  in  the  model  will  be 


assumed  to  have  acquisition  and  tracking  times  uniformly 
distributed  between  a  minimum  and  maximum  value. 

Summary  of  the  Defensive 
System  Structure 

The  WW  and  strike  aircraft  attempt  to  penetrate 
an  enemy  defensive  structure  comprised  of  four  types  of 
SAMs,  one  AAA  weapon  system,  and  an  EW  radar  network.  The 
enemy  weapon  systems  attempt  to  engage  the  aircraft  as  soon 
as  possible  and  evaluate  the  probability  of  destroying  the 
target.  System  characteristics  such  as  maximum  detection 
range  and  minimum  detection  altitude  limit  the  distance  at 
which  an  aircraft  can  be  detected.  After  detection  the 
radars  track  the  targets  for  a  specified  time.  The  dura¬ 
tion  of  this  acquisition  and  tracking  time  depends  on 
whether  or  not  the  site  received  track  information  on  the 
aircraft  from  an  EW  site.  During  this  time  period  the 
site  calculates  the  aircraft's  heading,  velocity,  and  both 
the  slant  range  and  bearing  from  the  site. 

The  sites  attempt  to  launch  a  missile  or  fire  on 
the  aircraft  so  as  to  maximize  the  P^.  P^  calculations 
for  SAMs  depend  on  the  range  to  intercept,  R^,  the  air¬ 
craft's  RCS  at  intercept,  and  the  weapon's  lethal  radius. 
AAA  P^  depend  on  the  TOE  of  each  round,  the  dispersion 
pattern  of  rounds  about  the  target,  the  target's  vulner¬ 
able  surface  area,  and  the  number  of  rounds  fired  at  the 
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aircraft  per  engagement.  If  the  weapon's  estimated  is 
greater  than  a  5  percent  threshold  value  the  site  will 
launch  at  or  fire  on  the  target. 

The  model  also  provides  a  method  for  evaluating 
launch  and  impact/intercept  times  for  maneuvering  targets. 

At  the  end  of  acquisition  and  tracking  the  model  calculates 
a  range  to  intercept,  ,  at  the  launch  time.  When  the 
simulation  time  reaches  this  scheduled  launch  time  the  model 
recalculates  and  compares  this  new  R^  to  the  one  calculated 
at  the  end  of  acquisition  and  tracking.  If  the  two  values 
compare  within  10  percent,  the  site  launches  its  missile. 

If  the  two  differ  a  new  launch  time  is  scheduled.  In  a 
similar  manner,  the  Ri  at  impact /intercept  is  calculated 
and  compared  to  the  one  computed  at  launch.  If  the  two 
differ  by  more  than  10  percent,  the  impact/ intercept  is 
rescheduled . 

Summary 

The  chapter  described  each  of  the  three  major  ele¬ 
ments  comprising  the  7ZBA  penetration  air  battle  and 
developed  the  analytical  methodology  of  the  interaction 
between  the  three  elements.  In  the  next  chapter  this 
methodology  will  be  translated  into  a  simulation  model. 
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III.  Simulation  Model 


Overview 

In  the  previous  chapter  the  components  and  charac¬ 
teristics  of  the  WW  defense  suppression  system  were  described 
and  the  system  structure  defined.  Once  this  has  been  done 
the  model  can  be  computerized.  In  this  chapter  the  struc¬ 
tural  model  is  translated  into  the  SLAM  simulation  language. 
First,  simulation  models  are  reviewed  with  emphasis  on  the 
two  basic  timekeeping  orientations.  Next,  the  SLAM 
language  is  introduced  and  the  interfaces  with  the  struc¬ 
tural  model  are  covered.  Finally,  the  simulation  model  is 
presented  in  a  logical,  sequential  manner. 

Simulation  and  Combined 
Simulation  Models 

Simulation  modeling  was  chosen  as  the  methodology 
to  analyze  the  MW  defense  suppression  problem  because  there 
were  no  analytical  methods  available  that  could  represent 
the  dynamic  interactions  and  extreme  complexities  of  WW 
operations.  In  addition,  it  would  be  very  difficult  if 
not  impossible  to  conduct  an  experiment  with  the  necessary 
system  components  of  the  WW  defense  suppression  mission. 
Simulation  modeling  offers  a  viable  alternative  for  ana¬ 
lyzing  the  WW  system. 
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Simulation  modeling  is  a  technique  for  studying 
problems  in  which  a  model  of  the  problem  is  constructed 
resembling  the  system  under  investigation.  After  the 
model  is  developed  experiments  are  conducted  over  the  time 
period  of  interest  simulating  the  operation  of  the  system. 
Data  are  then  gathered  to  estimate  the  characteristics  of 
the  problem. 

Models  can  be  classified  as  either  discrete  change 
or  continuous  change  systems.  The  basic  difference  between 
the  two  systems  is  the  manner  in  which  system  time  is 
modeled.  In  discrete  change  simulations  system  variables 
change  only  at  specified  points  in  simulation  time.  These 
points  are  commonly  called  event  times.  An  example  of  a 
discrete  change  system  is  a  bank  where  the  number  of  cus¬ 
tomers  in  the  bank  changes  only  when  a  customer  arrives  or 
departs .  In  continuous  change  simulations  the  system  vari¬ 
ables  continuously  with  time.  An  example  is  an  airplane 
in  flight.  Because  the  variables  change  continuously  over 
time  differential  equations  are  required  to  define  the  rela¬ 
tionships  between  the  variables. 

Real  systems  are  neither  discrete  nor  continuous 
but  a  combination  of  both.  In  a  combined  discrete- 
continuous  simulation  there  are  three  types  of  interactions 
that  can  occur: 

1.  A  discrete  event  may  cause  a  discrete  change 
in  the  value  of  a  continuous  system  variable. 


2.  A  discrete  event  may  cause  the  relationship 


governing  a  continuous  system  variable  to  change  at  a  par¬ 
ticular  time. 

3.  A  continuous  state  variable  achieves  some  pre¬ 
determined  value  (threshold  value)  which  may  cause  a  dis¬ 
crete  even  to  occur  or  be  scheduled  (Ref  8:47) . 

The  conceptual  framework  for  the  combined  model  is 
that  the  system  cam  be  described  in  terms  of  entities  (such 
as  airplanes  or  radar  systems) ,  their  attributes  (which  are 
characteristics  of  entities  such  as  velocity  and  altitude 
for  an  airplane  or  minimum  and  maximum  effective  range  for 
a  radar  system) ,  and  state  variables  (which  are  the  con¬ 
tinuous  system  variables  that  change  with  time  such  as  an 
airplane's  position  in  flight).  The  behavior  of  the  system 
is  simulated  by  computing  the  values  of  the  state  variables 
at  small  time  steps  and  by  computing  the  values  of  attri¬ 
butes  of  entities  at  event  times  (Ref  12:72).  In  combined 
simulations,  events  can  occur  at  a  scheduled  point  in  time 
or  when  the  system  reaches  a  certain  state.  An  example  of 
the  former  is  when  an  aircraft  is  detected  by  a  radar  sys¬ 
tem.  Its  exact  position  can  only  be  determined  after  a 
given  time  interval,  representing  the  radar's  tracking 
and  acquisition  time.  In  a  simulation,  this  event  would 
be  scheduled  to  occur.  An  example  of  an  event  occurring 
when  the  system  reaches  a  certain  state  is  when  a  WW 
searches  for  a  threat  radar  to  attack.  As  the  WW  proceeds 
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across  the  defense  area,  its  position  continuously  changes. 
As  its  position  changes,  its  line-of-sight  distance  and  the 
distance  to  the  nearest  available  radar  that  the  WW  can 
attack  are  compared.  As  the  line-of-sight  distance  crosses 
the  distance  to  the  nearest  available  radar,  a  state  event 
is  defined  to  have  occurred:  the  radar  has  been  detected 
by  the  WW.  The  possible  occurrence  of  a  state-event  must 
be  tested  at  every  simulation  time  advance. 

SLAM  and  Structural 
Model  Interfaces 

The  WW  defense  suppression  system,  as  described  in 
the  structural  model  presentation  in  Chapter  II,  consists 
of  three  entities — WW  aircraft,  attack  aircraft,  and  threat 
radars.  The  attributes  of  the  aircraft  can  be  considered 
the  following:  call  sign,  velocity,  altitude,  ARM  configura 
tion  (WW  only) ,  turn  rate  (WW  only) ,  radar  cross-section, 
and  heading.  The  attributes  of  the  radar  are  its  sequen¬ 
tial  number,  minimum  effective  altitude,  maximum  effective 
range,  associated  EW  radar  (if  any),  power  configuration, 
and  position.  The  state  variables  for  the  system  are  the 
attack  aircraft's  position  and  the  WW's  position,  heading, 
velocity  components,  distance  to  the  attacked  radar,  line- 
of-sight  distance,  relative  and  absolute  relative  bearing 
to  the  attacked  radar.  Note  that  the  concept  of  a  state 
variable  is  dependent  on  the  viewpoint  of  the  modeler. 


As  previously  noted,  continuous  modeling  involves 
characterizing  a  system's  behavior  through  a  set  of  time- 


dependent  equations  (Ref  12:62).  The  WW's  position  and 
heading  are  described  by  a  set  of  such  equations.  Figures 
16  and  17  depict  how  the  model  solves  for  the  aircraft's 
position  over  each  small  time  increment.  In  Figure  16, 
an  aircraft  is  located  at  point  A  with  an  x  and  y  coordin¬ 
ate  (x^,  y^)  and  a  heading  H.  The  aircraft's  velocity 
vector,  v,  at  point  A  can  be  resolved  into  two  components, 
V  and  V  .  An  increment  of  time  At  later,  the  aircraft 

^  y 

is  located  at  point  B  (x2,  y2) .  The  model  solves  for  the 
value  of  the  new  coordinates  as  follows: 


x2  =  Xj  +  v  At 


X2  *  X1  +  (V  cos  0)  At 


Similar  evaluations  of  the  y  components  of  heading  allows 
the  aircraft  to  "fly"  across  the  battle  area. 

The  SLAM  simulation  language  was  chosen  for  this 
problem  because  of  its  power  and  flexibility  to  incorporate 
a  combined  simulation  model  necessary  for  the  MW  defense 
suppression  mission.  The  many  dynamic  interactions  and 
inherent  complexities  of  the  WW  system  required  a  language 
that  was  capable  of  representing  model  processes  in  a 
simple  and  direct  manner. 


B(x2'y2) 


Fig.  16.  Continuous  Model  Diagram 


Vy-V  sin  0 


Fig.  17.  Continuous  Tins  Velocity  Components 
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In  the  WW  structural  model,  aircraft  fly  through 
the  defense  area  enroute  to  the  target  area  (attack  air¬ 
craft)  or  fly  through  the  area  in  pursuit  of  threat  radars 
to  attack  (WW) .  SLAM  portrays  this  process  of  entities 
flowing  in  a  system  by  employing  the  concept  of  a  network 
structure.  The  network  consists  of  special  symbols  called 
nodes  and  branches.  These  symbols  model  elements  in  the 
WW  system  where  an  activity  (branch)  or  an  event  (node) 
could  be  realized.  The  network  is  essentially  a  pictorial 
representation  of  the  system.  The  use  of  a  network  struc¬ 
ture  to  describe  the  system  structure  is  especially  helpful 
to  the  modeler  in  getting  the  "whole  idea"  concept  of  the 
model. 

Aircraft  (entities),  therefore,  are  modeled  as 
flying  through  the  defense  area  where  processes  occur. 

The  processes  that  aircraft  undergo  are  attack  by  threat 
radars  or,  in  the  case  of  the  WW,  attack  on  the  WW  by  threat 
radars  and  the  WW  attack  on  the  threat  radars. 

SLAM  incorporates  discrete  change  simulation  by 
allowing  the  modeler  to  define  an  event  and  the  changes 
that  occur  in  the  system  at  discrete  points  in  time.  The 
time  where  the  state  of  the  system  changes  is  called  an 
event  time  and  the  associated  logic  for  processing  the 
change  is  called  an  EVENT  in  the  SLAM  notation.  The  WW 
attack  profile  logic  is  modeled  using  events,  specifically 
events  1  through  13.  See  Appendix  A. 


The  interface  between  the  network,  discrete  event, 
and  continuous  models  is  crucial  in  modeling  discrete- 
continuous  change  systems.  Continuous  change  variables 
interact  with  discrete  change  variables  through  a  state- 
event.  This  mechanism  is  nothing  more  than  an  event  which 
is  triggered  by  a  continuous  state  variable  reaching  a  pre¬ 
determined  value.  The  predetermined  value  is  commonly 
called  a  threshold  value. 

The  SLAM  DETECT  node  is  the  primary  interface 
between  the  network  and  continuous  portions  of  the  combined 
model  and  is  used  extensively  in  the  WW  attack  profile. 

The  duration  that  the  WW  is  engaged  in  certain  activities 
is  keyed  to  the  release  of  the  appropriate  DETECT  node. 

For  example,  when  the  WW  initially  detects  a  threat  radar 
a  turn  direction  is  computed  for  the  WW  as  it  begins  its 
ranging  routine  on  the  threat  (EVENT  1  in  the  simulation 
model) .  The  duration  or  time  the  WW  is  in  that  turn  is 
keyed  to  the  release  of  a  detect  node  (DETECT  node  W17 
for  WW1) .  When  the  continuous  state  variable,  in  this 
example  the  absolute  relative  bearing  to  the  threat  radar 
(state  variable  SS(49)  for  WW1) ,  reaches  the  predetermined 
value  of  75  degrees  then  the  duration  of  the  first  ranging 
routine  is  completed  and  the  WW  reaches  the  next  network 
node  (EVENT  2  for  the  example) .  See  Appendix  A  for  the 
SLAM  network. 
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SLAM  incorporates  numerous  modeling  mechanisms  to 
assist  the  modeler  in  translating  the  system  into  the 
simulation  model.  These  mechanisms  are  invaluable  tools 
for  not  only  making  the  computerization  of  the  model 
straightforward  but  also  for  helping  crystalize  a  system's 
complex  interactions  into  clear  and  simple  programming. 

The  more  important  mechanisms  are  presented  in  the  follow¬ 
ing  paragraphs  while  the  remainder  are  discussed  in  the 
simulation  model  itself. 

A  simulation  model  is  built  for  a  specific  purpose. 
Once  the  model  is  built  system  variables  are  changed  to 
test  various  conditions  of  the  model.  The  ease  in  which 
these  variables  can  be  changed  is  a  test  of  the  language's 
flexibility.  SLAM  incorporates  the  concept  of  a  global 
variable  (XX (i)  in  SLAM  notation)  to  provide  this  flexi¬ 
bility.  System  variables  such  as  aircraft  altitude  or 
velocity  can  be  defined  in  the  simulation  model  as  global 
variables.  The  global  variables  can  be  changed  to  another 
value  by  the  modeler  in  subsequent  experimental  runs  of 
the  model.  For  example,  WW  altitude  is  defined  in  the 
simulation  model  equal  to  global  variable  XX  (3)  (see 
Appendix  A) .  For  the  first  5  experimental  runs  WW  altitude 
is  to  60  meters  (XX (3) *60,  see  line  172  of  the  SLAM  program 
in  Appendix  A).  For  the  following  5  experimental  runs, 
altitude  is  changed  to  200  meters  (XX (3) *200,  see  line  178 
of  the  SLAM  program) .  The  ease  and  efficiency  of  changing 
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system  variables  in  this  manner  demonstrates  the  power  and 
flexibility  of  SLAM. 

Entities  are  generated  in  the  model  and  routed 
through  the  network  by  CREATE  nodes.  At  these  nodes,  the 
time  of  an  entity's  first  creation,  the  time  between  entity 
creations,  and  the  number  of  entity  creations  are  specified. 
For  the  simulation  model  two  CREATE  nodes  are  used,  one 
for  generating  WWs  and  the  other  for  attack  aircraft. 

Entities  have  unique  characteristics  associated 
with  them  which  SLAM  denotes  as  attributes.  In  the  simula¬ 
tion  model  WW  aircraft  have  two  explicit  attributes: 
call  sign  (211  or  212)  and  ARM  configuration.  These 
attributes  flow  with  their  respective  entities  in  the  net¬ 
work.  SLAM  designates  attributes  at  ASSIGN  nodes.  In  the 
model  ASSIGN  node  WAS  prescribes  attributes  for  the  WW 
and  ASSIGN  node  AAS  for  the  attack  aircraft. 

Simulation  Model 

The  WW-threat  environment  system  structure  was 
modeled  into  the  SLAM  network  and  is  shown  in  Appendix  A. 

The  network  can  be  considered  to  be  in  two  parts:  a  WW 
attack  network,  where  the  WW  hunts  for  threats  to  attack; 
and  a  radar  attack  network,  where  threat  radars  attack 
aircraft  flying  through  the  defense  area. 

Two  WWs  enter  the  engagement  through  the  CREATE 
node  with  the  time  between  arrivals  specified  as  global 


76 


variable  XX  (8).  The  WWs  enter  the  conflict  area  50  km 
prior  to  the  FEBA  to  allow  sufficient  time  to  search  for 
threats.  The  HW's  attributes  are  specified  in  the 
ASSIGN  node  WAS.  Attribute  one  is  the  WW's  tail  number 
(211  or  212)  and  attribute  two  is  its  ARM  allocation. 

The  first  digit  of  attribute  two  indicates  the  number  of 
AGM-  Bs  and  the  second  digit  the  number  of  AGM-As .  Each 
WW  enters  the  area  ±500  meters  of  the  centerline  that 
divides  the  area  into  equal  north  and  south  components. 

The  attack  force  aircraft  are  created  at  the  second 
CREATE  node  with  time  between  creations  equal  to  XX (7). 

A  total  of  ten  attack  aircraft  are  created.  As  with  the 
WWs,  the  attack  force  enters  the  area  50  km  prior  to  the 
FEBA.  Attribute  one  for  the  attack  force  is  assigned  at 
node  AAS  and  represents  the  call  sign  of  the  aircraft 
(1  to  10).  The  attack  force  also  enters  ±500  meters  of 
the  FEBA  centerline. 

WW-Attack  Profile 

Continuous  modeling  concepts  are  used  to  simulate 
the  WW-attack  profile  of  the  model.  This  essentially 
involves  characterizing  the  behavior  of  the  WW  to  a  system 
of  equations.  As  the  status  of  the  WWs  change  in  the 
system,  the  equations  that  describe  the  WW  also  change. 

After  the  WWs  are  created  they  are  separated  at 
node  G1  based  on  their  tail  number  and  go  to  either 
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GWl  or  GW2.  (All  aircraft  are  routed  to  node  RATK  for  the 
radar-attack  phase  of  the  model.  This  network  will  be 
described  later  in  the  chapter.)  Each  WW  is  then  routed 
through  a  sequence  of  events  determined  by  whether  certain 
conditions  are  detected  by  key  model  parameters.  These 
conditions  roughly  correspond  to  the  phares  in  the  struc¬ 
tural  model.  The  network  phases  are  identical  for  each  WW. 

At  node  Gl,  WWl  begins  an  activity.  The  duration 
of  the  activity  from  node  Gl  to  node  RNl  is  specified  as 
REL(WIS).  REL  is  a  release  specification  and  it  is  used 
throughout  the  WW-attack  portion  of  the  network.  When  the 
node  upon  which  the  release  specification  is  realized,  in 
this  case  node  W1S,  then  the  activity  is  completed.  Thus, 
the  time  it  takes  for  WWl  to  get  from  Gl  to  RNl  is  keyed 
to  the  release  of  node  W1S.  WlS  is  a  DETECT  node  whose 
condition  for  release  is  that  the  range  of  WW1(SS(37)) 
decreases  to  less  than  WWl's  radar  line-of-sight  (SS  ( 41 ) ) . 
WWl  is  held  in  the  activity  between  Gl  and  RNl  until  WlS 
is  detected.  Continuous  state  equations  describe  the  WW's 
position  as  time  advances  and  are  explained  in  the  para¬ 
graph  on  subroutine  STATE.  At  each  new  iteration  of  time, 
the  WW's  position  is  changed  and  ranges  to  all  candidate 
threat  radars  are  evaluated.  Finally,  when  the  WW's  range 
is  less  them  the  line-of-sight  range,  WlS  is  detected  which 
in  turn  releases  REL (WlS) .  WWl,  along  with  its  attributes 
proceeds  to  node  RNl. 


At  RN1  event  1  is  called  (line  31) .  Event  1  turns 
the  WW  in  the  shortest  direction  to  begin  its  triangula¬ 
tion  routine  on  the  radar.  Turn  rate  for  the  WW  is  spe¬ 
cified  as  RATE  and  is  set  to  four  degrees  per  second 
(simulating  a  60°  bank,  level  turn).  The  WW's  working 
designator,  SS(I+12),  is  set  to  1  to  key  subroutine  STATE 
from  computing  new  threat  parameters  when  the  WW  has  just 
started  to  engage  a  radar.  The  radar  the  WW  is  attacking 
has  its  14th  attribute  reset  to  the  WW's  call  sign,  replaced 
in  Pile  2,  and  then  copied  into  the  WW's  attack  file. 

WW1  now  starts  the  activity  from  RN1  to  RE1.  The 
release  specification  for  this  activity  is  W17.  W17  detects 

an  absolute  relative  bearing  (SS(49))  of  75  (degrees)  and 
sends  WWl  to  event  2.  Event  2  rolls  the  WW  out  of  its  turn 
(RATE*0 ) .  The  next  activity  is  begun  and  is  keyed  to 
detect  node  W10  which  is  realized  when  the  WW's  absolute 
relative  bearing  increases  to  105.  Upon  detection,  WWl  is 
sent  to  event  3  (lines  447-454)  which  now  turns  the  WW 
back  into  the  threat  (RATE«±4  depending  on  turn  direction) . 
WWl  now  waits  until  detect  node  WIT  is  realized  which  occurs 
when  the  aircraft  is  within  10  degrees  relative  bearing  to 
the  threat.  When  WIT  is  released  event  4  is  called  (lines 
455-464)  which  reduces  the  WW's  turn  rate  to  2  degrees  per 
second.  Following  event  4,  WWl  begins  an  activity  until 
detect  node  W1B  is  realized.  This  occurs  when  WWl's 


heading  (SS(29))  and  relative  bearing  to  the  threat  (SS(45)) 
are  identical  (aircraft  boresighted  on  the  radar) . 


Event  5  (lines  466-529)  begins  by  selecting  which 
ARM  to  fire.  If  the  WW  is  too  close  to  fire  an  ARM  (range 
to  threat  is  less  than  minimum  range  of  the  ARM)  then  the 
WW  is  sent  back  into  the  network  for  a  "repositioning 
routine"  (events  10  and  11).  Otherwise,  ARM  firing  range 
is  determined  based  on  minimum  and  maximum  ranges  of  the 
ARM  and  the  distance  from  the  WW  to  the  threat.  Flight 
time  of  the  ARM  is  computed  based  on  ARM  firing  range  and 
velocity.  P^  is  determined  from  a  random  sample  draw. 

The  time  for  the  WW  to  get  to  the  ARM  release  point  is  cal¬ 
culated  (RLWW)  and  is  loaded  into  XX (15)  (XX (IR) ) .  XX (16) 
is  set  equal  to  the  ARM's  flight  time  (TOF) .  The  WW's 
turn  rate  is  set  to  zero.  The  WW  now  re-enters  the  network. 

Back  in  the  network,  WWl  can  take  one  of  two  paths. 
If  the  WW  is  too  close  t6  fire  an  ARM  (range  to  threat  less 
than  the  minimum  ARM  range,  SS(38)  less  than  800)  then  it 
begins  a  repositioning  routine  to  place  it  in  firing  range 
and  takes  the  branch  to  node  WR2  with  a  release  specifica¬ 
tion  of  WR2.  (Repositioning  routine  network  will  be  dis¬ 
cussed  later  in  the  chapter.)  Otherwise  the  WW  takes 
the  branch  to  event  6  with  a  duration  specified  as  XX  (15) , 
the  time  to  get  WWl  to  ARM  release  point. 

At  event  6  the  WW  launches  an  ARM.  Each  of  the 
WW's  four  ARMs  has  its  own  TOF  designator.  For  WWl,  ARM 
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number  one's  TOF  is  set  to  XX(85),  number  two's  to  XX(86), 
number  three's  to  XX(87),  and  number  four's  to  XX(88). 

Thus  the  capability  exists  in  the  model  to  have  a  WW 
launch  all  four  ARMS  simultaneously.  The  TOF  for  each  ARM 
is  used  to  evaluate  its  respective  at  impact  time.  This 
is  accomplished  by  calling  the  appropriate  ENTER  node.  For 
WWl  ENTER  node  1  is  called  with  the  attributes  of  the  radar 
being  passed  along  to  the  ENTER  node.  After  the  WW  is  sent 
back  to  the  network ,  a  check  is  made  to  see  if  there  are 
any  ARMs  remaining.  If  the  WW  is  out  of  ARMs,  then  its 
working  designator  is  set  to  nine  to  key  subroutine  STATE. 

WWl  goes  back  to  the  network  and  can  branch  in  one 
of  two  directions.  If  it  no  longer  has  any  ARMs  then  the 
WW  is  sent  to  node  WGH  to  begin  a  "go  home  routine." 
Otherwise,  after  a  five-second  delay  to  simulate  release 
of  the  ARM,  the  WW  goes  to  event  7.  Here  its  working 
designator  is  set  to  zero  to  key  subroutine  STATE  to  select 
a  new  radar  to  attack.  After  a  0.1  second  delay  in  the  net¬ 
work,  which  allows  STATE  to  select  a  new  threat,  event  8 
is  reached. 

Event  8  enters  the  WW  back  into  the  network  depend¬ 
ing  on  its  distance  to  the  new  radar.  If  the  distance  to 
the  threat  is  less  than  line-of-sight  range  (for  WWl, 

SS  (37)  less  them  SS(41)),  then  the  triangulation  routine 
is  started,  the  WW's  working  designator  set  to  one,  a  turn 
direction  is  computed,  and  the  network  is  re-entered.  If 
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the  distance  to  the  threat  is  beyond  the  WW's  line-of -sight 
range  then  no  activities  are  performed  by  the  WW  in  event  8. 
When  the  WW  returns  to  the  network  from  event  8  it  can 
take  one  of  two  branches.  If  the  distance  to  the  threat 
was  less  than  its  line-of-sight  distance,  then  the  WW 
branches  to  node  RE1,  with  a  release  specification  of  W1C, 
absolute  relative  bearing  passing  75.  (Node  RE1  is  event  2 
where  the  remainder  of  the  triangulation  routine  is  per¬ 
formed.)  If  the  WW  does  not  detect  the  new  threat,  which 
was  selected  by  subroutine  STATE  as  the  closest  threat  to 
the  WW,  then  the  WW  branches  to  node  GWl,  where  it  essen¬ 
tially  begins  the  process  of  searching  for  threats  to 
attack  once  more. 

As  described  previously  in  the  paragraph  on 
event  6,  when  the  WW  launches  an  ASM,  the  appropriate 
ENTER  node  is  called.  The  network  for  ENTER  nodes  one  and 
two  simulate  the  flight  path  of  the  ARM.  At  ENTER  node 
one,  the  ARM  can  take  one  of  four  paths  to  node  WRK  depend¬ 
ing  on  whether  it  is  the  first,  second,  third,  or  fourth 
ARM  launched  by  WW1  (XX (17) -1,2, 3,  or  4).  The  duration 
of  the  activity  for  the  ARM  to  get  from  the  ENTER  node  to 
WRK  is  equal  to  the  TOF  of  the  ARM  which  was  calculated 
in  event  5  and  set  equal  to  the  appropriate  global  vari¬ 
able  (XX(85,86,87,  or  88)). 

At  node  WRX  event  9  occurs  which  evaluates  the 
of  the  ARM.  If  the  WW  kills  the  radar,  then  the  following 
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sequence  of  events  takes  place.  First,  if  the  radar  was 
an  EW  radar  (ATRIB (1) *1)  then  all  threat  radars  associated 
with  it  have  their  5th  attribute  set  to  10  to  key  tracking 
and  acquisition  time  computations,  which  are  described  in 
the  radar-attack  portion  of  the  model.  Next,  the  destroyed 
radar  is  removed  from  both  the  current  radar  file  (File  1) 
and  the  WW  available  file  (File  2) .  All  events  associated 
with  this  radar  are  removed  from  the  event  calendar. 
Finally,  XX (49)  is  set  to  attribute  one,  the  radar's 
sequential  number  in  the  model. 

If  the  radar  is  not  killed,  then  its  14th  attribute 
is  reset  to  zero  to  indicate  that  the  radar  is  once  again 
available  for  engagement  by  a  WW. 

From  event  9  the  ARM/WW  will  go  to  node  WKL  if  the 
radar  was  killed.  It  then  enters  an  assign  node  where 
XX (50)  counts  the  number  of  radars  killed  by  the  WW. 

If  after  the  WW  completes  the  triangulation  routine 
and  finds  itself  too  close  to  fire  an  ARM  due  to  the 
minimum  range  of  the  ARM  (event  5)  then  a  repositioning 
routine  is  started.  This  routine  begins  in  event  5  which 
initially  turns  the  WW  away  from  the  site  after  the  range 
to  the  radar  was  determined  to  be  too  small  to  fire  the 
ARM.  From  event  5  the  network  is  entered  and  WW1  is  sent 
to  node  WWR  after  waiting  until  node  W1R  is  released,  which 
occurs  when  the  absolute  relative  bearing  to  the  site  is 
179.  This  position  places  the  radar  at  the  WW's  6  O'clock 
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position  as  the  WW  attempts  to  get  sufficient  distance  to 
fire  its  ARM.  Node  WWR  is  event  10  which  rolls  the  WW 
out  of  its  turn.  From  event  10,  WW1  enters  the  network  with 
an  activity  duration  whose  release  specification  is  W1D. 
Node  W1D  is  realized  when  WWl's  range  to  the  radar  passes 
18,000  meters.  WW1  proceeds  to  event  11  where  it  is  turned 
back  into  the  site  to  re-attack.  When  WWl  is  bore sighted 
on  the  site  (W1E) ,  event  5  is  reached  once  again  where 
an  evaluation  for  an  ARM  launch  is  made. 

When  the  WW  fires  its  last  ARM  it  is  sent  back 
across  the  FEBA  via  node  WGR,  event  13,  and  node  WHM. 

Global  variable  XX (48)  counts  the  WWs  as  they  reach  the 
home  point. 

The  network  logic  for  WWl  and  WW2  is  identical. 

Each  WW,  however,  has  its  own  separate  network  where  only 
activities  that  pertain  to  it  can  occur.  The  logic  for 
the  events  that  describe  the  WW-attack  scenario  are  the 
same  no  matter  if  the  aircraft  is  WWl  or  WW2. 

Radar-Attack  Profile 

In  the  radar-attack  po’ tion  of  the  network  threat 
radars  engage  both  attack  force  and  WW  aircraft.  As  each 
aircraft  enters  the  network  the  range  to  the  closest  radar 
that  isn't  trackingan  aircraft  is  calculated.  If  the 
range  is  less  than  the  maximum  effective  range  of  the 
threat  radar,  then  a  sequence  of  events  is  started  that 
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simulates  the  profile  for  a  threat  system  to  launch  a  SAM 
or  shoot  a  AAA  at  the  target  aircraft. 

All  aircraft,  after  they  are  created,  proceed  to 
node  RATK,  where  event  14  is  called.  This  event  simulates 
the  radar  search  phase  of  the  system.  If  the  aircraft 
has  reached  the  target  area  then  global  variable  XX (55) 
is  set  to  the  call  sign  of  the  aircraft  to  key  the  network 
and  the  event  is  terminated.  Similarly,  if  there  are  five 
threats  already  engaged  with  the  aircraft,  then  radars  are 
not  allowed  to  search  for  this  particular  aircraft  and  the 
event  is  terminated.  If  these  two  conditions  are  not 
encountered  then  subroutine  SEARCH  is  called. 

Subroutine  SEARCH 

Subroutine  SEARCH  simulates  the  actual  threat 
radars  searching  for  aircraft  as  they  fly  through  the 
threat  environment.  When  an  aircraft  enters  SEARCH  every 
radar  is  evaluated  to  determine  if  it  can  engage  the  air¬ 
craft.  Four  checks  are  first  made  to  determine  if  the 
radar  can  engage  the  aircraft.  If  the  radar  is  an  EW 
radar  or  if  the  aircraft's  altitude  is  below  the  radar's 
minimum  altitude  or  if  the  radar  is  already  engaged  or  if 
the  radar  is  not  operating;  if  any  of  these  conditions  are 
met,  then  the  radar  cannot  engage  the  aircraft.  If  the 
radar  passes  these  four  checks,  then  the  distance  is  greater 
than  the  multi-path  range  then  the  radar  will  be  unable 
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to  track  the  aircraft.  Finally,  if  the  distance  from  air¬ 
craft  to  radar  is  greater  than  the  maximum  effective  range 
of  the  radar,  then  this  too  excludes  the  radar  from  track¬ 
ing  the  aircraft.  If  the  radar  passes  all  of  these  checks 
then  the  model  allows  the  radar  to  start  the  engagement 
sequence.  The  radar's  8th  attribute  is  set  to  the  air¬ 
craft's  call  sign,  the  radar,  along  with  its  attributes, 
are  filed  in  the  radar  aircraft  file  (LF) ,  the  radar,  with 
its  new  8th  attribute,  is  copied  back  into  File  1,  and  an 
acquisition  and  tracking  time  is  computed.  This  time,  TRC, 
is  computed  from  a  uniform  distribution  based  on  the  threat 
radar's  minimum  and  maximum  acquisition  time  unless  the 
threat  radar's  associated  EW  radar  was  killed  previously 
by  a  WW.  In  this  case,  TRC  is  set  to  the  radar's  maximum 
acquisition  time.  Discrete  event  15  is  scheduled  at  the 
end  of  the  tracking  and  acquisition  time  which  will  make 
the  initial  calculations  of  the  weapon  P^. 

Radar-Attack  Profile — continued 

Subroutine  SEARCH  returns  to  event  14  where  the 
network  is  entered.  From  event  14  the  aircraft  takes  one 
of  three  paths.  If  the  aircraft  reached  the  target  area 
then  it  will  go  to  node  TGT.  If  the  aircraft  was  killed 
by  a  threat  weapon  system  in  discrete  event  17,  which  will 
be  described  later  in  the  chapter,  then  the  aircraft  will 
take  the  path  to  node  AKL.  If  the  above  two  conditions 
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are  not  met,  then  the  aircraft  will  proceed  back  to  node 
RATX  after  a  one-second  activity  duration.  Thus  radars 
will  discretely  search  for  the  aircraft  as  they  fly 
through  the  threat  system  in  one-second  intervals.  If 
the  aircraft  either  reaches  the  target  area  or  is  killed 
then  the  screening  process  stops  for  that  aircraft. 

Discrete  events  15,  16,  17,  18,  and  19  apply  to 
the  radar-attack  portion  of  the  system.  These  events  are 
FORTRAN  code  using  SLAM  subroutines  that  evaluate  the 
radar's  attack  profile  on  the  engaged  aircraft. 

Event  15,  which  is  scheduled  from  subroutine 
SEARCH,  occurs  at  the  end  of  the  tracking  and  acquisition 
time,  TRC.  Here,  the  radar  first  determines  if  it  has  a 
chance  of  getting  a  shot  at  the  aircraft.  Subroutine 
PROB  is  called  which  computes  the  of  the  weapon,  the 
estimated  range  to  intercept  (RI) ,  time  of  weapons  launch 
(TL) ,  and  time  of  impact  (TI).  If  the  P^  is  below  5  per¬ 
cent,  then  the  radar  is  disengaged  from  the  profile  and 
event  19  is  scheduled,  which  releases  the  radar  from  the 
aircraft.  Otherwise,  the  radar-attack  profile  is  continued. 
If  TL  is  equal  to  the  current  simulation  time,  TNOW,  then 
the  threat  radar  "fires"  its  associated  weapon  and  event 
17  is  scheduled  to  occur  at  impact  time.  If  TL  is  greater 
than  TNOW,  then  event  16  is  scheduled  at  estimated  launch 
time. 
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Event  16  occurs  at  launch  time,  TL.  Subroutine 


PROB  passes  new  estimates  of  P^,  RI,  TL,  and  TI  back  to 
the  event.  Once  again,  if  the  P^  is  below  5  percent  then 
event  19  is  scheduled  to  disengage  the  radar  and  aircraft. 

If  the  P^  is  acceptable  event  17  is  scheduled  to  occur  at 
the  new  time  of  impact. 

Event  17  is  realized  at  weapon's  impact  time. 
Subroutine  PROB  is  called  one  last  time  for  the  final  evalua¬ 
tion  of  P^  and  RI .  If  the  aircraft  turned  away  from  the 
site  when  the  weapon  was  airborne  then  RATIO  is  computed 
to  determine  how  much  its  position  has  changed.  If  RATIO 
is  greater  than  1.1  then  event  17  is  rescheduled  with  time 
of  occurrence  set  how  long  it  takes  the  weapon  to  get  to 
the  new  intercept  range.  If  the  aircraft  did  not  turn  or 
did  turn  but  the  turn  was  not  appreciable,  then  a  kill 
determination  is  made  by  using  a  random  sample  draw.  If 
the  aircraft  is  killed,  then  XX (54)  is  set  to  the  call  sign 
of  the  aircraft  to  key  the  network.  Global  variable 
XX(59)  or  XX(57)  is  incremented  as  appropriate  and  all 
events  associated  with  the  destroyed  aircraft  are  removed 
from  the  event  calendar.  In  any  case,  event  18  is  scheduled 
to  occur  in  30  seconds,  representing  the  delay  time  of 
disengaging  the  radar  from  the  aircraft. 

Event  18  frees  the  appropriate  radars  to  re-engage 
new  targets.  If  the  aircraft  was  not  killed,  ATRIB(6) 
not  equal  to  1,  then  only  that  radar  that  fired  at  the 
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aircraft  is  released.  If  the  aricraft  was  killed,  then  all 
radars  associated  with  it  are  disengaged. 

Subroutine  PRQB 

This  subroutine  calculates  the  probability  of  kill 
for  the  threat  weapon  system,  the  range  to  intercept  the 
target  aircraft  (RI) ,  the  time  of  weapons  launch  (TL) , 
and  the  time  of  weapons  impact  (TI).  First  the  attributes 
of  the  aircraft  are  decoded  (lines  979-992) ,  the  multi-path 
range  calculated,  and  the  appropriate  radar-attack  file 
obtained.  The  distance  from  radar  to  aircraft  is  compared 
to  the  multi-path  range  and  the  maximum  effective  weapon 
range.  If  the  radar -aircraft  distance  is  greater  than 
either  of  the  above  then  is  set  to  zero  and  the  sub¬ 
routine  terminated. 

For  the  SAM  threats  the  following  is  accomplished. 
The  relative  bearing  between  the  radar  and  the  aircraft 
is  computed  (ANG,  lines  1014-1025) .  Based  upon  whether 
the  aircraft  is  ahead  or  behind  the  threat  the  appropriate 
values  of  TL,  TI,  and  RI  are  calculated.  Aircraft  radar 
cross-section  is  obtained  from  File  19  based  on  the  aspect 
angle  between  the  threat  and  the  aircraft  (lines  1056-1061) . 
Next,  P^  is  calculated.  First,  if  RI  is  greater  than  the 
maximum  effective  radius  of  the  weapon,  P^  is  set  to  zero. 
Otherwise,  CEP  is  computed  based  on  whether  the  aircraft 


is  ECM  (attack  force)  or  non- jamming  (WW) .  Using  CEP, 

Pk  is  calculated. 

For  the  AAA  threats  the  following  applies.  "G" 
forces  are  determined  based  on  whether  the  aircraft  is  a  WW 
(turning  or  straight  and  level)  or  an  attack  force  (in  or 
outside  the  expected  area  of  AAA  concentration) .  P^ 
calculations  are  given  in  lines  1092-1100. 

Subroutine  INTLC 

Subroutine  INTLC  initializes  variables  to  their 
starting  values.  In  addition  the  radar  file.  File  1,  and 
the  WW-attack  file.  File  2,  are  built.  Radar  cross- 
section  data  is  built  and  stored  in  File  19. 

Subroutine  STATE 

Subroutine  STATE  contains  the  equations  that  define 
the  continuous  set  of  variables  for  the  model.  These  equa¬ 
tions  explicitly  characterize  the  model  as  it  changes  over 
time.  This  time-dependent,  event-dependent  portrayal  is 
fundamental  to  the  concept  of  the  system  as  a  combined 
simulation  model. 

Rate  equations  for  the  attack  force  aircraft  con¬ 
sist  of  just  one  differential  equation.  Because  the 
attack  aircraft  never  change  heading  and  their  heading 
is  090  degrees  in  the  model's  coordinate  system,  the 
velocity  of  the  attack  aircraft,  DD(I) ,  is  defined  to  be 
a  constant,  XX (14). 


Rate  equations  for  the  WW  are  only  slightly  more 
complicated.  The  WW  will  have  velocity  components  in  both 
the  x  and  y  direction.  In  addition,  because  the  WW  changes 
heading  (direction)  as  it  works  its  way  through  the  threat 
environment/  its  heading  rate  of  change  is  expressed  as 
a  differential  equation.  Line  39  gives  the  differential 
equation  for  this  heading  rate  change  as  equal  to  RATE. 

The  next  two  equations  correct  the  WW's  heading  to  the 
limits  of  0°  to  360°.  Lines  45-46  define  the  velocity  of 
the  WW  in  the  x-  and  y-direction.  Next,  the  WW  working 
designator  is  evaluated.  If  the  WW  is  engaged  in  an  attack 
profile,  working  designator  between  one  and  eight,  then  the 
WW  does  not  enter  the  radar  selection  phase  of  STATE 
(lines  62-75)  but  instead  proceeds  to  "101"  to  calculate 
the  remaining  STATE  variables.  If,  on  the  other  hand,  the 
working  designator  is  equal  to  nine,  indicating  that  the 
WW  is  out  of  ARMs,  then  a  designated  go  home  point  is 
chosen  (lines  53-57).  If  the  WW's  working  designator  is 
zero,  signifying  that  the  WW  is  searching  for  threats  to 
attack,  then  the  radar  selection  routine  is  performed  (lines 
62-75) .  In  this  routine  the  closest  radar  that  the  WW  can 
attack  from  File  2,  the  WW-attack  file,  is  selected,  its 
sequential  number  placed  in  "L"  and  its  coordinates  set  to 
XX(I+11)  and  XX (1+12).  Line  76  calculates  the  distance 
from  the  designated  radar  to  the  WW.  Line  78  is  the  line- 
of -sight  range  to  the  WW.  The  relative  bearing  is 
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calculated  in  lines  79-88  based  on  the  position  and  heading 
of  the  WW  with  respect  to  the  radar.  The  absolute  relative 
bearing  is  calculated  in  lines  89-90  and  corrected  for  the 
limits  of  0°  to  180°. 

Appendix  F  lists  a  summary  of  the  events  and 
subroutines . 

In  this  chapter  the  structural  model  was  trans¬ 
lated  into  the  simulation  model  using  the  SLAM  language. 
Basic  concepts  of  simulation  modeling  were  presented  and 
the  interfaces  of  SLAM  and  the  structural  model  discussed. 
Once  the  simulation  model  has  been  developed  tests  are 
ready  to  begin.  The  next  chapter  demonstrates  the  capa¬ 
bility  of  the  model  to  analyze  WW  operations  through  experi¬ 
mentation. 


IV.  Experimentation 


Overview 

The  purpose  of  the  research  effort  was  to  develop 
a  methodology  that  could  be  used  to  analyze  WW  defense  sup¬ 
pression  operations.  Once  the  structural  model  has  been 
translated  into  the  simulation  model  then  the  simulation 
model  can  be  executed  and  an  experiment  conducted  to  demon¬ 
strate  the  capability  of  the  model  to  analyze  WW  operations 
In  this  chapter  experimentation  of  the  simulation  model  is 
performed.  First,  the  data  collection  phase  is  presented: 
WW  variables  that  may  be  analyzed  in  the  model  are  dis¬ 
cussed  and  two  are  selected  for  the  experiment;  experi¬ 
mental  design  is  set  forth;  sample  size  for  the  experiment 
is  determined.  In  the  next  phase  the  experiment  is  con¬ 
ducted  and  the  results  are  discussed.  In  the  final  phase, 
the  model  is  validated  in  a  three-step  process. 

Phase  One:  Data  Collection 

The  simulation  model  was  constructed  with  the  capa¬ 
bility  to  analyze  the  WW  defense  suppression  mission  firmly 
in  mind.  To  that  end,  SLAM  global  variables  were  used  so 
that  various  system  variables  could  be  changed  quickly  and 
effortlessly  to  analyze  different  aspects  of  the  WW  mission. 
Table  VI  lists  these  variables  and  their  associated  global 


TABLE  VI 


VARIABLES  THAT  CAN  BE  CHANGED /EVALUATED 


Variable 

Global  Variable 

ww 

Altitude 

XX  (3) 

Velocity 

XX  (5) 

Interval /Spacing 

XX  (8) 

Number  in  System 

XX  (6) 

ARM  Configuration 

XX (81-82) 

ARM  Launch  Range 

XX  (4) 

Entry  Position 

SS  (25,29) * 

Attack  Aircraft 

Altitude 

XX (13) 

Velocity 

XX (14) 

Interval /Spacing 

XX  (7) 

Number  in  System 

XX  (9) 

Entry  Position 

XX  (61-70) 

Radar  Systems 

Dispersion  on  LOC 

XX  (52) 

•State  Variable — see  text. 
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variable.  (The  entry  positions  of  WW1  and  WW2  are  changed 
with  state  variables  SS(25)  and  SS(29),  respectively.) 
Through  the  use  of  global  variables  any  of  the  system 
variables  listed  in  the  table  can  be  changed  and  experiment 
conducted  to  analyze  their  effect  on  the  output  data. 

To  demonstrate  the  capability  of  the  simulation 
model  two  variables  were  selected  for  the  experiment:  WW 
altitude  and  WW  tactic.  These  two  variables  were  suggested 
by  WW  experts  at  George  AFB  (Refs  10;  11)  as  having  sig¬ 
nificant  impact  on  WW  defense  suppression  operations.  The 
two  altitudes  chosen  for  the  experiment  were  60  meters  and 
200  meters.  At  a  relatively  higher  altitude  the  WW  can 
"see"  farther  with  its  radar  homing  and  warning  equipment 
and  will  be  able  to  begin  an  attack  on  threat  radars  sooner 
than  if  it  was  at  a  lower  altitude.  But  at  the  same  time, 
because  the  WW  is  at  a  higher  altitude,  threat  radars  will 
be  able  to  detect  the  WW  sooner  and  begin  their  attack  on 
the  WW  earlier.  The  two  tactics  chosen  for  the  experiment 
were  the  WW  leading  the  attack  force  into  the  threat  area 
by  30  seconds  and  accompanying  the  attack  force  into  the 
threat  area.  If  the  WW  leads  the  attack  force  into  the 
area,  it  may  be  possible  to  attrite  enough  radars  so  that 
more  attack  aircraft  can  get  to  the  target  area.  But  the 
MWs  will  have  increased  their  exposure  time  to  the  threats 
in  this  tactic. 
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In  the  SLAM  program  WW  altitude  is  specified  at 
line  172  (XX(3)  -  60)  and  line  178  (XX(3)  -  200).  Tactic  1, 
WWs  leading  the  attack  force  by  30  seconds  is  specified  in 
the  SLAM  program  in  Appendix  A.  This  is  denoted  by  the 
"30"  in  the  attack  aircraft  CREATE  node  entry  at  line  18 
signifying  that  the  time  of  the  first  attack  aircraft 
generated  in  the  system  is  at  30  seconds  simulation  time. 
Tactic  2  is  not  shown  in  Appendix  A  but  is  identical  to  the 
Appendix  A  program  except  that  the  attack  aircraft  CREATE 
node  entry  in  line  18  is  changed  from  "30"  to  "5"  signify¬ 
ing  that  the  first  attack  aircraft  follows  the  first  WW  in 
the  simulation  by  5  seconds. 

Experimental  Design.  Experimental  design  provides 
a  way  of  deciding  which  system  variables  to  experiment  so 
that  the  results  can  be  obtained  in  the  most  efficient  man¬ 
ner,  efficiency  equated  to  least  time  and  money. 

(In  the  terminology  of  the  experimental  design  the 
input  variables  (for  the  WW  model— altitude  and  tactic) 
are  called  factors  and  the  different  conditions  of  the  fac¬ 
tors  (WW  altitude  of  60  m  and  200  m)  are  called  levels.) 

The  design  chosen  for  the  experiment  was  a  full 
factorial  design  (Ref  8:372).  The  full  factorial  design  is 
one  in  which  all  levels  of  a  given  factor  are  combined  with 
all  levels  of  every  other  factor  in  the  experiment.  The 
advantages  of  the  factorial  design  are  as  follows: 


1.  Maximum  efficiency  in  the  estimation  of  the 
effects  of  the  variables. 

2.  Correct  identification  and  interpretation  of 
factor  interactions  if  they  exist. 

3.  The  effect  of  a  factor  is  estimated  at  several 
levels  of  other  factors,  and  thus  the  conclusions  reached 
hold  over  a  wide  range  of  conditions. 

4.  Ease  of  use  and  interpretation  (Ref  14:165). 

The  main  disadvantage  of  a  full  factorial  experi¬ 
ment  is  that  the  number  of  runs  required  to  test  all  levels 
of  all  factors  may  become  prohibitively  excessive.  The 
number  of  runs  required  for  this  experiment  did  not  pre¬ 
clude  the  use  of  the  full  factorial  design.  This  will  be 
discussed  in  Sample  Size  Determination. 

Measure  of  Merit.  The  measure  of  merit  chosen  to 
evaluate  the  experiment  was  the  number  of  aircraft  surviving 
each  run  and  reaching  the  target  area.  This  permitted  an 
analysis  of  WW  defense  suppression  effectiveness  within 
the  context  of  the  model's  threat  environments. 

Size  Determination .  A  simulation  model  is 
an  abstraction  of  a  real  system.  Experiments  are  conducted 
on  the  model  in  order  to  draw  inferences  on  the  real  system. 
Because  the  model  is  an  abstraction  or  approximation  of  the 
real  world,  the  model  must  be  executed  a  number  of  times 
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before  valid  inferences  can  be  made  on  the  output  data. 
Determining  the  number  of  runs  for  the  experiment  is  called 
sample  size  determination.  Determining  how  large  a  sample 
to  use  in  an  experiment  depends  on  the  size  of  the  risk 
you're  willing  to  take  on  the  inferences  and  variability 
that  is  present  in  the  model. 

The  sample  size  that  was  chosen  for  the  experiment 
was  determined  by  Stein's  method  (Ref  7:482).  A  trial 
experiment  of  five  simulations  was  conducted.  Each  factor 
in  the  trial  was  set  at  its  lowest  level:  WW  altitude  at 
60  and  WW  tactic  leading  the  attack  force.  The  objective 
of  the  trial  experiment  was  to  be  95  percent  confident  that 
the  sample  mean  would  be  within  one  aircraft  of  the  true 
mean.  The  following  equation  gives  the  required  number 
of  runs. 


where  M 


c 

s2 

ha/2 

"n-1 


M 


2 

s 


c 

required  number  of  runs, 

maximum  units  wrong, 

sample  variance  in  trial,  and 

t-statistic  for  (1-a)  confidence  level  with 
(n-1)  degrees  of  freedom. 
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The  trial  experiment  produced  the  following  results. 


\ 


Run 
Run  1 
Run  2 
Run  3 
Run  4 
Run  5 

Thus, 


c 


Number  of  Aircraft  Surviving 
_ (out  of  10) _ 

3 
1 
1 
5 

4 

2.76  x  3.36  _  9  2 
1 


The  results  of  the  experiment  indicated  that  the  minimum 
number  of  runs  would  be  10.  This  resulted  in  a  total  of 
five  replications  per  cell  for  the  experimental  design. 
See  Table  VII. 


TABLE  VII 

EXPERIMENTAL  DESIGN 

WW  ALT— 60  WW  ALT— 200 

xxxxx  xxxxx 

xxxxx  xxxxx 

x  *  replication. 


Tactic  1 — WW  Leading 
Tactic  2 — WW  Accompanying 
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Data  Analysis 


Once  the  factors  and  levels  for  the  experiment  were 
chosen,  the  number  of  runs  determined,  and  the  experimental 
design  set,  the  experiment  was  performed.  The  results  of 
the  experiment  are  analyzed  in  this  phase  of  the  experi¬ 
mentation. 

To  determine  if  WW  altitude  or  tactics  effected 
attack  force  survivability  a  two-way  analysis  of  variance 
(ANOVA)  statistical  procedure  was  conducted.  The  output 
of  the  statistical  analysis  along  with  the  results  from  each 
experimental  run  are  listed  in  Appendix  G. 

The  null  hypothesis  for  the  experiment  is  the 
following. 

H  :  Neither  altitude  nor  tactic  effected 
aircraft  survivability. 

The  test  statistic  for  the  experiment  is  the  F-statistic. 

The  null  hypothesis  is  rejected  with  the  F-statistic  is 
greater  than  the  following. 


F. 05, 1,16 


4.49 


The  ANOVA  shows: 


F  -  1.14 

Thus  the  null  hypothesis  cannot  be  rejected  at  the  .95 
confidence  level. 
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Although  the  results  of  the  experiment  show  an 
inconsequential  relationship  between  the  two  WW  factors  and 
attack  force  survivability,  the  experiment  must  be  viewed 
within  the  broad  perspective  of  the  developed  model. 

The  following  factors,  viewed  separately  or  in  combination 
with  one  another  may  have  contributed  to  the  experimental 
results. 

1.  Threat  environment  density.  The  number  of 
threat  radars  (78)  in  the  threat  area  was  too  great  for  a 
force  of  only  2  WW  to  suppress.  The  threat  radars  were  not 
modeled  to  stop  radiating  in  anticipation  of  an  ARM  impact, 
real  or  imagined.  Whether  radars  would  stop  radiating  in 
response  to  a  WW  attack  in  an  all-out  encounter  between 
NATO  and  Warsaw  Pact  is  open  to  conjecture.  The  simulation 
model  allows  for  this  concept  by  defining  the  threat  radar's 
ninth  attribute  as  the  radiation  attribute.  See  Appendix  E 

2.  WW  tactic.  The  WW  tactic  required  the  WW  to 
boresight  on  the  threat  before  an  ARM  could  be  launched. 

In  addition,  the  WW  attacked  every  threat  regardless  of  the 
attack  positioned  the  WW. 

3.  Ranging  routine.  Threats  were  required  to  be 
accurately  located  before  an  ARM  was  launched.  This  meant 
that  the  WW  had  to  complete  the  entire  ranging  routine 
before  ARM  launch.  Another  tactic  that  a  WW  might  use  in 
a  radar-rich  threat  environment  is  to  pre-emptively  launch 
ARMS  in  anticipation  of  threat  radar  emissions.  Although 


this  would  decrease  the  WW's  exposure  time  in  the  threat 
area  it  is  not  certain  to  what  extent,  if  any,  it  would 
decrease  the  HW's  effectiveness. 

4.  Self -protection  jamming.  The  WW  did  not  use 
self -protection  jamming  either  before  or  after  attack  by  a 
threat  radar.  If  the  WW  did  use  jamming  its  survivability 
might  be  increased. 

Needless  to  say  there  are  many  factors  that  influ¬ 
ence  and  are  critical  to  WW  defense  suppression  operations. 
A  thorough  and  comprehensive  investigation  of  these  many 
factors  is  easily  realized  with  this  WW  simulation  model. 
Manipulations  of  various  parameters  such  as  altitude,  air¬ 
speed,  timing  tactics,  radar  shutdown,  and  jamming  can  be 
accomplished  with  minor  adjustments  of  the  defining  simula¬ 
tion  variables. 

Validation 

Law  and  Kelton  (Ref  8:338)  elaborate  on  a  three- 
step  approach  to  validation  that  was  first  presented  by 
Naylor  and  Finger. 

1.  Develop  the  model  with  high  face  validity. 

2.  Test  assumptions  of  the  model  empirically. 

3.  Determine  how  representative  the  simulation 
data  Ate . 

This  three-step  approach  will  be  used  for  the  WW 
simulation  model  in  order  to  establish  its  validity. 
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Model  Face  Validity.  A  model  that  has  high  face 
validity  is  one  which  seems  reasonable  to  people  who  are 
knowledgeable  about  the  system.  To  establish  the  model's 
face  validity  the  following  was  accomplished. 

1.  Before  the  WW  model  was  developed,  experts  on 
WW  employment  from  the  37th  Tactical  Fighter  Wing  at 
George  AFB  were  interviewed.  Crew  members  and  instructors 
who  are  familiar  with  the  WW  and  its  varying  missions  were 
consulted.  In  addition,  as  the  model  was  being  built, 
these  experts  were  asked  to  confirm  various  aspects  on  pos¬ 
sible  WW  tactics  and  current  operational  concepts  to  make 
sure  that  the  model's  assumptions  were  realistic  (Ref  11). 

2.  In  order  that  the  model ' s  threat  environment  be 
represetnative  of  a  typical  enemy  threat  array,  defensive 
systems  experts  were  consulted  (Ref  4) .  The  threat  scenario 
that  resulted  developed  from  these  consultations. 

Empirical  Tests.  In  order  to  ensure  that  the  model 
behaved  as  it  was  intended  to  behave,  quantitative  tests 
were  run.  These  tests  consisted  of  the  following. 

1.  WW  Profile.  To  determine  the  validity  of  the 
WW  attack  portion  of  the  model,  a  WW  was  followed  as  it 
proceeded  through  the  threat  array  hunting  for  radars  to 
attack.  The  logic  and  action  points,  as  well  as  the  proba¬ 
bility  of  kill  routine  for  the  WW,  were  hand-calculated 
and  then  compared  to  model  output.  The  results  are  listed 


in  Appendix  I.  The  model's  output  compares  favorably  with 
the  hand-calculated  figures. 

2.  Radar-Attack  Profile.  As  with  the  WW  portion 
of  the  model,  the  logic  and  subroutines  of  the  radar- 
attack  profile  are  validated  by  following  an  attack  force 
aircraft  as  it  flies  through  the  threat  array.  The  results 
are  listed  in  Appendix  I  and  these  too  compare  favorably. 

3.  Extreme  Values.  The  model  was  tested  at 
extreme  values  of  various  system  parameters  to  determine 
if  the  model  behaved  as  it  should.  When  attack  aircraft 
velocity  was  decreased  to  50  m/sec,  all  aircraft  were 
killed.  Similarly,  when  attack  aircraft  altitude  was  set 
at  zero,  all  but  one  aircraft  survived,  the  lone  kill  being 
recorded  by  AAA.  Additionally,  WW  altitude  was  decreased 
to  zero  and  WWs  were  able  to  launch  a  total  of  seven  ARMs 
before  being  killed  by  AAA. 

Simulation  Output  Data .  The  best  test  for  a  Simula' 
tion  model  would  be  to  establish  that  the  model ' s  output 
data  closely  follows  the  output  data  one  would  expect  from 
the  real  system  being  modeled.  Because  the  developed  WW 
simulation  model  had  no  equivalent  system  with  which  to 
test  it,  a  modified  Turing  test  (Ref  8:341)  was  conducted. 
The  object  of  the  Turing  test  was  to  compare  the  model's 
output  data  against  a  hypothetical  system  by  experts  who 
are  knowledgeable  with  the  system..  To  that  end  three  WW 


experts  were  asked  to  predict  the  outcome  of  the  model's 
threat  scenario  given  the  same  input  data  as  the  model's. 
The  WW  experts  consisted  of  two  crew  members  and  a  WW 
systems  project  officer,  all  of  whom  are  familiar  with  the 
capabilities  and  concept  of  operation  of  the  F-4G  WW.  The 
predictions  of  the  panel  of  experts  agreed  consistently 
with  the  model  output  data. 

In  this  chapter  experimentation  of  the  simulation 
model  was  presented.  This  included  selecting  the  data  and 
methodology  for  the  experiment,  performing  the  experiment, 
and  analyzing  the  results.  In  addition,  the  three-step 
process  of  model  validation  was  discussed.  In  the  next 
chapter,  conclusions  and  recommendations  are  presented. 


Conclusions  and  Recommendations 


The  objective  of  this  thesis,  as  stated  in  Chapter  I 
was  to  develop  a  methodology,  through  a  simulation  model, 
for  evaluating  the  WW  defense  suppression  mission.  Conclu¬ 
sions  resulting  from  this  objective  break  down  into  two 
categories: 

1.  Primary — those  based  on  the  SLAM  modeling  of 
the  FEBA  air  battle,  and 

2.  Secondary — those  based  on  the  model's  experi¬ 
mental  results. 

Primary  Conclusions 

This  thesis  developed  a  model  for  evaluating  com¬ 
bined  air  operations  and  tactics  near  the  FEBA.  The  model 
expanded  on  Leek  and  Schmidt's  initial  work  by  including 
WW  defense  suppression  support  for  an  ingressing  fighter 
strike  force.  The  model  required  the  flexibility  to  allow 
an  experimenter  the  latitude  for  evaluating  alternative 
procedures  and  tactics  yet  retain  the  accuracy  required  to 
portray  the  actual  combat  area.  The  SLAM  model  does  cap¬ 
ture  the  dynamic  nature  of  the  FEBA  air  battle.  The  com¬ 
bined  modeling  approach  allows  the  experimenter  to  view  the 
interactions  between  the  aircraft  and  threats  on  a  moment- 
by-moment  basis.  SLAM's  flexibility  offers  an  experimenter 


the  opportunity  of  varying  the  state  and  global  variables 
to  design  meaningful  experiments  and  evaluate  effects  on 
WW  tactics  and  procedures.  Appendix  C  lists  these  variables 
and  gives  an  idea  to  the  type  of  experimental  design  that 
can  be  undertaken . 

Experiment  Conclusions 

The  model's  present  structure  was  developed  with 
the  objective  of  obtaining  the  highest  P^  per  engagement 
for  both  the  WW  and  enemy  threat  systems  independent  of 
other  action.  it  represents  the  most  restrictive  case 
for  both  the  major  system  elements  (WWs  and  threats) 
and  required  the  most  computer  logic  and  subroutines.  At 
the  other  end  of  the  scenario  spectrum  lies  the  area  of 
maximum  self -protection  for  both  elements.  It  requires 
little  logic  (stay  low  for  the  WW  and  do  not  radiate  for 
the  threat  radars) .  In  between  these  two  extremes  lies 
the  area  of  actual  operations  and  an  area  where  experimenta¬ 
tion  can  be  done:  achieve  a  desired  level  of  P^  on  a  threat 
while  achieving  some  desired  level  of  self -protection. 

Consider  the  WW.  Its  stated  objective  for  a 
defense  suppression  mission  is  to  eliminate  the  enemy's 
radar  controlled  weapon  threat.  The  existing  model  simu¬ 
lates  a  methodology  which  attempts  to  maximize  this  proba¬ 
bility  of  destroying  the  threat.  The  WW's  ranging  routine 
duration  (30  degree  of  relative  bearing  change)  represents 
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the  method  of  most  accurately  determining  the  threat's 
location  given  the  WW  must  operate  in  the  low  altitude 
arena.  The  attack  phase  requirement  for  boresighting  on 
the  threat  again  represents  the  method  required  for  achiev¬ 
ing  the  highest  P  . 

Jv 

•  The  same  can  be  said  of  the  threat's  portion  of 

i 

the  model.  The  logic  attempted  to  maximize  the  single 
t  engagement  P^.  For  both  the  SAM  and  AAA  systems,  the  shot 

was  delayed,  if  possible,  until  the  aircraft  came  abeam 
the  site,  thus  increasing  the  maximum  radar  cross-section 
while  launching  at  the  minimum  range  (resulting  in  the 
minimum  CEP  and  maximum  P^) . 

To  get  from  the  restrictive  case  of  the  existing 
model  to  a  likely  operational  area  requires  small  deletions 
from  and  changes  to  the  existing  model  whereas  going  from 
the  total  self-protection  end  of  the  spectrum  to  the  opera¬ 
tional  area  requires  major  additions  to  the  program.  For 
this  reason  the  restrictive  case  is  easier  for  an  experi¬ 
menter  to  use  and  modify.  The  experimental  conclusions  and 
1  recommendations  are  made  with  this  in  mind  and  represent 

I  areas  where  operators  in  the  field  indicate  a  need  exists 

for  additional  study  and  the  analyst  can  readily  adapt  the 
model. 

The  model's  experiment  of  varying  WW  tactics  (timing 
with  the  strike  force)  and  altitude  did  not  affect  strike 
force  survivability.  It  must  be  emphasized  that  this 
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conclusion  is  drawn  based  on  the  constraints  and  limita¬ 


tions  imposed  in  the  model.  In  addition,  the  experiment's 
measure  of  merit  was  strike  force  aircraft  reaching  the 
target  and  not  WW  survivability.  The  results  appear  con¬ 
sistent  for  an  area  where  the  ground  threats  enjoy  a  7:1 
t  numerical  advantage  over  the  aircraft. 

Specific  limitations  and  constraints  imposed  by 
the  model  on  the  WW  and  affecting  its  performance  are 
explained  below.  In  general,  the  WW  tactics  employed  in 
the  model  increased  the  WW  exposure  time  in  the  threat 
environment . 

1.  The  WW  continued  to  attack  threats  until  it 
depleted  all  weapons.  In  the  dense  threat  environment 
this  longer  exposure  time  to  the  threats  increased  the 
site's  P^  against  the  WW.  {Restrictive  case  discussed 
above. ) 

2.  The  WW  did  not  employ  an  ARM  turn  launch 
option  (off  the  bore sight  axis) .  The  increased  time 
required  to  turn  and  boresight  on  the  threat  again  increased 
the  WW' s  exposure  time.  (Again,  restrictive  case  discussed 
above.) 

3.  The  WW  did  not  immediately  abort  the  mission 
scenario  if  the  slant  range  to  the  site  was  less  than  the 

i 

minimum  ARM  range.  (Restrictive  case.) 

4.  The  WW  ARM's  against  a  site  did  not  vary 
for  launch  altitude,  although  degradation  for  lower  launch 
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altitude  does  exist.  (Restrictive  case.) 

Model  limitations  imposed  on  the  threat  radars  and 
affecting  the  experiment's  outcome  are  listed  below: 

1.  The  SAMs  did  not  have  a  minimum  launch  range. 
Actual  SAMs  have  both  a  minimum  and  maximum  launch  range. 

By  excluding  the  minimum  range,  the  site  could  allow  the 
aircraft  to  fly  over  the  site  before  launching  a  missile 
and  thus  the  model  calculated  an  inflated  based  on  this 
lower  launch  range  to  impact. 

2.  The  sites  neither  jammed  the  WW's  RHAW  receivers 
nor  stopped  radiating  if  detecting  a  possible  WW  attack. 

3.  The  threat  radars  were  modeled  as  totally  con¬ 
centrating  on  the  attacking  aircraft.  Combined  arms  opera¬ 
tions  near  the  FEBA  may  hinder  the  site ' s  operation  when 
other  ground  threats  in  the  area  attack  near  the  SAM  sites. 

4.  The  threat  network  included  neither  infrared 
(IR)  and  visual  detection  means  nor  small  arms  fire. 

Recommenda  t ion  s 

Based  on  the  experimental  conclusions,  the  following 
recommendations  are  made  (here,  the  model  moves  from  the 
restrictive  to  operational) : 

1.  When  possible,  ARM  launch  should  occur  at  the 
maximum  launch  range. 

2.  If  the  WW  finds  itself  in  a  position  inside 
the  minimum  launch  range,  the  aircraft  should  immediately 
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abort  the  attack  and  not  attempt  to  maneuver  to  achieve  a 
permissible  launch  condition. 

3.  An  effort  should  be  made  to  decrease  the  time 
required  to  complete  the  ranging  phase  of  the  mission. 

4.  WW  crews  should  always  operate  at  the  lowest 
possible  altitude. 

5.  An  analysis  should  be  continued  in  developing 
a  longer  range  ARM. 
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VI.  Recommended  Follow-on  Study 

This  thesis  concentrated  on  developing  a  model  for 
experimenting  with  combined  air  operations  near  the  FEBA 
/  and  did  not  concentrate  on  the  actual  experimentation.  The 

model  achieved  an  operational  status  and  is  now  ready  to 
•  be  used  for  experimentation.  The  following  indicate  pos- 

sible  follow-on  study  areas: 

1.  The  WW  preemptive  ARM  launch  option  should  be 
investigated.  Preemptive  launch  requires  no  mission 
ranging  phase.  Based  on  intelligence  estimates  the  crew 
launches  its  ARM  into  a  concentrated  area  of  weapons  from  a 
maximun  range. 

2.  The  use  of  WW  self-protective  ECM  and  its 
effect  on  WW  survivability  must  be  analyzed.  This  may  be 
limited  by  the  WW’s  onboard  RHAW  equipment  and  the  possible 

¥ 

interference . 

3.  The  number  of  ARM  launches  the  WW  attempts  per 
mission  may  affect  its  survivability .  For  example,  con¬ 
tinuous  attacks  in  a  dense  threat  radar  environment 
resulted  in  the  WW  being  destroyed  on  all  experiment  runs. 

By  limiting  the  attacks  to  a  smaller  number  per  mission, 
followed  by  the  WW  withdrawing  to  safe  airspace  may  increase 
the  number  of  surviving  WW. 
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4.  A  modification  of  the  model  would  allow  the 
WW  to  employ  the  ARM  launches  from  its  turn  mode,  thus 
decreasing  the  time  required  to  perform  each  mission 
scenario. 
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(ARM  launch  and  Impact,  PK  determined) 

WHO  WEASEL  1  -  ATTACK  PROFILE  NETWORK 
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Appendix  B 


SLAM  State  Variables 


Corres¬ 

ponding 


ss 

Use 

DD 

Definition  of  DD 

1-10  Attack  aurora  ft  (1-10)  x-position  ] 

L-10 

Velocity 

21 

WWL  x-position 

21 

WWL  velocity 

in  x-directicn 

22 

WW2  x-position 

22 

WW2  velocity 

in  x-directicn 

25 

WWl  y-positicn 

25 

WWl  velocity 

in  y-direction 

26 

WW2  y-position 

26 

WW2  velocity 
in  y-direction 

29 

WWL  heading 

29 

WWl  heading  rate 

change 

30 

WW2  heading 

30 

WW2  heading  rate 

change 

33 

WWL  working  designator 

— 

34 

VM2  working  designator 

— 

37 

Distance  from  WWL  to  attacked  threat 

— 

38 

Distance  from  WW2  to  attacked  threat 

— 

41 

WWl  line-of -sight  distance 

— 

42 

WW2  line-of-sight  distance 

— 

45 

WWL  relative  bearing  to  threat 

— 

46 

WW2  relative  bearing  to  threat 

— 

49 

WWL  absolute  relative  bearing  to  threat 

— 

50 

WW2  absolute  relative  bearing  to  threat 

— 

XX 


Use 


I- 2  Not  used 

3  WW  altitude 

4  Actual  firing  range  of  ABM 

5  WW  velocity 

6  Counter  for  number  of  VWs  created 

7  Time  between  creations  for  attack  force  aircraft 

8  Time  between  creations  for  W 

9  Counter  for  nuiber  of  attack  aircraft  created 

II- 12  Not  used 

13  Attack  force  altitude 

14  Attack  force  velocity 

15  Time  for  VW1  to  get  to  ABM  release  point  after  completion  of 

triangulation  route 

16  ABM  time  of  flight  for  WWl 

17  AIM  counter  for  VW1  (number  of  ABMs  launched) 

18  Not  used 

19  Time  for  WW2  to  get  to  ABM  release  point 

20  ABM  tine  of  flight  far  WW2 

21  ABM  counter  far  WW2 

22-31  Not  used 

32  x-positian  of  threat  being  attacked  by  WWl 

33  x-position  of  threat  being  attacked  by  WW2 

34-35  Not  used 

36  y-position  of  threat  being  attacked  by  WWl 

37  y-position  of  threat  being  attacked  by  WW2 

38-47  Not  used 

48  Counter  far  WWs  reaching  "heme" 

49  Flag  far  threat  killed  by  a  WW 

50  Counter  for  threats  killed  by  a  WW 

51  Not  used 

52  Standard  deviation  far  road  network 
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XX  Use 

53  Not  used 

54  Flag  for  WWl  killed  by  a  threat 

55  Flag  for  attack  aircraft  reaching  target 

56  Counter  for  attack  aircraft  reaching  target 

57  Counter  for  attack  aircraft  killed  by  threats 

58  Flag  far  VW2  killed  by  threat 

59  Counter  far  WWs  killed  by  aircraft 

60  Not  used 

61-70  y-position  of  attack  aircraft  1-10 

81  AFM  load  of  WWl 

82  ABM  load  of  WW2 

83-84  Not  used 

85  Flight  time  ABM  #1  WWl 

86  Flight  time  ABM  #2  WWl 

87  Flight  time  ABM  #3  WWl 

88  Flight  tine  ABM  #4  WWl 

89-92  Flight  time  ABM  #1  through  #4  for  WW2 
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Appendix  D 
SLAM  File  Structure 


File 


All  radars  in  threat  scenario 
Radars  that  WWs  can  attack 
WWl  working  threat  file 
WWl  kill  file — radars  that  WWl  has  killed 
WW2  working  threat  file 
WW2  kill  file — radars  that  WW2  has  killed 
Radar  aircraft  files;  file  7  contains  those  radars  engaged 
with  aircraft  1,  file  8  with  aircraft  2,  etc. 

Radar — WWl  file 
Radar— WW2  file 
Data  far  radar  cross-section 


Attribute  Listin' 


Attribute 

1 

2 


3 

4 

5 


6 

7 

8 
9 


14 


Use 

Call  sign 

m—  211,  212 

Attack — 1  to  10 
ABM  Load  (VJW  only) 

Use 

Sequential  ntxcber  (1-85) 

Type 

1—  AAA 

2—  SAM-A 

3—  SAM-B 

4—  SAM-C 

5—  -SAM-O 

6 —  EW 

x-positicn 
y-positian 
Associated  EH 
0 — None 

10 — Associated  EH  killed  by  HH 
Threat  maximum  effect  range  (meters) 

Threat  minimum  effective  altitude  (m) 

Call  sign  of  aircraft  being  attacked  by  threat 
Radiating 
0— No 
1— Yes 

Call  sign  of  HW  attacking  threat 
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Event/ 

Subroutine 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


Type 
WW  only 

WW  only 

WW  only 

WW  only 

WW  only 

WW  only 

WW  only 
WW  only 
WW  only 

WW  only 


Function  and  Cements 

Starts  WW  ranging  routine;  realized  from 
DETECT  nodes  WlS  or  W2S 

Rolls  WW  out  of  turn;  realized  frati 
DETECT  nodes  W17  or  W27  - - 

Starts  WW  bade  into  threat  it  is  attack¬ 
ing;  realized  from  DETECT  nodes  W10  or 
W20 

Reduces  WW  heading  rate  as  it  starts  to 
line  qp  for  ARM  launch;  realized  from 
DETECT  nodes  WIT  or  W2T 

Determines  ABM  firing  parameters  or  if 
too  close  to  threat  starts  repositioning 
routine;  event  realized  from  DETECT 
nodes  W1B  or  W2B 

ARM  fired  at  this  event;  ENTER  nodes  1 
or  2  are  called  to  determine  ABM  PK; 
event  realized  at  the  end  of  ABM 
release  time — XX  (15/19) 

Gets  the  WW  back  into  the  hunting 
stage  5  seconds  after  ARM  release 

Enters  WW  back  into  network  based  on 
position  and  threat  it's  to  attack 

Realized  fxcm  WW  ENTER  node  1/2; 
determines  PK  of  the  ABM;  removes  radar 
from  appropriate  file  (as  necessary) 

Repositioning  routine  event;  realized 
from  event  5  after  time  duration  of 
DETECT  node  W1R  or  W2R 
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Event/ 

Subroutine 

Type 

11 

WW  only 

12 

WW  only 

13 

WW  only 

14 

All  acft 

15  All  acft 


16  All  acft 


17  All  acft 


Function  and  Garments 

Repositioning  routine;  realized  from 
DETECT  node  W1D  or  W2D 

Starts  VM  hone  after  it  fires  its  last 
ASM;  previous  event — event  6  (ARM  launch 
event) 

Rolls  VJW  out  after  its  initial  turn 
hone;  from  DETECT  node  W1P  or  W2P 

Radar  search  event;  occurs  every  1 
second  as  long  as  the  aircraft  is  in 
the  FEBA  area  and  not  killed;  calls 
subroutine  SEARCH;  evaluates  if  aircraft 
has  reached  target  area;  allows  only  a 
maximum  of  5  threats  to  work  aircraft; 
this  is  a  network  event 

Discrete  event;  called  at  the  end  of 
tracking  and  acquisition  tine  by  function 
SCHDL  in  subroutine  SEARCH;  based  cn 
evaluation  of  PK  for  the  threat,  events 
16,  17,  or  19  are  scheduled 

Discrete  event;  called  at  time  of 
launch  (TL)  from  event  15  by  9CH3L  func¬ 
tion;  based  on  PK,  event  17  or  19  is 
scheduled 

Discrete  event;  call  at  time  of 
weapons  inpact  (TI)  from  event  15  ot  16 
by  SCHDL  function;  event  18  is  scheduled 
30  seconds  from  this  event's  time 
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Event/ 

Subroutine 


Type 

All  acft 


Function  and  Ccmnents 


18 


19 


SEARCH 


All  acft 


All  acft 


Discrete  event;  called  from  event  17; 
releases  appropriate  radars  if  aircraft 
was  killed/not  killed,  removes  aircraft 
from  network  if  killed 

Discrete  event;  called  from  event  15  or 
16  when  PK  is  below  certain  value; 
frees  radars  to  search  for  new  aircraft 
to  attack 

Discrete  event;  called  by  event  14; 
radars  are  paired  vp  with  aircraft  based 
on  correct  values  of  radar  /aircraft; 
event  15  is  scheduled  based  on  calcula¬ 
tion  of  tracking  and  acquisition  time 
for  radar 

Discrete  event;  called  from  event  15, 

16,  and/or  17;  evaluates  four  factors — 
time  of  launch  (TL) ,  time  of  inpact  (TI) , 
range  to  intercept  (RI) ,  and  PK  of 
weapon  (FKR) 


FROB 


All  acft 


SLAM  Computer  Model 
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I 


1  RBNiCN17MM.TlfiflilGiS0.  T82M24, NEWER, 4547 

2  ATTACHtPROCFILt IDsA810171>SM=ASDAD. 

3  BEGIN, NOSFllE. 

4  GETiBM.IOsMAGGIE. 

5  ATTACH* PROCFILtSLAHPROC* I0=AFIT . 

4  BEGIN, SLAN„R*BBN£,PL*14#ff. 

7  CEN.ANOERSGN  AND  H£MNER.THESIS*2/9/8Zrl0rTES» 

8  UNITS*  19.14. 5#ff 

9  INITIAU2E.I.8W. 

If  C0NTINU0DS.3Z.Zfil.5ili; 

11  network; 

12  CREATE, H(8>, I, ,2,1! 

13  ACT/1! 

14  ! WEASEL  ASSIGN  NODE 


15 

MAS 

ASSIC«.II(4ls«(41+l.ATRIB(l)-n(4)+21fi 

14 

ASSIGN.  II*II(4)+8ftATRIB(2)-U(II)i 

17 

ACT/2. i. Cl. 

18 

CREATE. H 17), 31.. 11, 1! 

19 

ACT/3! 

Zf 

{ATTACK  FORCE  ASSIGN  NODE 

Zt 

AAS 

ASSIGN.  II  (9)  ill  19) +1 .  ATRIB(l)  *11  (9)  5 

22 

ACT/4.’ 

23 

Cl 

COON, 2! 

24 

ACT,,ATRIB(1) .EO.ZlliGNl! 

25 

ACT,,ATRIB(l).EQ.212iGW2! 

21 

act.mRatk; 

27 

.'NETWORK  FOR  Mil 

28 

CN1 

coonj-i; 

29 

ACTmSS(37).LT.$S(41),RN1! 

31 

ACT,REL(W1$1! 

31 

mi 

EVENT. 1.1! 

32 

ACT.RELIU171! 

33 

REl 

EVENT,/, l! 

34 

ACT.REUWlf)! 

35 

EVENT. 3.1! 

34 

ACT.REUN1T)! 

37 

EVENT.4.1! 

38 

ACT, REL (NIB)! 

39 

FRt 

EVENT, 5.1! 

41 

ACT,Ra«HR),SS137).LT.8fff.fi«W; 

41 

ACT, 11(15! ! 

42 

EVENT ,4. 15 

43 

ACTr.l.ATRIB(Z) .EQ.9.WCH? 

44 

ACT, 5! 

45 

EVENT, 7,1! 

44 

ACT..1! 

47 

EVENT, 8,1! 

48 

ACT , REL (NIC) »SS (37) .IT .SS (41) ,RE1! 

49 

ACT,,, GUI, 

SI 

ENTER, 1,1! 

C 
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51 

ACT » XX (35) » XX ( 17) .EQ.liURK! 

52 

ACT  >  XX (86) » XX (17) .E8.2.URK! 

53 

ACTiH(87)iIX(17).EQ.3>URKi 

54 

ACT i XX (88) t XX (17) . EQ . 4 1 URK  > 

55 

m 

EVENT. 9.1J 

56 

ACT . . ATRIB ( 1 ) .EQ. XX (49) > 

57 

UKL 

ASSIGN.  !I(5I»=1X  (5l)+t ! 

58 

act; 

59 

C0LCT.H(5«)  .RADARS  UU  KILL. 

61 

act; 

61 

tern; 

62 

m 

EVENT. Ilf  1! 

63 

ACTtREL(UlO); 

64 

EVENT.ll.l.? 

65 

acti&l(uie)..fri; 

66 

UGH 

EVENT. 12. l; 

67 

ACT.REL(UIP); 

68 

EVENT. 13. li 

69 

ACT.REL(UIJ); 

7« 

m 

ASSIGN. 11(48) =11 (48) +i; 

71 

act; 

72 

HFR 

COLCT *  11 (48) > UEASEL  HONE  FREE; 

73 

act; 

74 

tern; 

75 

iNETUORK  FOR  UU2 

76 

GU2 

cooN.i; 

77 

ACT>iSS(38).LT.SS(42)iRN2; 

78 

ACT.REKU2S)  i 

79 

RN2 

EVENT. l.l! 

81 

ACT.REL(U27); 

81 

RE2 

EVENT. 2. 1? 

82 

ACT* REL(U2f ) S 

83 

EVENT. 3. If 

84 

ACT.REMU2T)! 

85 

EVENT.4.1? 

86 

ACT.REUU2B); 

87 

FR2 

EVENT, 5. 1! 

88 

ACT.REKU2R)  >SS(38)  .LT.8fNJ.UR2! 

89 

ACT .11 (191! 

9« 

EVENT.6.1! 

9t 

ACT..ATRIB(2).E0  J.UG2! 

92 

ACT .5! 

93 

EVENT.7.1! 

94 

ACT. .1! 

95 

EVENT .8. ti 

96 

ACT.REKU2C)  .SS(38)  .LT.SS(42)  .RE2! 

97 

ACT...CU2! 

98 

UR2 

EVENT.ll.l; 

99 

ACT.REUU29); 

Iff 

EVENT.ll.l! 

138 


ill 

ACTiREL(V2£)((FR2( 

1IZ 

UG2 

EVENT ( IZ» 1 » 

113 

ACT(REL(U2P)( 

104 

EVENT# 13(1! 

115 

ACT,REL(V2J),,«HR! 

1ft 

ENTER?2(1! 

117 

ACT ,  XX  (89) ,  XX  (21 )  .£Q .  1 (UKZF 

118 

ACT, 11(91), 11(21). EQ. 2, UK25 

in 

ACT, XX  (91)  «XX(Z1>.  £0.3, UK2! 

in 

ACT. 11(92) (11(21). E0.4(UK2i 

in 

UK2 

EVENT, 9,1! 

112 

ACT..  ATRIBUl.EQ.  11(49).  VKL! 

113 

ENTER, 3.1. 

114 

TERN. 

115 

(DETECT  NODES  FOR  UN’S 

116 

(DETECT  FOR  UU1 

117 

VIS 

DETECT, SS(37) .IN, SS(41)! 

118 

TERN! 

119 

U17 

DETECT, SS(49), 1,75.1, 8! 

121 

TERN! 

121 

Uli 

DETECT, SS(49), IP, 115,1! 

122 

TERR! 

123 

WIT 

DETECT, SS (49), IN, 11,1! 

124 

TERR! 

125 

U1B 

DETECT, SS(29),1,SS(45) ,8, 

126 

TERR! 

12  7 

U1R 

DETECT, SS(49)d,  179, 1! 

128 

TERR! 

129 

U10 

DETECT, SS(37),l,18Nf,f; 

131 

TERR! 

131 

VIE 

DETECT, SS(29),l,SS(45),f! 

132 

TERR! 

133 

HIP 

DETECT, SS(29),1,SS(45),I! 

134 

TERR! 

135 

V1J 

DETECT, SS(37),  1, 1H4,I! 

136 

TERR! 

137 

VIC 

DETECT, SS<49), 1,75,1! 

138 

TERR! 

139 

(DETECT  NODES  FOR  IMZ 

141 

V2S 

DETECT,  SS(38),U,SS(42),f! 

141 

TERR! 

142 

V27 

DETECT, SS(5f), 1*73, «! 

143 

TERR! 

144 

V2f 

DETECT, SS(5f), IP, 115,1! 

143 

TERR! 

146 

V2T 

DETECT,SS(SI),lN,lf,f! 

147 

TERR? 

148 

V2I 

DETECT, S$(3f),X,S$(46),i; 

149 

TERRI 

151 

V2R 

DETECT, SS(5f>, 1,179,1! 
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1S1 

TERM! 

15Z 

U20  DETECT  iSS(381iXi18M8i8> 

153 

TERN! 

154 

U2E  DETECT, SS<38M.SS<46)  ,8! 

155 

TEW! 

156 

U2P  DETECT, SS(38),X,SS<46), 8! 

157 

TERN! 

158 

U2J  DETECT iSS( 38) ,1,1888,8! 

159 

TERN! 

168 

UZC  DETECT  >  SS ( 58 ) >  X  >  75  r  0 ! 

161 

TERN! 

162 

• 

» 

163 

1 NETWORK  FOR  ALL  AIRCRAFT 

164 

• 

f 

165 

RATK  EVENTiUili 

166 

ACT i * ATRIB(l) .EQ.XX (55! tTCT 

167 

ACT » 1 » »RATK! 

168 

TCT  ASSi5N»n(56)=n(56)+lfn(5 

169 

ACT* iXX(51) .EQ.18! 

178 

TERfl, i; 

171 

ENDNETUORK! 

172 

IMTLCt  XX (7) =5»  XX (8»  =38 1 XX (3) =681 

173 

SIMULATE! 

174 

SIMULATE! 

175 

simulate; 

176 

SIMULATE! 

177 

simulate; 

178 

INTLCiXX(3) =288! 

179 

SIMULATE! 

188 

SIMULATE! 

181 

SIMULATE! 

182 

SIMULATE! 

183 

SIMULATE! 

184 

FIN! 

'4 

1 
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t 

PROGRAM  MAIN ( INPUT i OUTPUT t T APE5= INPUT t TAPE6=0UTPUT i TAPE7) 

2 

DIMENSION  NSET (25818) 

3 

COMMON  QSET125008) 

♦ 

COMM/SCIMl/ATRlBdff)  >00(188)  iDDL(18f)iDTN0UiIIiHFA>HST0P 

5 

ltNCROR>MPRNT>IWRUN*l0tSET»NTAPE>SS(100)  >SSL(108)  >TNEIT> 

6 

tTN0H>U(188) 

7 

EQUIVALENCE (NSET ( 1 ) iQSET ( 1 ) ) 

8 

MCR0R=5 

9 

NPRHT-6 

If 

NTAP£*7 

11 

NNS£T=25080 

12 

CALL  SLAH 

13 

STOP 

14 

END 

IS 

SUBROUTINE  STATE 

16 

COmON/SCOHl  /  ATRIB  ( Iff)  >00(188)  iDDLdff)  >DTNOU*II>NFA.RSTOPi 

17 

t>NCR0R>NPRNTiNNRUN>NNSETiNTAPE>SS(i88) * SSL ( 188) iTNEXTi 

18 

(TNOUtHIlBB) 

19 

COHMON/UCOM1 /RATE (2) >L(2) 

2f 

DIMENSION  A ( 16) 

21 

C 

22 

c 

RATE  EQUATIONS  FOR  STRIKE  AIRCRAFT 

23 

c 

24 

JJ=XX(9) 

25 

IF(JJ.EQ.8)C0  TO  15 

26 

00  If  I*l>JJ 

27 

08(11*11(14) 

28 

11 

CONTINUE 

29 

c 

3f 

c 

RATE  EQUATIONS  FOR  IM  AIRCRAFT 

31 

c 

32 

15 

KK:H(6)+2f 

33 

IF (KK.EQ.2I) RETURN 

34 

00  28  I-Z1>KK 

35 

J»I-Z8 

36 

C 

37 

C 

HEADING  RATE  CHANCE  IS  EQUAL  TO  RATE 

38 

C 

39 

D0(l*8MRATE(J) 

41 

IF(SS(I+8).LT.8)SS(I+8)*SS(I+8)+368 

41 

IF($S(I+8I .CT.368 ) SS( 1+8) *SS( 1+8) -368 

42 

c 

43 

c 

VELOCITT  OF  UH  IN  I-OIRECTION  AND  T-DIRECTION 

44 

c 

45 

D ( I ) *XX (5) «COSO (98-SS ( 1+8) ) 

44 

BD(I+4)*H(5)*SIND(98-SS(I+8) ) 

47 

c 

48 

c 

BASED  ON  TIE  VALUE  OF  THE  UH  WORKING  OESICNATORi 

49 

c 

EITHER  COMPUTE  A  "LOOK  AHEAD"  RADAR  TO  ATTACK 

5« 

c 

OR  PROCEED  UITH  OTHER  STATE  VARIABLE  COMPUTATIONS. 
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c 


51 


52 

IF(S$(I+12).CT.I.AND.SStI+12).LT.9)C0  TO  111 

53 

IF(S$(I+12).Efl.9)TH£H 

54 

IKI+llhf.l 

55 

XX(I+l5)=-l#flf.l 

54 

CO  TO  111 

57 

ENDIF 

58 

C 

59 

C 

PRE-SELECT  RADAR  FOR  W  TO  ATTACK  BASED  ON  THE 

4« 

c 

CLOSEST  RADAR  TO  THE  UU. 

41 

c 

42 

SR2 -999999999. 9 

43 

DO  22  18=1.1*0(2) 

44 

CALL  COPT (10.2. A) 

45 

IF  (A  ( 14 )  .CT.DCO  TO  22 

44 

X=A(3) 

47 

7=A<4> 

48 

SR1=S8RT(<SS(I)-X)m2+<SS(I+4)-T)h2+XX(3)«2) 

49 

IF(SR1.LT.SR2)THEN 

71 

SR2-SR1 

71 

L(J)--A(1) 

72 

xxu+n)*x 

73 

mi+tsw 

74 

ENDIF 

75 

22 

CONTINUE 

74 

111 

SSII+14J=SttT((S$(I)-IX(I+li))«2+<$S(I+4)-niI+15» 

77 

t«2+XX(3J«21 

78 

SSa+2f)=4tl7.2+S8RT(XX(3>> 

79 

0*  (XI  (1+15)  -SS  ( 1+4)  >  /  (XU  t+1 1 )  -$$  ( I ) ) 

8f 

IF(D.CE.I.AND.(II(I+15)-SSn+4)).CE.I)T1€N 

81 

SS(I+24)*57.3+ATAN(l/B) 

82 

ELSEIFtD.LT.f  .AND.  (IX<I+15)-SS<I+4) )  .CE.DTHEN 

83 

SSII+24) =278+57 . 3+ABS (ATAN(D) ) 

84 

ELSEIFtD.LT.f. AW).  (H(I+15)-$S(l+4)).LT.f)  THEN 

85 

S$(I+24)=9#+57.3»A8S(ATAN(B)) 

84 

ELSE 

87 

SS(I+24)*27I-57.3«ATAN(B) 

88 

ENDIF 

89 

SS ( 1+28) * ADS ($S ( 1+8) -SS ( 1+24) ) 

91 

IF (SS(I+28) .CT . 18f ) SS ( I+Z8) *34f-SS( 1+28) 

91 

21 

CONTINUE 

92 

RETURN 

93 

END 
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94 

SUBROUTINE  INTLCO 

95 

COWtON/SCOHl/flTRIB(tli).DD(lM)  iDDLllll)  iDTNOU* II »HFAiHSTOP»MCLN 

94 

4NCRDR  •  NPRMT  t  NNRUN  ,  NNSET  •  NTAPEi  SS  ( 1M)  >  SSL  ( 1M )  iTNEXT » 

97 

(INOUiU(lM) 

98 

DIMENSION  RNC(12}»A(141, ARY (37,6) 

99 

COtWON/UCOHl /RATE  (21  >L(2) 

Iff 

DATA  A/ 14*1.1/ 

111 

DATA  RNC/12»3IHI.I/ 

112 

C 

113 

C 

DATA  IN  "ART"  IS  USED  TO  CALCULATE  THE  RADAR  CROSS 

114 

c 

SECTION  OF  THE  AIRCRAFT  BASED  ON  THE  ASPECT  THE  RADAR 

115 

c 

VIEUS  THE  AIRCRAFT.  “ART"  IS  PLACED  IN  FILE  19. 

116 

c 

117 

DATA((ARY(lA»JA)»UA=l,61,lAsl»37)/l»-2l.8,-19.3,-26.4»-26.4,-21. 

118 

45,-23.3, -25. 3, -29.4, -29.6,-23. 3, 11, -28. 1,-25. 4, -27.5, -27.5,-28.1 

119 

415,-26.5,-32.3, -27., -27., -26.5, 21, -19, 2,-25. 2, -28. 3, -28. 3, -19. 2, 

ill 

425,-26.1, -27.3, -27. 5, -27. 5, -26. 1,31,-24. 8,-28. 8, -19. 4,-19. 4, -25. 

111 

435,-31. 2,-31. 7, -22. 1,-22. 1,-31. 2, 41, -28.3, -27. 9, -26.6,-26. 6, -28. 

112 

445, -23.7, -25.3, -26.3,-24.3,-23.7, 51, -27., -28.3,-31. ,-31. ,-26.9, 

113 

455.-29. ,-28.9.-31.7, -31,7.-29. i6l»-28. 1.-34., -31., -31., -28. 1, 

114 

465, -29.8, -29.6, -31., -31., -29.8, 71, -21.8,-21.8, -21.6, -21.6, -21.8, 

115 

475, -14.6, -16.3, -15.5, -15.5, -16.6, 81, -13.3, -14.»-13.7, -13.7, -13.3 

116 

485, -13.9, -16.6, -14.1, -14.1, -13.9, 91, -5.6, -4.2, -5., -5., -5.4, 

117 

495,-11.8, -5.1, -11., -11., -1I.8.1H, -13.4, -13., -14.3, -14.3, -13.4, 

118 

4115, -24.2, -23., -21.8, -21.8, -24.2, 111, -21.5, -21.2, -21.3, -21.3, -21 

119 

4115, -23.8, -28.3, -26.4, -26.4, -23.8, 121, -22., -31.8, -27.7, -27.7, -22 

121 

4125, -25.7, -31.2, -27., -27., -25.7, 131, -24.5, -32.3, -27., -27., -24.5, 

121 

4135, -25.7, -29.9, -31., -31. ,-25.7, Ml, -26.1, -31.3, -32.6, -32.6, -26. 

122 

4145, -24, -31.9, -29.5, -29.5, -24., 15l,-24.8, -27.2, -29.5, -29.5, -24., 

123 

4155, -23.1, -23.5, -29.6, -29.6, -23.1, 161, -27., -31., -29.4, -29.4, -27. 

124 

4165, -25.3, -29.5, -24., -24., -25.3, 171, -21.1, -25.9, -23.5, -23.5, -21. 

125 

4175, -14.5, -25.4, -21. ,-21., -16.5, 1M, -15., -29., -21., -21. ,-15./ 

126 

DO  47  IM=  l » 37 

127 

A(1)=ART(!N,1) 

128 

A(Z)*ART(IH«2)+3I.I 

129 

A(3)*ART(IN,3)+3I.I 

131 

AM)  *ARY  (IN, 41+31.1 

131 

A(5)*ARY  (IN,  51+31.1 

132 

A(4)*ARTUN» 4J+3I.I 

133 

CALL  FILEH(19,A) 

134 

47 

CONTINUE 

135 

C 

136 

C 

INITIAL  VALUES  FOR  ATTACK  AIRCRAFT 

137 

C 

138 

DO  21  1*61,71 

139 

!!(I)  *UNFW(245H.I,255M.1, 11 

14« 

JJ*I-6I 

141 

SS(JJ1*I.I 

142 

21 

CONTINUE 

143 

C 

143 


s**  V 


144 

C 

INITIAL  VALUES  FOR  Ml 

145 

C 

144 

DO  22  1=21 »22 

147 

SS(I1=8.I 

148 

SS(I+4hUNFRH(245M.§.255M.M) 

149 

SS(I+8)=9#.I 

151 

SS (I+12J  =#.• 

151 

SS( 1+16) =1.1 

152 

SS(I+2f)=f.f 

153 

SS  (1+24) *1.1 

154 

SS(I+28)=8.0 

155 

00(1+41=1.8 

156 

DD(I+8)=0.0 

157 

22 

CONTINUE 

158 

C 

159 

C 

SET  VALUES  FOR  OTHER  VARIABLES 

161 

C 

161 

II (52) =4080.0 

162 

U (51 *247.1 

163 

11(131=311.1 

164 

H (14) *247.1 

165 

n<81)*22.» 

166 

11(821*22.1 

167 

11(83) *22.8 

168 

11(84) *22.1 

169 

RATE<1)=0. 

17* 

RATE (2) =1. 

171 

C 

172 

C 

CONSTRUCTION  OF  RADARS 

173 

C 

174 

DO  55  I=l»85 

175 

A(1)*I 

176 

IF (A(  1)  .CE. 49. AND. A(l)  .LE.5DC0  TO  ! 

177 

IF(A(l).CE.52.AND.A(l).LE.6f)C0  TO 

178 

IF(A(1».CE.61.ANB.A(1).LE.70>CO  TO 

179 

IF(A(1) .C£.7l.AND.A(l) ,LE.77)C0  TO 

181 

IF(A(1).CE.78)C0  TO  62 

181 

C 

182 

C 

TYPE  AAA  RADARS 

183 

C 

184 

A(2)=l.« 

185 

A(6) *2998.1 

186 

A(7)*f.« 

187 

A(9)*l.l 

188 

IF(A(1).L£.12)THEN 

189 

A(3)*4«f.f 

191 

A(4)=RN0RH(RNC(1) .11(52)  .1) 

191 

ELSEIF(A(1) .LE.24)TMEM 

192 

A(3)*15N.f 

193 

A(4)*RN0RH(RNC(2)t 11(52) >1) 
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194 

ELSEIF(A(1).LE.34)THEN 

195 

A(3)-Z5M.I 

194 

A(4):RN0M(RHC(3)'H(52)«l> 

197 

ELSE 

198 

A(3)*34#i.f 

199 

A(4)-RN0RH(RNG(4)iH(52M) 

2H 

EM)  If 

211 

CO  TO  43 

212 

C 

213 

C 

TYPE  SAH-A  RADARS 

214 

C 

215 

57 

A(2)*2.l 

m 

A(6)--5##5».f 

217 

A(7)*91.l 

218 

A(9)*l.l 

219 

IF(A(1)  .LE.5IITHEN 

210 

A (3) =45000.0 

211 

A(4) =RNORR(RMC< 10) . 11(52) .1) 

212 

ELSE 

213 

A(3)=80000.0 

214 

A(4) sRMJRH (RNC( 12) ill (52) » 1 ) 

215 

ENDIF 

214 

CO  TO  43 

217 

C 

218 

C 

TYPE  SAH-B  RADARS 

219 

C 

22f 

58 

A 12) =3.1 

221 

A<4) =74150.0 

222 

A(7) s305.f 

223 

4(91*1.1 

224 

IF(A(1).LE.54)THEN 

225 

A(3)MMM.I 

224 

A(4!=RHORfl(RAC(4) >11(52) tl) 

227 

ELSE 

228 

A(3)*25W.I 

229 

A(4)>RA0Rf1(RAC(8)tll(52)fl) 

230 

EXOIF 

231 

» 

CO  TO  43 

232 

C 

233 

c 

type  sah-c  mm 

234 

c 

235 

59 

4(2) *4.1 

234 

A (4) >22251.0 

237 

4(71*15.1 

238 

4(91*1.1 

239 

IF(A(1).LE.45)THEN 

241 

4(3)*25N.I 

241 

4(4)*RMM(RNC(3)tlX(52)fl) 

242 

ELSEIF (4(1). LE.48) THEN 

243 

4(3) *5400.0 
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244 

A (4) -RNORH (RMG  45) t  XX (52) 1 1 

245 

ELSE 

244 

A(3) =15009.0 

247 

A(4)*RN0RH(RNC(7).XX<52)i: 

248 

END  IF 

249 

CO  TO  43 

234 

C 

251 

C 

TYPE  SAH-D  RADARS 

252 

c 

253 

44 

A<2)*5.4 

254 

A(4)*14244.4 

255 

A(7)*45.4 

254 

A<9) =1 .4 

257 

IF(Ad).L£.72)TKEN 

258 

A<3) *2544.4 

259 

A  (4)  *RHORN  (RNC  (3)  »XX  (52) » : 

244 

ELSEIF (A(l) .LE.75) THEM 

241 

A(3) =35444.4 

242 

A(4)*RN0f»(RNC(9)  »XX(52) 

243 

ELSE 

244 

A(3) *44444.4 

245 

A(4)=RN0Rll(RNCdi).XX(52) 

244 

ENDIF 

247 

CO  TO  43 

248 

C 

249 

c 

EU/CCI  RADARS 

271 

c 

271 

42 

4(21*4.4 

272 

A(9)*1.4 

273 

IF(A(1).LE.79)TNEN 

274 

A<3) -5444.4 

275 

A (4) = RNORH (RNC (5) • XX (52) » 

274 

A (51*1.4 

277 

IF(A( 1 ) .£0. 79) THEN 

278 

A(4)=A(4)-14444.4 

279 

A(5)*2.4 

281 

ENDIF 

281 

ELSEIF(Ad) .LE. 84) THEN 

282 

A (31*14444. 4 

283 

A(4)*RMRH(RNC(4)fIX(52)i 

284 

A(5)*3.4 

285 

ELSEIF (A(l) .LE.82) THEN 

284 

A(3)*234M.« 

287 

A{4)*RN0RJKRNC(8)  ill  (52)  • 

288 

• 

A (5) *4.4 

289 

IF(A<l).Efl.82)THEN 

294 

A(4)*A(4)-14444.4 

291 

ACS) *5.4 

292 

ENDIF 

293 

ELSEIF (Ad)  .LE.841THEN 
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294 

A(3):35000.0 

295 

A<4):RN0M(RNC(8)  »H(52)  .1I+1000.0 

296 

A(5) =6.1 

297 

IF(A(1).E8.84)THEM 

298 

A(4)sA(4) -10000.0 

299 

A(5) -7.0 

3#f 

ENDIF 

311 

ELSE 

312 

A(3)M5M0.0 

313 

A  (4)  =RH0ftf1  (RNC  ( 10) » XX  (52)  1 1 J 

314 

A  (5) --8.0 

305 

ENDIF  — 

316 

A(3)sA(3)+500#0.0 

317 

CO  TO  64 

jVO 

63 

A(3):A(3)+50000.0 

3f9 

A (5) =0.1 

311 

IF (A<3) .LE.55000) THEN 

311 

IF(A(4).CT.RNC(5))THEN 

312 

IF  (DRAIGJ  (2)  .CT.  .5)  A(5) =1.0 

313 

ELSE1FIAI4) .LE.RNCI5) ) THEN 

314 

IF (BRAND (2) .CT. .5) A(5)=2.0 

315 

ENDIF 

316 

ELSEIF (DRAND(2) .CT . .5) THEN 

317 

A<5)*8.0 

318 

ENDIF 

319 

IF(A(2).EQ.3)THEN 

321 

IF(A(3) .Efl.6ffff.fi  THEN 

321 

A(S)*3.f 

322 

ELSCIF (AC4) .CT.RNC(8) ) THEM 

323 

A(5M4.f 

324 

ELSE 

325 

A(5)*5.f 

326 

ENDIF 

327 

ENDIF 

328 

64 

CONTINUE 

329 

CALL  FILEH(liA) 

331 

55 

CONTINUE 

331 

C 

332 

C 

IN  RADAR  FILE  *2"  -  ONLY  THOSE  RADARS  NEETINC  CERTAIN 

333 

0 

SPECIFICATIONS  ARE  PLACED  IN  THIS  FILE.  INI  CAN  ONLY 

334 

C 

ATTACK  RADARS  FROM  THIS  FILE. 

335 

C 

336 

NECT*HHFE(1) 

337 

11 

!F(NEIT.EQ.0>CO  TO  21 

338 

CALL  COPY (-NEXT. l.A) 

339 

IF! (A(2) .EQ.4.0R.AC2) .E0.5.0R.AI2) .EB.6) .AN0.AI4) .LT.5D000ITHEN 

341 

CALL  FILEN(2iA) 

341 

ENDIF 

342 

NEXT*NSUCR(NEXT) 

343 

CO  TO  If 

O' 
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344 

21  CONTINUE 

345 

CALL  PRNTFU) 

344 

CALL  PRNTFI2) 

347 

RETURN 

348 

END 

349 

SUBROUTINE  EVENT (II) 

358 

COHHON/SCOMl/ATRIB(188)  *00(188)  *  DDL  (180)  * DTNGU . 1 1 . MFA » HSTOP . NCLN 

351 

8NCRDR»NPRNT.NNRUN>NNSETiNTAPE.SSd§8)  >SSL(108)  .TNEXT, 

352 

(TNOUrXX (100) 

353 

DIMENSION  REV(4)»T0F(4)»PK(4) 

354 

C0NH0N/UC0M2/TL i TI i PKR iRI « SRR 

355 

COMNON/UCOH1/RATE(Z) *L(2) 

354 

DIMENSION  A (14) >C(14) ,F( 14) .6(14) >EA(14) ,EB(14) 

357 

I-ATRIB(l) 

358 

IF(I.CT .288) THEN 

359 

1*1-198 

348 

J*I-28 

341 

jy*(2tj)+l 

342 

IR= (4*J)+11 

343 

LF--I-4 

344 

END  IF 

345 

IF(I.LE.18)LF*I+4 

344 

C 

347 

C 

IF  A  RADAR  OR  AN  AIRCRAFT  IS  KILLED  IN  THE 

348 

c 

IMMEDIATE  PRIOR  EVENT.  EVENT  NODE  "3"  IS  CALLED 

349 

c 

TO  ENSURE  THE  ENTITY  IS  REMOVED  FROM  THE  PROGRAM. 

378 

c 

371 

IFdX.LT.  14ITHEN 

372 

IF(ATRIB(1)  .E0.IK54)  .OR.ATRIB(l)  .EQ.XX  (58) )  THEN 

373 

CALL  ENTER(3.ATRIB) 

374 

RETURN 

375 

374 

377 

378 

379 

EMBIF 

ENDIF 

GO  T0(1.2.3.4.5.4.7.8.9.18.tl. 12.13.14. 15. 

814.17. 18. 19). II 

388 

c 

381 

c 

EVENT  KIN  NODE  ONLY)  •  REALIZED  UHEN  A  MM'S  RANGE  TO 

382 

c 

RADAR  DECREASED  TO  DM  FIELD-OF-VIEM. 

383 

c 

384 

l 

PRINT*. 'EV  1  TNOH*  ».TNOM.»  ATRIB(l)*  ’.ATRIB(l) 

385 

c 

BASED  ON  POSITION  AND  HEADING  OF  MM  TURN  THE  MM  TO 

384 

c 

START  RANGING  ROUTINE. 

387 

c 

388 

IF($S(I+24) . CT . 8 . AND . SS ( I +24 ) .LT. 188) THEN 

389 

IF(SS(I*8) .CT.SS(I*24> .AND.SS(I*8) .LT. (SS(I*24)* 

391 

8188)1  THEN 

391 

REV(J)*1.8 

392 

IF(SS(I+28) .LT. 75. 8) THEN 

393 

RATE  (J)  *4.8 

394 

ELSE 

395 

RATE(J)*-4.8 

394 

EMBIF 

397 

ELSE 

398 

REV(J)*2.8 
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399 


IF<SS(I+28)  .LT.75.0THEN 
RATE(J) =-♦.# 

ELSE 

RATE(J)*4.I 
END  IF 
ENOIF 
ELSE 

IF(SS(I+8).CT.(SS{I+24M8I),ANB.SSCI+8).LT. 
tSS  ( 14-24) )  THEN 
REV (J) *2.1 

IF (SS (1+28) .LT. 75) THEN 
RATE(J)*-4.I 
ELSE 

RATE (J) *4.1 
END  IF 
ELSE 

REV(J)*l.f 

IF(SS(I+28).LT.75.i)THEN 
RATE U) *4.1 
ELSE 

RATE(J)*-4.I 

ENOIF 

ENDIF 

ENOIF 

C  SET  W  UORKING  DESIGNATOR  TO  1 
SS(I+12)*l.l 
C 

C  FOR  THE  RADAR  THAT  THE  W  IS  ATTACKING.  SET  ITS 
C  14TH  ATTRIBUTE  EQUAL  TO  THE  CALL  SIGN  OF  THE  UN. 

C 

AL*L(J) 

NCET=NFIND(1.2>lif*ALil,f) 

CALL  RH0VE(NGET.2.A) 

A(14)*ATRIB(1) 

CALL  FILEH(2.Ai 
CALL  FILEH(JU.A) 

RETURN 

C 

C  EVENT  2(UU  NODE  OUT)  •  REALIZED  WEN  IN  DETECTS  71  DECREES 

C  RELATIVE  BEARING  ON  THE  RADAR  SITE  IT  IS  ATTACKING.  THIS 

C  EVENT  ROLLS  THE  IN  OUT  OF  ITS  TURN. 

C 

2  PRINT* »'EV  2  TNOU*  '.THOU.1  ATRIB(l)*  NATRIB(l) 

RATEtJHf.f 

RETURN 

c 

C  EVENT  3<IN  NOSE  OUT)  •  REALIZED  WEN  IN  DETECTS  115 

C  DECREES  RELATIVE  BEARING  ON  TIC  RADAR  IT  IS  HORNING. 


449 

450 

451 

452 

453 

454 

455 
454 

457 

458 

459 
441 

441 

442 

443 

444 

445 
444 

447 

448 

449 
471 

471 

472 

473 

474 

475 
474 

477 

478 

479 
481 

481 

482 

483 

484 

485 
484 

487 

488 

489 
491 

491 

492 

493 

494 

495 
494 

497 

498 


C  THIS  EVENT  ROLLS  THE  DU  BACK  INTO  THE  SITE. 

C 

3  PRINT»r'EV  3  THOW=  MNOU»»  ATRIB(l)-  'iATRIB(l) 

IF(REV(J) .EQ. 1 ) RATE (J) =-4.i 

IF{R£V(J) .  £Q .  Z )  RATE  ( J )  =4 .  f 
RETURN 

WWHHHWHHWWWWW<W<HW4HHHH<<HIHHWW 

C 

C  EVENT  4CUU  NODE  ONLTI  -  REALIZED  WHEN  UU  DETECTS  10  DECREES 
C  RELATIVE  BEARING  ON  THE  SITE  IT  IS  ATTACKING.  THIS  EVENT 
C  REDUCES  THE  TURN  RATE  OF  THE  UU. 

C 

4  PRINTfi'EV  4  TNOU=  1 »TNOU» 1  ATRIB(l)-  'iATRlB(l) 
IF(REV(J).EQ.l)RATE(J)-*2.0 

IF (REV ( J) . £0.2) RATE (J) =2. I 
RETURN 

C 

C  EVENT  5(UU  NODE  ONLY)  •  REALIZED  UHEN  W  IS  BORESICHTED  ON 

C  THE  RAOAR  IT  IS  ATTACKING.  BASED  ON  THE  POSITION  OF  THE 

C  UHi  A  DECISION  IS  HADE  OF  IIHETHER  TO  FIRE  AN  ARH.  IF 

C  THE  UU  DECIDES  IT  HILL  FIRE  AN  ARH <  THEN  THE  PK  OF  THE 

C  ARH  IS  DETERMINED  AND  THE  FIRING  RANGE  OF  THE  ARH  CAL- 

C  CULATED. 

C 

5  PRINT*i'EV  5  TNOH>  STNOHi'  ATRIBd!  =  ’iATRIB(l) 

C  IF  IM  HAS  ACH-78Si  THEN  SELECT  THAT  TYPE  ARH  TO 
C  FIRE (ATRIB C2>  GREATER  THAN  II).  OTHERUISE  SELECT 

C  AN  ACM-45 

IF(ATRIB(2).CE.10)CO  TO  40 
C  RH45*HINIHUH  LAUNCH  RANGE  OF  AGH-45 
RH45>8000.0 

C  XX(ll>»HAXIHUfl  LAUNCH  RANGE.  IF  RANGE  TO  TARGET  IS 

C  GREATER  THAN  1500fH<  THEN  USE  15000H  AS  HAXIHUH  RANGE 
C  FOR  THE  MISSILE.  IF  THE  ACTUAL  RANGE  TO  THE  TARGET 

C  IS  LESS  THAN  THE  HININUH  RANGE  OF  THE  ARM*  THEM  DO 

C  NOT  CONFUTE  PK!  GO  TO  "REPOSITIONING  ROUTINE." 

XXU0)*SS(I+14) 

IF(XX(li).GT. 15000. l)XX(ll) >15000. I 
IF(S$U+14).LT.RH45)C0  TO  42 

C  XX(4)*ACTUAL  FIRING  RANGE  OF  ARK.  TOF  IS  FLIGHT  TIME  OF 
C  THE  MISSILE  FROH  THE  LAUNCH  POINT  TO  IMPACT  POINT. 
XX(4)*UNFRH(8H0.lfII(10)fl) 

T0F(J1«XI(*J 7350.0 
ATRII(2)*ATRIB(2)-1 
C  DETERMINE  PK  OF  ARM 

SAHPLE«DRAND(l) 

IF(SAMPLE.LT..8)THEN 

PK(J)>1.0 

ELSE 
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499 

PK(J):f.« 

SM 

END  IF 

511 

CO  TO  41 

512 

C 

CALCULATIONS  FOR  ACH-78 

513 

41 

RH78*8###J 

514 

n<l#)=SS(I+l6l 

5f5 

IF<n(lf).CT.25i#f.f)H(tf):25fW.I 

514 

IF (SS ( 1+16) .LT.RH78) GO  TO  42 

5#7 

XX (4) =UNFRH(8«#f .#>XX (10) i 1) 

5f8 

T0F(JhXX(4)/45l.l 

519 

ATR!B(2) SATRIB(2) -11.1 

511 

SAHPLE=DRAND(l) 

511 

IF (SAHPLE.LT. .85) THEN 

512 

PX(J)*1.I 

513 

ELSE 

514 

PK(J)sf .1 

515 

ENDIF 

516 

41 

CONTINUE 

517 

C 

RLMMTIHE  FOR  THE  Mi  TO  GET  FROM  ITS  PRESENT  POSITION 

518 

C 

TO  ITS  ARM  FIRING  RANGE. 

519 

RLMis((SS(I+t6)*COS(ASIN(IX(3)/SS(I+16)  > ) -XX  (4) ) /XX  (5) ) 

521 

n(IR)<RLW 

521 

H(IR4l)*TOF(J) 

522 

RATE (J) *1.1 

523 

REV(J) =#.l 

524 

RETURN 

525 

C 

IF  THE  Mi  HAS  TOO  CLOSE  TO  FIRE  THE  ARH  THEN  START  UH 

526 

C 

TURNING  FOR  "REPOSITIONING  ROUTINE." 

527 

528 

529 

53# 

42 

RATE(J)*-4.f 

IF (REV(J) .EQ.l)RATE(J)=4.i 

RETURN 

531 

C 

532 

C 

EVENT  6 (Ml  NODE  OUT)  •  REALIZED  UHEN  THE  Mi  GETS  TO  ARN 

533 

C 

FIRING  RANGE.  "ENTER  NODE  1/2"  CALLED  FROH  THIS  EVENT  IN 

534 

C 

ORDER  TO  EVALUATE  THE  PK  OF  THE  ARN  AT  IMPACT  TINE. 

535 

C 

536 

6 

PRINT*t'EV  6  TNOH>  MNOH.'  ATRIKl)*  ’iATRIB(l) 

537 

IFU.EQ.DTHEN 

538 

n(17)>n(17)«l.f 

539 

P*XZ(17)+84 

S4fl 

II(P)*«(IR*1) 

541 

ELSEIF (J.EQ.2) THEN 

542 

II(21)»H(21)*1.I 

543 

P*II(21)*89 

544 

W(P)»H(IR*1) 

545 

ENDIF 

546 

C 

FIND  THE  RADAR  IN  FILE  2  THAT  THE  Ml  IS  AnACXINC. 

547 

c 

CALL  ENTER  NODE  1/2  WITH  THE  ATTRIBUTES  OF  THE  RADAR 

548 

c 

LOADED  INTO  THE  ATRIB  ARRT. 
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AL*L(J) 

NCTsNFINBU»2»li0.AL»0i 
IFCJ.Ea.UTHEN 
CALL  COPY (NCTt2fEA) 

CALL  ENTER(ltEA) 

ELS£IF(J.EQ.Z)THEN 
CALL  C0PT(NCT.2.EB) 

CALL  ENTER<Z»EB) 

EMDIF 

CHECK  TO  SEE  IF  UU  OUT  OF  ARKS.  IF  IT  IS  THEN  SET 
UU  WORKING  DESIGNATOR  EQUAL  TO  9. 
IF(ATRIB(2).EQ.0)SS(I+1Z)=9.0 
RETURN 


EVENT  7 (UU  NODE  ONLY)  -  THIS  EVENT  OCCURS  5  SECONDS  AFTER 
THE  UU  LAUNCHES  ARH.  UU  WORKING  DESIGNATOR  IS  SET  TO  0 
TO  ALLOW  TIE  UU  TO  BEGIN  SEARCHING  FOR  A  NEW  RADAR  TO 
ATTACK. 


PRINT*. »EV  7  TNQU= 

SS(I+12)*0.0 

RETURN 


ATRIB(1)=  ’iATRIB(l) 


EVENT  8(UU  NODE  ONLY)  -  THIS  EVENT  ENTERS  THE  UU  BACK  INTO 
THE  NETWORK  BASED  OR  ITS  POSITION  RELATIVE  TO  THE  RADAR 
SITE  IT  IS  ATTACKING. 

PRINT*, 'EV  8  TNOU=  '.TNQU,'  ATRIB(l)*  '.ATRIBdl 
IF(SS(I+14).LT.SS(I+Z0))THEN 
AL*L(J) 

NCET*NFIND(1,2,1,0,AL,0) 

CALL  RHQVEOtCET.Z.A) 

A(14)=ATRIB(1) 

CALL  FILER(Z.A) 

CALL  FILER(JUiA) 

ssa+iz)»i.f 

IF(SS(I+24).CT.0.AND.SS(I+Z4).LT.180)THEN 
IF (SSI 1*8) .CT.SS(I+Z4) .AN0.SS<I+8) .LT.(SS(I+Z4)+ 

8180)1 THEN 

REV(J)*1.0 

IF(SS(I+Z8) .LT. 75) THEN 
RATE (J) *4.0 
ELSE 

RATEM)*-4.0 

EMDIF 

ELSE 

REV(J)*2.0 

IF (SSI 1+28) .LT.75)THEN 
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RATE(J) 

ELSE 

RATE(J)  =4.1 
ENOIF 
END  IF 
ENDIF 
ELSE 

IF  (SSI  1+8)  .CT.  (SS(I+Z4) -18#)  .AND.SS(I+8)  .LT.SSU+Z4) )  THEN 
REV (J) =2.1 

IF(SS(I*Z8).LT.75)THEN 

RATE(J)--4.f 

ELSE 

RATE(J)=4.I 

ENOIF 

ELSE 

REV(J)=l.fl 

IF(SS(I+Z8) .LT.751THEN 
RATE (J) =4.1 

ass 

RATE(J)=-4.I 

ENDIF 

ENOIF 

ENOIF 

RETURN 

c 

C  EVENT  9<HU  NODE  ONLT)  -  EVENT  REALIZES  FROM  UU  ENTER 

C  NODE.  DETERMINES  PK  OF  ARM  AT  THE  END  OF  ARM'S  FLIGHT 

C  TIME. 

C 

9  PRINT*. »EV  9  THOM*  '.THOU.'  ATRIi(l)=  ATRIB(l) 
LF*ATRI»(14I-194 
J=ATRIB(14)*Z1I 
JU*(Z*J)+1 

PRINT*. 'ATU4I*  f . ATRIB(  14) . fLF=  ’.LF.'J*  ».J 
NFZSNFINB(1. Z.l.f . ATRIBC1) if) 

IF (NFZ.EQ.I) RETURN 
NF«tFIND(l.  JH.  l.f  .ATRIB(l)  *•) 

IF(NF.E8.I)RETURN 

NFHFINDU.l.l.f.ATRIlUl.f) 

PRINT*. 'NF1*  '.NF1 
IF(NFl.EO.f) RETURN 

C  IF  THE  IN  KILLS  TIE  RADAR.  THEN  THE  RADAR  IS  CHECKED 

C  TO  SEE  IF  IT  HAS  A  EH/CCI  RADAR.  IF  IT  HAS.  THEN  ALL 

C  ASSOCIATED  RADARS  THAT  USE  THIS  EH/ttI  FOR  ACQUISITION 
C  AND  TRACKING  HAVE  THEIR  5TH  ATTRIBUTE  SET  TO  If. 

C  THE  "KAO*  RADAR  IS  REMOVED  FROM  FILE  l.FILE  2  AND  THE 

C  UEAS&’S  HORNING  FILE.  H(49)  IS  SET  EQUN.  TO  THE 
C  RADAR’S  SEQUENTIAL  NUHSER(ATRIKl)  FOR  THIS  EVENT). 

C  TIC  RADAR  IS  REMOVES  FROM  THE  EVENT  CALENDAR. 
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649 

651 


IF(PK(J).EQ.1)TH£N 

PRINT** *RAOAR  »»ATRIB(1>.'  KILLED  BY  W  ' >ATRIB ( 14) 

CALL  RNQVE(ffl.l.A) 

PRINT*! 'A(l) =  'iA(l) 

CALL  RHOVE(NF.JU.A) 

IF(A(2).EQ.6)THEN 
DO  66  ITslfWWU) 

CALL  COPT(IT.l.C) 

IF(A(5).EQ.C(5))THEN 
CALL  IWOVEUT.l.C) 

C(5)'ll.l 
CALL  FILEHU.C} 

ENDIF 

66  CONTINUE 
ENDIF 

CALL  FILEH(JU+t»A) 

CALL  RH0VE(NF2.2.A) 

H(49)*ATRIB<1) 

91  NREHsNFIND(liNCUjRi3!0iATR!B(l)»i) 

IF(NREH.CT.0)THEN 
CALL  ULINK(WEHiNCLNR) 

CO  TO  91 
ENDIF 

C  IF  THE  RADAR  IS  NOT  KILLED.  IT  IS  REMOVED  FROM  THE  HU'S 
C  WORKING  FILE  AND  ITS  14TH  ATTRIBUTE  IS  RESET  TO  0  IN 

C  FILE  2. 

ELSEIF (PK ( J) .EQ.0) THEN 
CALL  RHOVE(NF.JU.A) 

CALL  RN0VE(NF2.2>A) 

A(14)*0.0 
CALL  F!LEH(2iA) 

ENDIF 

RETURN 

C 

C  EVENT  104UU  NODE  ONLY)  -  REPOSITIONING  ROUTINE. 

C  EVENT  REALIZED  WHEN  W  TURNS  179  DECREES  ABSOLUTE 

C  RELATIVE  BEARING  FROM  SITE  IT  IS  ATTACKING. 

C  EVENT  ROLLS  W  OUT  OF  TURN. 

C 

10  PRINT*. 'EV  10  TNOU*  '.TNQM.»  ATRIB(l)*  ’.ATRIB(l) 

RATED)  *0.0 

RETURN 

imH«mH***M***M*m**tHHmHIHimHHW»«*IH»MH»m 

c 

C  EVENT  11 (UN  NOSE  ONLY)  -  REPOSITIONING  ROUTINE. 

C  EVENT  REALIZE1  WEN  HU  DETECTS  REQUIRED  RANGE  FROM 
C  SITE  IT  IS  ATTACKING.  EVERT  ROUS  UU  BACK  INTO  TURN. 

C 

11  PRINT*. ’EV  11  TNOU*  '.TNOU.'  ATRIB(l)*  '.ATRII(l) 
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700 
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733 
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735 
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738 
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742 

743 

744 

745 

746 

747 

748 


RATE! J) =-4.0 

IF (REV ( J) .EQ. 1 1 RATE ( J) =4.0 
RETURN 


c 

C  EVENT  12 (UN  NODE  ONLY)  •  UU  SENT  HONE. 

C  THIS  EVENT  SENDS  THE  UU  HONE  AFTER  IT  FIRES  ALL 
C  OF  ITS  ARHS. 

C 

12  PRINT*! *£V  12  TNOH=  »»TNOH»*  ATR1B(1)=  »»ATRIBC1) 

RATE (J) =4.1 

REV(J)-1.0 

RETURN 

m*************************************** ************************** 

C 

C  EVENT  13 (UU  NODE  ONLY)  -  UU  SENT  HONE. 

C  THIS  EVENT  ROLLS  UU  OUT  AFTER  ITS  INITIAL  TURN  HONE. 

C 

13  PRINT*! ’EV  13  TNQU*  ’,TNQU.'  ATRIB(l) =  'iATRIB(l) 

RATE(J)=0.0 

RETURN 


C 

C  EVENT  144ALL  AIRCRAFT)  -  THIS  EVENT  APPLIES  FOR  ALL 

C  AIRCRAFT.  IT  SIHULATES  RADARS  SEARCHING  FOR  AIRCRAFT. 

C  A  CHECK  IS  FIRST  HADE  TO  ENSURE  THAT  NO  MORE  THAN 
C  5  RADARS  ARE  ENGAGED  U1TH  THE  AIRCRAFT.  IF  THERE  ARE 

C  5  RADARS  ALREADY  UORKINC  THE  AIRCRAFT  THEN  THE  EVENT 

C  IS  BT-PASSED.  OTHERUISE  SUBROUTINE  SEARCH  IS  CALLED. 
C 

14  IF(SS(I).CT.145000.AND.ATRIB<1).LT.200)THEN 

XX (55) SATRIB(1) 

RETURN 
END  IF 

IF<NNQ(LF).E0.5)R£TURN 
CALL  SEARCH 
RETURN 


C 

C  EVENT  15 (ALL  AIRCRAFT)  -  EVENT  REALIZED  AT  THE  END  OF 

C  TRACKING  AND  ACQUISITION  TINE.  SUBROUTINE  PROB  IS  CALL- 

C  ED  TO  DETERMINE  PKi  TKTINE  OF  INTERCEPT)!  TLtTIME  OF 
C  LAUNCH) .  BASED  ON  THESE  VALUE  THE  FOLLOUINC  CAN  OCCUR. 

C  1)  EVENT  16  IS  SCHEDULED  AT  ESTIMATED  LAUNCH  TIME. 

C  2)  EVENT  17  IS  SCHEDULED  AT  ESTIMATED  IMPACT  TINE. 

C  3)  EVENT  19  IS  SCHEDULED  TO  OCCUR  IN  30  SECONDS  BAS- 

C  ED  ON  A  LOU  VALUE  OF  PK. 

C 

15  PRINTli'EV  15  THOU*  »iTN0H.'  ATRIB(l)*  '.ATRIBUI 
CALL  PROB 
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751 


IF<PKR.LT. .fS) THEN 
CALL  SCHDL(19i30.#iATRIB) 
RETURN 
ENDIF 

NPl=NFIND(lilil»iiATRIB(3) »i) 
PRINT*. »NP1=  '.NPl 
IF (NP1 .EQ.0) THEN 
CALL  SCHDL (19.30.1. ATRIB) 
RETURN 
ENDIF 

CALL  COPT(NPl.l.A) 

ATRIB(4)=TL 
ATRIB(5)=TI 
ATRIB (6) =PKR 
ATRIB (7) = A (6 ) 

ATRIB (8) =A (3) 

ATRIB (9) =A(4) 

IF (TNOU.EQ.TL) THEN 
ATRIB ( 10) =RI 
T=TI-TNOH 

CALL  SCHDL  (1 7.  T»  ATRIB) 

ELSE 

T=TL-TMO« 

CALL  SCHDLU6.T.  ATRIB) 

ENDIF 

RETURN 


C 

C  EVENT  16(ALL  AIRCRAFT)  -  EVENT  REALIZED  AT  ESTIHATED 
C  MISSILE  LAUNCH  TINE.  SUBROUTINE  PROB  IS  CALLED  TO  DE- 
C  TERNINE  NEU  PK  AND  TI.  BASED  ON  THESE  VALUES  EITHER 

C  EVENT  19  IS  SCHEDULED (PK  TOO  LOU)  OR  EVENT  17  IS 

C  SCHEDULED  AT  MISSILE  I HP ACT  TINE. 

C 

16  PRINT*, »EV  16  TNOU=  '.THOU.'  ATRIB(1)=  ’iATRIB(l) 

CALL  PROB 
IF(PKR.LT..I5)THEN 
CALL  SCHDL(19i3l.ltATRIB) 

RETURN 

ENDIF 

ATRIB(6):PKR 
ATRIBL5) *TI 
ATRIB ( If) *RI 
TTsTI-TNOU 

CALL  SCHDLU7.TT, ATRIB) 

RETURN 


I I 111 H************************** *************** *1***************** 


C 

C  EVENT  17 (ALL  AIRCRAFT)  -  EVENT  REALIZED  AT  MISSILE  IMPACT 

C  TINE.  SUBROUTINE  PROB  IS  CALLED  TO  DETERMINE  PK  AND  NEU 
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799 

801 

801 

80Z 

803 

804 

805 
804 

807 

808 

809 

810 
811 
812 
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815 
814 

817 

818 
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820 
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822 

823 

824 

825 
824 

827 

828 

829 

830 

831 

832 

833 

834 

835 
834 

837 

838 

839 

840 

841 

842 

843 

844 

845 
844 

847 

848 


C 

c 

c 

c 

c 

17 


55 


TI.  IF  THE  AIRCRAFT  HAS  TURNED  AUAT  FROM  THE  RADAR 
SITE ("RATIO"  GREATER  THAN  1.1)  THEN  ANOTHER  ITERATION 
IS  HADE  AND  EVENT  17  IS  SCHEDULED  AT  THE  TIHE  IT  TAKES 
FOR  THE  AIRCRAFT  TO  HOVE  TO  THE  NEW  TI.  OTHERWISE.  THE 
KILL/NOT  KILL  IS  EVALUATED  BASED  ON  HONTE  CARLO  DRAW. 

PRINT*. 'EV  17  TNOW*  ».TNOW»»  ATRIB(l)*  »»ATRIB(1) 

PRINT*. 'ATRIB(3)=  ».ATRIB(3) 

CALL  PROS 

T*  ( (ATRIB  (5)  -ATRIB(4)  )*SRR)  /ATRIBU0) 

RATIO=SRR/ATRIB ( 10) 

IF (RATIO.CT .1.1) THEN 
ATRIB(4)*PKR 
CALL  SCHDL(17.T.ATRIB) 

RETURN 

ENDIF 

SAHPLE*DRAND(l) 

rr/CAMptr  it  ATQIBfAHTUni 

PRINT*.' AIRCRAFT  '.ATRIBU).'  KILLED  BT  RADAR  >.ATRIB(3) 
IF(ATRIB(1).GT.200)THEN 
11(591*11(59) *1.0 
IF (ATRIB ( U . EQ.21 1 ) XX (54) -ATRIB ( 1 ) 

IF(ATRID( 1 ) .£0.212) XX (58) *ATRIB ( 1 ) 

ENDIF 

IF(ATRIB(1).LT.200)XX(57)*XX(57)+1.0 
NR£H*NFIND(l.NCLNR.1.0iATRIB(l) .0) 

IF(NREH.GT.0)THEN 
CALL  ULINK(NREN.NCLNR) 

CO  TO  55 
ENDIF 

ATRIB (4)*1.0 
ENDIF 

CALL  SCHDL(18i30.0.ATRIB) 

RETURN 


w 

C  EVENT  18 (ALL  AIRCRAFT)  •  EVENT  OCCURS  30  SECONDS  AFTER 
C  HISSILE  IHPACT.  BASED  ON  WHETHER  THE  AIRCRAFT  WAS  KILLED/ 
C  NOT  KILLED  RABAR(S)  ARE  FREED  TO  BEGIN  NEW  SEARCH. 

C 

18  PRINT*.'EV  18  TNOW*  '.TNON.'  ATRIBdl*  '.ATRIB(l) 

IF (ATRIB <4),NE.l) THEN 

C  IF  THE  AIRCRAFT  WAS  NOT  KILLED  THEN  ONLT  FREE  THAT  RADAR 
C  THAT  WAS  SHOOTING. 

NC*NFIND ( 1 . LFi 1 .0. ATRIB (3) .0) 

IF(NC.EQ.0)GO  TO  200 
CALL  RNOVE(NGfLF.A) 

NC1*NFIND  (Id.  1 .0.  ATRIB  (3)  .0) 

PRINT*. 'NCI*  '.NG1 
IFOIGl.EOJIGO  TO  210 
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851 
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896 
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898 


CALL  RHOVE(NCl.l.A) 

A(8)s0.0 
CALL  FILEH(l.A) 

2M  RETURN 
ELSE 

C  IF  THE  AIRCRAFT  WAS  KILLED  THEN  FREE  ALL  RADARS  TRACK- 
C  INC  IT. 

112  NR*NFIND(1»LF.8.0.ATRI8U).0> 

IF (NR.CT.i) THEN 
CALL  RHOVE(NR.LFiA) 

CD  TO  102 
ENDIF 

103  NR1*NFIND(1.1.8.0.ATRI8(1).0) 

IF (RR1 .CT.0) THEN 
CALL  RHOVE(NRl.l.A) 

A(8) =0.0 
CALL  FILEH(l.A) 

CO  TO  103 
ENDIF 
ENDIF 
RETURN 

c 

C  EVENT  191ALL  AIRCRAFT)  -  EVENT  IS  SCHEDULED  FROM  EVENTS 
C  15  OR  16  UHEN  THE  CALCULATED  PK  IS  DETERMINED  TOO  LOU 
C  FOR  THE  RADAR  TO  CONTINUE  TRACKING. 

C 

19  PRINT*.'EV  19  TN0H='.TM0H»’AT(l)s'.ATRIB(l).'AT{3)s,*ATRIB(3J 
NPSNFIND(1.LF*1*0.ATRIB(3)»0) 

PRINT*. »NP=  ».NP 
CALL  RHOVE(NP.LF.A) 

NP-NFIND(l*lil*0*ATRlB<3)*0) 

PRINT*. 'NP*  ».NP 
CALL  RHGVE(NP.l.A) 

A  (81-0. 

CALL  FILEH(l.A) 

RETURN 

END 


C 

C  SUBROUTINE  SEARCH  IS  CALLED  BT  EVENT  14.  THIS  SUBROUTINE 
C  SIMULATES  THE  SEARCHING  BT  RADARS  FOR  AIRCRAFT.  BASED  ON 
C  VARIOUS  PARAHTERS  THE  RADAR  HILL  DETECT  THE  AIRCRAFT  AND 
C  SCHEDULE  TRACKING  AND  ACQUISITION  TIIC. 

C 

SUBROUTINE  SEARCH 

COIVIM/SCOHl/ATRIB(10f).DD(lfl).DDL(100).BTNOH.II>HFA.HSTQP.NCLNR. 
tNCRDR.NPRNT.NNRUN.«SET.NTAPE.SSaN)  .SSL(lff)  .TNEXT. 

(THOU. II (101) 
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COmON/UCOKl/RATE  (2)  .L(2) 

C0RH0N/UC0N3/SAB(5.12) 

DIHENSION  A ( 16) .B(16) iBA( 16) 

I:ATRIB(1) 

IFd.LT.2M)  THEN 
ALT*XX(13) 

TP0S=XX(6i+I> 

LF*I+6 

ELSE 

1=1-191 

J=I-2I 

YP0S=SS(I+4J 

ALT=XX(3) 

LF=I-4 

ERDIF 

SRHP= ALT /S 1ND ( 0 . 25 ) 

LL=WFEd) 

CALL  COPY (-LL> liBAl 

IF  THE  RADAR  IS  AN  EU/CCI 
IF (BA (2) .EQ.6) CO  TO  27 

OR  IF  THE  AIRCRAFT'S  ALTITUDE  IS  BELOH  RADAR  COVERAGE 
IF(ALT.LT.BA(7))C0  TO  27 

OR  IF  THE  RADAR  IS  ALREADY  ENGAGED 
IF(BA(8).CT,i)G0  TO  27 

OR  IF  THE  RADAR  IS  SHUTDOHN  »»  CO  TO  THE  NEXT  RADAR 
IF(BA(9).E0.f)C0  TO  27 
X=BA(3) 

Y*BA(4) 

SR=SQRT(  (X-SS(I)  J«2+(Y-YP0S)h2+ALT«2) 

IF  THE  RANGE  TO  THE  AIRCRAFT  IS  LESS  THAN  NULTI-PATH 
IF(SR.CT.SRNP)GO  TO  27 

OTHERHISEt  IF  THE  RANGE  IS  LESS  THAN  THE  NAXIHUN  RANGE 
IF(SR.LT.BA(61)THEN 
BA(8)=ATRIBd) 

CALL  FILER (LF.BA) 

CALL  RROVE(-LL.l.A) 

CALL  FILERd.BA) 

PRINT*»’Ad)=  ’.Ad) i'  BAd)=  ’.BA(1) 

IT*BA(2) 

DETERMINE  ACQUISITION  AND  TRACKING  TIRE  BASED  ON  UNIFOI 
DISTRIBUTION  UNLESS  THE  RADAR'S  ASSOCIATED  EU/CCI  IMS 
KILLED(BA<5)=1I)  •  IN  THAT  CASE  USE  ONLY  THE  HIGH  TIRE 
TL*SAN(IT»7) 

TH*SAN(IT.8) 

TRMJNFRROL.TH.l) 

IF (BA (5) .EQ.lf)TRC*TH 

Bd)*ATRIBdl 

8(3)*BA(1) 

CALL  SCHSLdSiTRC.B) 
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IF(NM(LF) . E0.5) GO  TO  26 
27  LL=NSUCR(LL) 

IF (LL.NE.0) CO  TO  25 
26  CONTINUE 
RETURN 
END 


C 

C  SUBROUTINE  PROB  -  THIS  SUBROUTINE  CALCULATES  THE  PK'S 

C  FOR  THE  APPROPRIATE  STSTEHi  AS  NELL  AS  THE  TINE  OF  LAUNCH 

C  (TL)  AND  TIME  OF  IMPACT (TI). 

C 

SUBROUTINE  PROB 

COHHON/SCOHl/ATRIBUH)  >DD(10f>  iDDL(IM)  .QTNOU,II>MFA.MSTOP»NCLNR» 
»CRDR>NPRNTiNNRUN>NNSET>NTAPEiSS(lN)  >SSL(10f )  *THEXT» 
tTNOUflX(lM) 

C0HH0N/UCQH2/TL  t  T I » PKR » RI » SRR 
COHHQN/UCOMI/RATEIZ) *L(2) 

DIMENSION  A(16)«B(16hBA(16) 

COMHON/UCOH3/SAH(5> 12) 

DATA((SAH(IA»JA)fJAslil2) fIA=l»5)/ 

Uif(l>lififi6i25ifi-47.6ilt3fi 

l2.52E-3»961l»671i7.3£-27i6.36E-17i671»18i5lf592»-46.Zi56.4» 

113* 

t5.62E-6»25ii*232*7.3E-27t6.36E-L7*232> 12*26. 759* -51.4*43.6* 

(ISi 

«7.lE-7,22Mf58»4.l3E-27fl.i2E'14*58fl7»38.599,-53.4*26.2, 

(3f* 

13.25E-7*189«*2S»8.09E-27i4>84C-16t2S* If * 23*525»-S2.4* 22* 30/ 
I=ATRIId) 

IFd.LT.2M)  THEN 
YP0S=H(6f+I) 

ALT*II(13) 

VEL*n(14) 

HDC=9f.f 

ELSE 

1=1-19# 

J*I-19f 

HDG*SSd+8) 

TP0S*S$U+4) 

ALT=IX(3) 

VEL=H(5) 

FNDIF 

SRHP=ALT/SIND(f.Z5) 

IF(I.LE.l#)LF=I+6 

IFd.CT.ll)LF*I-4 

C  GET  THE  FILE  OF  THE  AIRCRAFT  BEING  ATTACKED. 
NG=NFIND(l*LFil«l*ATRIB(3) *•> 

CALL  COPT(NC»LF*A) 
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m 

X=fl(3) 

lfif 

r=fl<4> 

INI 

SR^SORTI  (X-SS(I)  )h2+(MP0S)«2+ALTh2) 

1N2 

AKCNS57 . 3*68 IN (ALT/SRI 

1113 

SRfcSR*COSD(ANCN) 

1M4 

C 

IF  THE  RANGE  IS  GREATER  THAN  HULTI-PATH  OR  GREATER  ' 

ins 

C 

MAXIMUM  RANGE  OF  THE  THREAT  SET  PK  EQUAL  TO  0. 

1006 

IF(SR.GT.SRMP.0R.SR.CT.A(6) ) THEN 

1M7 

PKR*0.0 

1M8 

RETURN 

1M9 

END  IF 

101# 

SRR-SR 

1111 

c 

IF  THE  THREAT  IS  A  AAA  CO  TO  AAA  ROUTINES) 

1112 

IF(A(2).E8.1)C0  TO  99 

1113 

c 

DETERMINE  ASPECT  ANGLE  OF  THREAT/AIRCRAFT. 

1014 

D=  tT-TPOS) / (X-SS(I) 1 

1015 

YT=T-YPOS 

1016 

IF  CD-GE.0. AND . TT .GE.i) THEN 

1017 

FANG:57.3<ATAN(1/D) 

1018 

ELSEIF (D.LT.0. AND. TT .GE.I) THEN 

1019 

FANC=270+57.3*ABS(ATAN(D)) 

1020 

ELSEIF (D.LT.I. AND.TT.LT.il  THEN 

1021 

FA0C=90+57 .3*AB$ ( ATAN (D) ) 

1022 

ELSEIF (D. GE.I. AND. TT.LT.0) THEN 

1023 

FANC=270-57 . 3*ATAN (Dl 

1024 

END  IF 

1025 

ANC=ABS(FANC-HDC) 

1026 

IF(ANC.C7. 180) ANC-360-ANC 

1027 

c 

IF  THE  ANCLE (ANC)  IS  LESS  THAN  90  DECREES  THEN  THE 

1028 

c 

AIRCRAFT  IS  IN  FRONT  OF  THE  THREAT.  CALCULATE  PARA¬ 

1029 

c 

METERS  TO  OBTAIN  TI  AND  TL. 

1030 

IFIANC.LT. 90) THEN 

1031 

RC-SR«SINO(ANC) 

1032 

AB*SR*COSD(ANC) 

1033 

TIKAC-AB/VCL 

1034 

TMF0sRC/SAM(A(2)i9) 

1035 

IF(TIHAC.CE.THFO)THEN 

1036 

TL=THOH+TIMAC-TMFO 

1037 

TI*THFO+TL 

1038 

RI*RC 

1039 

ELSE 

1040 

RI*SORT(  ( (THFO-  TINACJ  «VEL)  «2+RC»2) 

1041 

TL*TNQU 

1042 

TI*TL* (RI/SAH (A (2) »9) ) 

1043 

ENDIF 

1044 

ELSE 

1045 

ANCZM80-ANC 

1046 

RC*SQRT  ( (SRNtSIND  (ANC) )  «2+ALT«2) 

1047 

VRsVEL/ (SAH(A(2) >9) ) 

1048 

E»t.0-(VR«2) 

162 


184? 

1151 

1151 

lf52 

1153 

1154 
lf55 
1856 

1157 

1158 

1159 
1861 
1861 
1862 

1863 

1864 

1865 

1866 

1867 

1868 
1869 
1878 

1871 

1872 

1873 

1874 

1875 

1876 

1877 

1878 

1879 
1888 
1881 
1882 

1883 

1884 

1885 

1886 

1887 

1888 
1889 
1898 

1891 

1892 

1893 

1894 

1895 

1896 

1897 

1898 


F»- (2*$RN*VR*C0SD (AMC) ) 

G-SQRT({FM2)+(4*E*(SRN«2)1) 

RI*(C-F) / (2*E1 
RI*S0RT(R1«2+ALTH2) 

TL*TNQW 

TI=TL+(Rl/SAH{fl(2) »9) ) 

END  IF 

TRCS=188>57.3*(ASIM(RC/Rin 

IJ*A(2) 

1CET-NFIN0(1> 19» 1 f if TRCSi 2.51) 

CALL  C0PY(ICETfl9»B) 

JR*IJ+1 

SI6NA*B(JR) 

IF(RI.CE.A(6))TH£N 

PKR*8.B 

ELSE 

!F(ATR1B(1).LT.288)TH£N 
AJSBDUsSAH (A (2) i 18) +28*AL0Cli (RI ) -SIGHA 
AJS*1B«(AJSBDH/18) 

CEP-SORT  ( ( ( (SAN (A<2) » 1 )  *AJS*  <RI«2) ) +SAH  (A(2)  >2)  *AJS)  ♦ 
tSAN(A(2)i3))l 

ass 

SICN2MB**(SIGNA/18) 

CEP=SORT  ( ( (SAN  ( A  (2)  1 4)  *  (RI«6)  /SICN2) + (SAN  (A  (2) » 5)  * 
t  (RI«4)  /SIGN2) +SAN  ( A  (2)  r  6) ) ) 

ENDIF 

PKR-1.8-(.5«(((SAH(A(2)»11)/CEP)«2))) 

IF(PKR.LT..B2)PKR=8.B 

ENDIF 

RETURN 

99  IF(ATRIB(11.CT.2BB)THEN 
IF  (RATE(J)  .NE. BITTEN 
G-2.8 

ass 

6*1.3 

ENDIF 

ELSE 

IF(SS(I1  .GT.45888.8.AND.SS(I)  .LT.658M.8)  THEN 
6*2.8 
ass 

6*1.3 

ENDIF 

ENDIF 

SRX*SR/1888.8 

SI6AAA2*(28»$RK)«2.8 

VF*938.8»EIP(-.4965*SRK) 

TOFL* (2814.46/VF) *2. 166 
TI«TN0B*T0Fl 

DA*(2.»3.14159*SI6AAA2)+5.1727 

PKSS*(5.1727/DA)*£XP((-,5»<9.8*C*(T0FL«2.))m2.B)/DA) 
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1199 

11M 

mi 

1112 

1113 

1114 
HIS 
lift 
1117 
UK 
1119 
111! 
1111 
1112 

1113 

1114 
HIS 
Hit 


PKR*l.f-((l,l-PK$S)*«5l.) 

TL-THOM 

RETURN 

END 

MHUHHHUmmmtUHWHHIHHHIHHWmmiHHHW 
mHHinim»inimmniinninmmnimmiinnimnin 
SUBROUTINE  OTPUT 

C0*0N/SC0H1/ATRIB<1H)  .MHIH)  .OBL(IH)  .DTNQNiIIiRFAiRSWNCLNRi 
UICRI)R>NPf9ITfNNRUNiNNSETiNTAPEiSS(lfl)  iSSL(IN)  iTNEITi 
tTNOUtH(lff) 

PRINT* 

PRINT* 

PRINT** 'ATTACK  AIRCRAFT  SURVIVING*  MX(5t) 

PRINT*. 'ATTACK  AIRCRAFT  KILLED  *  MX (57) 

PRINT*. ’RADARS  KILLED  BT  UK  *  MX  (51) 

PRINT*! 'UK  KILLED  *  MX  (59) 

RETURN 

END 
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Appendix  H 


The  SPSS  program  for  the  two-way  ANOVA  was  run  with 
the  data  contained  on  the  following  page.  The  first  column 
is  the  number  of  aircraft  surviving  and  reaching  the  target 
area  for  the  run,  the  next  column  the  level  for  the  factor 
of  altitude  (l=60m,  2*200m) ,  and  the  last  column  the  level 
for  the  factor  of  tactic  (l=ahead  of  the  attack  force,  2- 
with  the  attack  force) .  The  overall  results  for  the  ANOVA 
are  listed  on  the  page  following  the  data.  There  were  no 
significant  main  effects  or  interactions. 


1  VOCELBACK  COMPUTINC  CENTER 

2  NORTHUESTERN  UNIVERSITY 

3 

4  SPSS--  STATISTICAL  PACKAGE  FOR  THE  SOCIAL  SCIENCES 

5 


6 

UMtm  ANALYSIS  OF  VARIANCE* 

*  *  *  * 

<  *  * 

7 

SURVIVE 

8 

BY  ALT 

1:68Ni2s288H 

9 

TACTIC 

i=AHEAD»2=UITH 

11 

HHHHH 

*  t  # 

*  *  *  *  * 

*  *  «  * 

*  *  * 

11 

12 

13 

SUN  OF 

MEAN 

SICNIF 

14 

SOURCE  OF  VARIATION 

SQUARES 

OF 

SQUARE 

F 

OF  F 

15 

16 

MAIN  EFFECTS 

4.988 

2 

2.458 

1.114 

.353 

17 

ALT 

2.458 

1 

2.458 

1.114 

.387 

18 

TACTIC 

2.458 

l 

2.458 

1.114 

.317 

19 

28 

2-UAY  INTERACTIONS 

.458 

1 

.458 

.285 

.657 

21 

ALT  TACTIC 

.458 

1 

.458 

.285 

.657 

22 

23 

EXPLAINED 

5.358 

3 

1.783 

.811 

.586 

24 

25 

RESIDUAL 

35.288 

16 

2.288 

26 

27 

TOTAL 

48.558 

19 

2.134 

28 

29 

31 


2f  CASES  HERE  PROCESSES 


WW  Attack  Profile  Validation 


This  appendix  contains  the  validation  of  the  WW 
attack  profile.  As  the  WW  proceeds  through  the  threat 
environment  hunting  for  radars  to  attack,  the  logic  associ¬ 
ated  with  its  profile  and  the  necessary  computations  that 
change  the  WW's  behavior  in  the  attack  are  hand-calculated 
and  compared  to  computer  output  data.  A  typical  WW  attack 
is  followed  as  the  WW  moves  from  one  logic  event  to  the 
next. 

To  hand-calculate  the  required  parameters  an  encounter 
geometry  is  defined  in  the  encounter  diagram.  Position  1 
(xl»  y^  depicts  the  WW  in  the  FEBA  geometry  as  it  begins 
a  turn  and  Position  2  (x2,  y2)  when  it  rolls  out.  The  WW's 
flight  distance  from  Position  1  to  2  is  given  by  the  follow¬ 
ing  formula: 

S  -  R9 

where  s  *  distance  WW  flies  from  1  to  2, 

R  »  WW's  radius  of  turn,  and 

8  *  number  of  degrees  WW  turns  from  1  to  2 
(in  radians) . 
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I 1 

c  P 

ill 


I 

it  §* 

a  a  *  j 
<£««* 

cLi 


The  radius  of  turn  R  for  an  aircraft  is  given  by: 


where  R  *  radius  of  turn, 

v  *  aircraft  velocity,  and 
TR  ■  aircraft  turn  rate  in  radians  per  second. 

Thus  for  the  WW  with  a  velocity  of  247  m/sec  and  a  turn 
rate  of  4  degrees /sec,  R  is  given  by: 


» 


t 


247 


4  x 


3538  m 


For  a  turn  rate  of  2  degrees  per  second,  R  =  7076  m. 

By  knowing  the  WW's  position  when  it  starts  its 
turn,  the  number  of  degrees  of  heading  change,  and  the 
radius  of  turn  R,  then  Position  2,  the  rollout  point,  can 
be  determined. 

WW1,  call  sign  211,  is  followed  on  its  attack  pro¬ 
file  as  it  hunts  for  threats  to  attack,  detects  a  threat, 
performs  a  ranging  routine,  launches  an  ARM,  and  begins  a 
new  engagement  sequence.  Output  data  from  the  computer 
simulation  is  listed  at  the  end  of  this  appendix. 

WW1  enters  the  threat  scenario  at  time  TNOW  «  0, 
x  and  y  coordinates  at  (0,  25110),  heading  090,  altitude 
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60  m.  Its  radar  horizon  or  line-of- sight  range  is  given 
by  equation  (2) : 

h  -  4117.3  ^60  «  31892 

The  closest  radar  that  WW  can  attack  is  radar  64  located 
at  (52500,  25139).  (See  Appendix  J  for  radar  coordinates.) 
WW1  will  proceed  across  the  FEBA  until  the  distance  to  radar 
64  is  detected  to  be  less  than  the  radar  horizon  distance 
at  which  time  event  1  will  be  called.  From  the  computer 
output  event  1  is  realized  as  TNOW  =  84.  WWl's  position  at 
that  time 


x  *  247  x  84  -  20748 
y  -  25110 

The  distance  to  radar  64  is  given  by  equation  (14) : 

SR  *  ^ (20748  -  52500)  2  +  (25110  -  25139)2+  602 

-  31752 

Event  1  occurs  at  the  correct  time  and  position. 

Event  1  starts  the  WW  ranging  routine  turning  the 
WW  away  from  radar  64  at  a  turn  rate  of  4  degrees /second. 
At  event  1  WWl's  relative  bearing  to  radar  64  is  given  by 
equation  (11) : 


.__-l  F  52500-207481 
R  ■  tan  ffgvjnr 'vw  i  * 


-  89.9° 


Event  2  occurs  at  TNOW  *  102  and  the  absolute  rela¬ 


tive  bearing  to  the  site  should  be  75  degrees.  The  heading 
change  from  event  1  to  event  2  is  given  by  equation  (10): 

0  «  (102-84)  x  4°/sec  -  72° 

Its  new  heading  at  TNOW  -  102  should  be: 

H  *  90°  +  72°  *  162° 

From  the  encounter  geometry,  the  position  of  WW1  at  event  2 
should  be: 

x  =  20748  +  3538  cos  18°  =  24113 
y  *  25110  -  [3538  +  3538  sin  18°]  *  22665 

Relative  bearing  to  radar  64  at  this  point  is,  by  equa¬ 
tion  (11) : 


p  ,  *  -1  f 52500-241131  = 

R  tan  L2$140-22665j  85 

The  absolute  relative  bearing  is  given  by  equation  (12) : 

AB  -  (162  -  85)  «  77° 

These  figures  all  agree  with  the  output  data. 

Event  3  occurs  when  the  absolute  relative  bearing 
is  105  degrees.  The  output  data  lists  TNOW  ■  159  for 
event  3.  The  position  of  WW1  at  this  time  is  given  by 
equations  (35)  and  (46): 
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'V 


x  *  24110  =  247  sin  162  x  (159-102)  -  28461 


y  =  22670  +  247  cos  162  x  (159-102)  -  9280 
The  relative  bearing  to  radar  64,  by  equation  (11)  is: 


R 


tan 


-1 


r52500-2846ll 
L25140-9280  J 


56.6° 


The  absolute  relative  bearing,  given  by  equation  (12)  ,  is: 
AB  =  (162  -  56.6)  -  105.4° 


Hand-calculations  confirm  output  data. 

Event  4  occurs  when  the  absolute  relative  bearing 
is  10  degrees.  Event  3  should  turn  the  WW  back  towards  the 
site  at  4  degrees /second.  Event  4  occurs  at  TNOW  *  186. 

The  heading  change  is  given  by  equation  (10) : 

6  =  (186  -  159)  x  4° /sec  «  108° 

WWl's  heading  at  event  4  should  be: 

H  »  162  -  108  »  54° 

From  the  encounter  geometry  WWl's  position  is  calculated: 

x  -  28460 + 3538 (cos  18  +  cos  54)  -  33908 
y  -  9279  +  3538  (sin  54  -  sin  18)  =  7510 

Relative  bearing  to  radar  64  at  this  point,  by  equation 
(11): 
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R 


1  152500-33908 
125140-7510  . 


6.5 


Absolute  relative  bearing  is  given  by  equation  (12) : 

AB  *  (54  -  46.5)  =  7.5° 

All  calculations  agree  with  the  output  data. 

Event  5  occurs  when  WW1  is  bore sighted  on  radar  64 
(absolute  relative  bearing  zero) .  Event  4  turns  the  WW 
into  the  threat  at  a  rate  of  2  degrees /second.  When 
event  5  occurs  at  TNOW  «  190,  heading  change  is  given  by 
equation  (10): 


9  -  (190  -  180)  x  2°/sec  *  8° 

WWl's  heading  should  be: 

H  -  54  -  8  -  46° 

From  the  encounter  geometry  the  WW  position  is: 

x  -  33910  +  7076  (sin  44  -  sin  36)  «  34618 
y  -  7510  +  7076  (cos  36  -  sin  46)  -  8146 

Relative  bearing  from  equation  (11)  is: 


Q  _  «.  -1  T 52500-346181 

R  ■  tan  I  ... ac  I 


46.4 


The  absolute  relative  bearing,  equation  (12) ,  is: 


AB  «  (46  -  46.5)  =  0.5° 

These  calculations  agree  with  the  output  data. 

Event  5  calculates  the  ARM  release  time  and  posi¬ 
tion,  the  time  of  flight  for  the  ARM  (TOF) ,  and  draws  a 
random  sample  to  evaluate  the  probability  of  kill  of  the 
ARM.  WWl's  distance  to  radar  64  at  event  5  is  given  by 
equation  (14) : 

SR  *  ^ (34660  -  52500) 2  +  (8144  -  25140)  2  +  622 
=  24639 

The  ARM  release  point  for  the  WW  is  taken  from  a  uniform 
distribution  of  low  value  equal  to  minimum  range  for  the 
ARM  and  high  value  equal  to  the  maximum  ARM  range  or  dis¬ 
tance  to  the  threat,  whichever  smaller  (as  long  as  ARM  is 
within  range  of  the  threat) .  Computer  output  sets  the 
release  point,  XX (4),  equal  to  24597.  The  time  for  WWl  to 
reach  this  point  from  its  present  position  is: 

Time  *  Distance  to  Release /WW  Velocity 

_  (24639  -  24597)  _  „  ,,  ___ 

*  - x  i  -j  1  ■'  m  u .  i  /  sec 


The  TOF  for  the  ARM  is: 

TOF  ■  ARM  Distance/ARM  Velocity 
-  24597/450  -  54.67 
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WWl  will  kill  radar  64  if  the  random  sample  drawn  is  less 
than  or  equal  to  0.85,  the  ARM  PK.  Event  5  indicates  that 
the  random  sample  drawn  is  0.6045.  Calculations  in  event  5 
concur  with  computer  data. 

Event  6  occurs  when  WWl  launches  the  ARM  at 
radar  64.  Hand  calculations  predict  this  should  happen  at 
event  5  time  plus  the  time  for  WWl  to  get  from  event  5's 
position  to  the  ARM  release  point: 

TNOW (Event  6)  *  190  +  0.17  =  190.17 

This  agrees  with  the  output  data. 

Event  9  simulates  ARM  impact  and  should  occur  when 
the  ARM  launched  by  WWl  reaches  radar  64.  This  time 
should  be: 

TNOW (Event  9)  -  TNOW (Event  6)  +  TOF 

-  190.15  +  54.67  -  244.82 

Since  the  sample  drawn  in  event  5  was  less  than  0.85,  the 
threat  should  be  killed,  as  in  fact  it  is  according  to  the 
output  data. 

After  WWl  fires  the  ARM  it  should  be  released  5 
seconds  later  to  start  another  engagement.  Output  data 
concurs  with  this. 

As  evidenced  by  the  above  comparisons  between  com¬ 
puter  output  data  and  hand-calculations,  the  WW  attack 
profile  is  validated. 
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1  ** INTERNED I ATE  RESULTS** 

2  EV  1  THOU=84  ATRIB(l) =21i  RADAR  64  IS  ENGAGED 

3  EV  2  TN0U-1A2  ATRIB(l)=2tl 

4  EV  3  TM0U=159  ATRIB(1)=2U 

5  EV  4  TNOU-186  ATRIB<1)=211 

6  EV  5  TNOtMVI  ATRIB)  11*211 

7  11(111=24435  XX(4) =24597 

8  T0F*54.&7  SANPLE=.6845 

4  pv-i  pt  uu-  tS4S 

If  EV  4  TNOU=198.154  ATRIB(1)=2U 

11  EV  7  TN0M=195.154  ATRIB(1)=2U 

12  EV  8  TN0U=195.254  ATRIB(1)=2U 


13 

EV  2  TNOU=2f7  ATRIB(1)=211 

14 

EV  9  TNOtl-244.8  ATRIB(1)=44 

15 

RADAR  64  KILLED  BT  Mi  211 

14 

■ 

17 

HTABLE  NUMBER  1** 

18 

TINE  XPOS  YPOS  HDNC 

RNGE  RLBR 

ABRLB 

19 

f  1  25111  91 

8  8 

8 

2f 

84  21751  25111  9f 

31758  89.96 

.839 

21 

112  24111  22671  162 

28498  85 

77.8 

22 

159  2844#  9279  162 

28888  54.6 

185.4 

23 

186  33911  7511  54 

25628  44.5 

7.5 

24 

198  34448  8144  44 

24448  44.4 

.4 

25 

195  35558  9882  46 

23488  44.4 

.4 

24 

287  34488  11458  359 

18928  75.5 

76.5 
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SAM  and  AAA  Probability  of  Kill  Validation 

Event  15  (at  the  end  of  acquisition  and  tracking) , 

event  16  (scheduled  missile  launch  time) ,  event  17  (scheduled 

>  missile  impact  time) ,  and  event  18  (freeing  the  threat  radar 

30  seconds  after  impact)  represent  the  discrete  events  asso- 

K  ciated  with  the  defensive  network.  Validating  this  por- 

« 

tion  of  the  model  involves  computing  the  variables  required 
to  estimate  the  P^  and  comparing  them  with  a  computer  print¬ 
out  of  model  derived  parameters.  A  set  of  calculations 
for  both  SAM  and  AAA  systems  will  be  developed. 

The  calculations  list  the  variable,  followed  by  the 
computer  variable  for  the  quantity  when  the  two  are  differ¬ 
ent.  For  example, 

SR  cos  9 
t^.  (TIMAC)  -  —2* - 

*  m 

where  "tAC"  is  the  notation  for  the  time  required  for  the 
♦  aircraft  to  reach  point  "C"  in  the  thesis  and  "TIMAC"  is 

the  computer  variable  of  the  same  quantity. 

SAM  Calculations 

* 

Strike  aircraft  "1"  was  in  threat  "54"'s  field  of 
view  (FOV)  at  TNOW  -  32.  Aircraft  "1”  was  created  at 
TROW  »  30  and  Subroutine  SEARCH  assigned  site  "54"  to  it  at 
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SAM  Probability  of  Kill  Computer  Printout 


> 


(a)  Subroutine  SEARCH  detects  an  aircract 

RADAR  START  T  K  AC  <  I  f^G  A~r*  1* 

AT  TIME  12. 

TRC*  23.26fU661>g«:g  _  _ 

(b)  Event  15:  At  the  end  of  Acquisition  and  Tracking 


%% 
a  \ 

’v  <3l 
%% 


EV  13  TNDW=  «5,26^»ei8?8t5  ATFIBd)*  1. 


RADAR  NUMBER  54. 

YPCS*  2- =91. 6.6242727  ALT*  11*.  VEL  =  2^7. 

HOG*  9°.  SFMP*  T li  <•€  ,  392  3  *1  L G*  7 

Xs  S  *  ' '  .  Y=  26  561.6?  3  73  3  3  SP=  i'3796.‘i3iZ33S  D  =  . 
ANG*  2,  32d-r2'.c2  *3=  1-6-,  6-  ‘233*6  TIM\C*  2idfcl 
TMFO*  2.:ib64  7193'£2  TLO  2 -*  .  219  ‘‘♦Si  36  Tit  272.^15- 
RI=  19  36  ''  21345  T?.3S=  8'. 99337  h7*8s  SIGMA*  25. 

PKP*  .9;  2  8  2.' 3i- 12  61  C£P=  21 . A~8  36 1 1 96& 

EV  15  TfOWs  ;=,  537561  St-5  iS  A'-'TSti)*  1. 


:  266 48 3 433 &7S 3 
4  2:  S7  7  5 
3  2:  33 


(c)  Event  16:  Missile  Launch  Time 


EV  13  7M3  W=  27  .29964  633  3  ATRI9<1>*  1. 

YP05*  24=91.56242727  ALT*  31  .  VEL*  2t7. 

HOG*  9  .  SFMP*  FI-.  46.992  I?*1?.  LG*  7 

X*  5  1  •  5  -  Y*  25361.83  373  .  9  SR*  20  96.4*1135199  3*  3.^  49  345  523  47 
ANG*  71. 34565453916  3C*  19  94.447188044  TIMAC*  2.6h‘-5%11;  .’753 


TMFO*  2.  6 2773r  155E 26  FLO  27*  ,  3165  '7  2 832  TI*  272.94*  2374337 
RI=  199'  .447188144  JRCS*  9  9,991376  47465  SIGMA*  25. 

PKP*  .9^2326  4  2332  3EP*  2  .  ♦9077932368 


(d)  Event  17:  Missile  Impact  Time 


» 


% 


•EV  17  TN  CM*  272.94^2374387  ATRIB(l)*  1. 

VPOS*  24591,66242727  ALT*  31  .  VEL*  247. 

HOG*  9'.  SRMP*  71:46.  392  3rG9  LG*  7 

X*  530^t.  Y*  26561. 8S  973  ';9  SR*  1994  .4  78355044  D*  -272.53351929"; 

ANG*  9C.2C354178149  *5  =  1994.465769825  TUAC*  . J5125742321L7 

THFO*  .8,98*32333339  TLO  2^2.9*42374  38  7  TI*  275.  7164  1  "64 

RI*  2i-’4 .772147843  T? CS*  1*8.5697372733  SIGMA*  3.7 

PKR*  • 1114t 71311733  :EP*  l "S. 6195 3 ’ 8 36 8 

AIRCRAFT  1.  KILLED  BY  RADA  8  54. 


(e)  Event  18:  Freeing  the  radar  30  seconds  after  impact 

EV  18  TNDW*_3C2.94/.237-,39  7  a^RIBCD*  i, 
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this  time.  Threat  "54"  is  a  type  SAM  B.  Specific  data  on 
both  the  threat  and  the  aircraft  are  specified  below: 


Threat  54  SAM  B 
x  coord  60,000  m 
y  coord  26,561  m 
Max  FOV  74,150  m 
Min  detec  mg  305  m 


Aircraft  1  Strike 
initial  x  coord  0 
y  coord  24,592  m 
x  velocity  247  mps 
Altitude  310  m 
Heading  090° 


Aircraft  1 ' s  y  coordinate  does  not  change  during 
the  mission.  Its  x  coordinate  at  any  time  can  be  determined 
as  follows: 

x  *  247 (TNOW-TCREATE)  (59) 


where 


x 
247 
TNOW 
TC REATE 


the  aircraft's  current  x  position, 

the  aircraft's  x  velocity, 

the  current  simulation  time,  and 

the  time  the  entity  was  created  by  the 
model. 


Aircraft  l's  x  position  at  TNOW  "  32  seconds: 
x  »  247(32-30)  »  494  m 

The  Slant  Range,  SR,  to  the  site  can  be  computed  by  solving 
equation  (14) : 

SR  -  y (26, 56 1-24, 592) 2 +  (60 ,000-4942)  + 3102 
-  59,540  m 
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Thus,  Aircraft  1  is  inside  site  54' s  maximum  detection 
range  of  74,150  m.  To  be  detected,  the  aircraft  must  be 
above  the  minimum  multipath  angle,  a,  of  .25  degrees. 
Solving  equation  (15)  for  a,  the  angle  above  the  horizon 
becomes: 


,  310  v 

01  =  3111  (5T75  4  0) 

=  .298° 

Aircraft  1  is  inside  the  maximum  detection  range  of  threat 
54,  above  the  minimum  tracking  altitude  of  305  m,  and 
above  the  .25  degree  minimum  multipath  angle.  The  model 
assumes  threat  54  has  detected  aircraft  1.  It  assigns  an 
acquisition  and  tracking  time  of  23.265  seconds.  This  in 
between  the  minimum  and  maximum  times  of  12  and  26  seconds. 
(See  Table  IV.) 

The  next  event  that  occurs  is  event  15  at  the  end 
of  acquisition  and  tracking.  From  the  computer  printout, 
event  15  occurs  at  TNOW  *  55.265  seconds  or  (32  +  23.265). 
At  this  time,  the  site  recomputes  a  SR,  and  determines 
numerous  quantities  including  a  missile  firing  time,  TLO, 
a  missile  impact  time,  TI,  the  range  from  the  site  to  the 
target  at  impact,  RI,  the  probability  of  kill,  PKR,  the 
radar  cross  section  at  impact,  RCS,  the  aspect  angle  at 
impact,  TRCS,  and  the  circular  error  probable,  CEP.  These 
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and  additional  values  will  now  be  calculated  and  compared 
to  the  computer  printout. 


The  aircraft's  x  coordinate  at  TNOW  =  55.265  seconds 
is  calculated: 

x  -  247(55.265  -  30) 

=  6240  m 

Solving  for  the  SR: 

SR  *  -y(26561  -  24592)  2  +  (60,000  -  6240)  2 
*  53796  m 

Because  the  (xA  <  Xg) ,  the  aircraft  is  located 
prior  to  the  closest  approach  point  (see  Figure  11) . 

Solving  equation  (11)  for  8  yields: 

e  (ANGLE)  -  tan'1  [§§§§£§«§!]  -  2.1" 

Solving  equations  (25)  and  (26)  determines  the  firing 
conditions : 


tAC(TIMAC) 


[53796  cos  (sin-1  (gfyjg) )]  cos  (2.1) 


247 


*  217.6  sec 
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tjg,  {TMPO) 


j [53796  cos  (sin-1  (——)  sin (2 . 1)  ] 2+3102p 

759 


■  2.6  sec 


Since  tAC  >  the  firing  will  be  delayed 

(217.6-2.6)  seconds  or  the  missile  launch  time  becomes: 


(TLO)  «  TNOW  +  (217.6  -2.6)  *  270.3  seconds 
Missile  impact  occurs  2.6  seconds  later  or: 

(TI)  *  270.3  +  2.6  *  272.9  seconds 

For  this  encounter,  the  intercept  occurs  at  90 
degrees  (ANG-89.99  degrees).  The  corresponding  aRCS  (SIGMA) 
for  a  SAM  B  at  a  90  degree  aspect  angle  is  25  dB  (see 
Table  V) . 

Aircraft  1  is  wet  (jamming) .  The  terms  A,  B,  and 
C  for  the  CEP  evaluation  (equation  (17))  are: 

A  «  5.62  *  10-6 
B  -  2500 
C  -  232 


The  term  K  for  equation  (20)  is  -51.4  dB.  Solving  equations 
(20)  and  (21)  for  J/S  yields: 

J/SdB  -  65.97  -  25  -  51.4  -  -10.4 


t  /e 


t n-10 .4/10  „ 


nanc 


From  equation  (17)  CEP  can  be  evaluated: 

CEP  -  ^(5.62xl0“6) (.0905) (1988) 2+ (2500) (.0905)+232 
■  21.45  m 

Finally,  solving  equation  (16)  for  the  SAM  B's  P^  yields: 

Pk(PKR)  «  1  -  ,5(43.6/21.45)2  _  >942 

Thus  event  15  schedules  a  missile  launch  time  (event  16) 
for  TNOW  ■  270.3  since  the  P^  value  is  greater  them  the 
threshold  value  of  .05.  The  computed  R^  is  stored  in  the 
model  as  attribute  10  (ATRIB  10) .  The  model  will  evaluate 
the  ATRIB  10  and  a  computed  value  for  R.  in  event  16  to 

•L 

determine  if  the  aircraft  has  maneuvered  since  event  15. 

At  TNOW  -  270.3  aircraft  l's  x  coordinate  is  deter¬ 
mined  as  follows: 

x  «  247(270.3-30) 

The  algorithm  goes  through  the  same  calculations  as  event 
15  to  determine  a  R^.  It  then  compares  this  R^  with  ATRIB 
10.  If  the  Ratio  is  less  than  1.1,  the  model  assumes  the 
aircraft  has  not  turned,  where  Ratio  is  defined  below: 


Ratio 


z 2! 


The  algorithm  goes  through  the  same  calculations  as 

event  15 


x  -  247(270.3  -30)  *  59354  m 
Solving  equation  (14)  for  SR  gives: 


SR  -^(26561-24592) 2 +  (60000-59354) 2 
-  20953  m 


The  R^  can  be  evaluated  as  follows: 


tan  a  « 


4 1 

Ax 


.26961-24592. 

60000-593^4 


8.9° 


ay 


i  cos  8 . 9 
new 


1991  m 


The  model  now  calculates  Ratio: 


Ratio  —  —  1.0015 


Since  the  Ratio  is  less  than  1.1,  the  model  assumes 
the  aircraft  has  not  maneuvered  since  event  15.  (Note, 
since  aircraft  1  is  a  strike  aircraft  and  can  not  turn, 
the  result  is  expected  and  the  small  difference  in  R^  is 
attributed  to  roundoff  error).  The  calculated  J/S,  CEP, 
and  PXR  remain  the  same  as  those  calculated  at  event  15. 
(See  computer  printout.)  Thus  the  missile  is  launched  at 
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this  tine  and  event  17  (missile)  impact  scheduled  for 
TNOW  +  TMFO,  or  270.3  +  2.6  *  272.9  seconds. 

Event  17  occurs  at  the  impact  time.  The  model 
again  checks  the  value  of  Ratio  at  impact.  The  x  coordinate, 
SR,  and  R^  calculations  are  depicted  below: 

x  *  247(272.9  -  30) 

-  59996  m 

SR  »  ^(26561  -  24592) 2  +  (60000  -  59996) 2  +  3102 
*  1994  m 

Ra  -  SR  -  1994  m 

At  this  time  the  model  selects  a  random  number  to 
determine  if  the  aircraft  has  been  destroyed.  If  the  random 
number  is  less  than  the  PKR  the  aircraft  has  been  destroyed 
and  the  computer  printout  annotated. 

Thirty  seconds  later,  event  18  occurs  and  frees  the 
site  to  search  for  another  target.  (See  computer  printout.) 

AAA  Calculations 

The  AAA  geometry  calculations  are  the  same  as  those 
in  the  SAM  portion  of  the  model.  The  difference  between 
the  two  algorithms  is  in  the  manner  the  is  evaluated. 

The  AAA  p.  calculations  is  determined  by  solving  equations 


AAA  Probability  of  Kill  Computer  Printout 


(a)  Subroutine  SEARCH  detects  an  aircraft 


RAOAR  2.  S7/-T  rtACKi.DG  aCcT 

AT 

TRC*  i“. 26292*15  '««• 


(b)  Event  15:  At  the  end  of  Acquisition  and  Tracking 


-v  l2  TM3W=  25  7.26  2C2F15  9  «rPTB(l)e  * . 
RADAR  NUMBER  2. 


i  uw  r\  c  i 

5*  2.SRK=  .7«3'9?'  i  312  73  SI GAU2*  2.2.P9J  97328 

a^-“32‘  32P336-7  T0?L:s  *921  >128*-83£*9  TI*  254.16^3; 
BKS>»  •  361  PKR=  ,  5  3t  d>-  2511 


(Note:  Event  15  calculated  an  immediate  launch, 
is  skipped  and  the  Program  schedules  Event  17.) 


9592 

Event  16 


(c)  Event  17:  Scheduled  Impact  Time  of  the  Round. 


- V  1 7  TNOW*  2f.+  .  18**3  79??2  A~P3B  Cl)  »  , 

G=  2.SRR-  .71355.  “3127?  SIGAAA2*  2  3.35.'  )7326 

VFa  652.  :  3263353 7  TDFL«  .  5  21  3126*  83649  TI*  255*1.  fc;£  8475 

PKSS=  •  361...  8 .125  353  i  PKR-  .15**5346«t251  J 


(d)  Event  18:  Freeing  the  radar  30  seconds  after  impact 

~EV  13  TNG W*  2«died3799"i  A^Rlidjs  *. 


COP* 


(50)  through  (58).  These  are  highlighted  below.  (See 
computer  printout.) 

The  is  determined  as  follows: 

Ri - ** - 

cos  [tan**1  (£y)  ] 

»  24752  -  24108 

cos [tan  1  (2475i-24i08)  1 


■  715  m 


Converting  to  kilometers: 


(SRX) 


Ri 

IMo 


.715  km 


From  equation  (51)  the  velocity  of  the  round  at  impact  is 
calculated  as  follows: 


Vf  (VF) 


930  exp  [-.4965 (.715) ] 


652 


m 

sec 


Equation  (53)  uses  the  Vf  to  determine  the  round ' s  TOF: 

TOF(TOFL)  -  -  2.166  -  .92 

Each  round's  dispersion  pattern  around  the  target,  a, 

becomes: 
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a  «  20R  -  (20)  (.715)  -  14.28  m 
m 

and  o2(SIGAAA2)  »  (14.28) 2  -  203.9  m2 

Evaluating  tha  P.  from  equation  (57)  yields: 

*ss 

p  (PKSS)  -  -  5,11  exD  f  l/C9.8(1.3)  (.92 

Pk  lFKSS)  2 it  (203. 9) +5. 17  exp  \  2\2tt  (203.9H-  5.1 
SS  ' 

-  .004 

This  leads  to  the  evaluation  of  the  overall  for  the 
engagement : 
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P.  (PKR)  -  1  -  (1-.004) 


17 
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****RADAR  positions*** 


I 

4 


y 


6 

T 


*> 


6. 

3 

4 

5 

RADAR  NO 

IPOS  TPOS 
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6 

• 

7 

1 

51480  34199 

8 

2 

50400  24118 

9 

3 

50408  32176 

If 

4 

58488  17963 

11 

5 

50488  32657 

12 

6 

58488  32641 

13 

7 

50488  21523 

14 

8 

50400  19788 

15 

9 

50488  29298 

16 

11 

58488  33852 

17 

tt 

58488  29662 

18 

12 

58480  22993 

19 

13 

51508  38960 

21 

14 

51500  21746 

21 

15 

51588  33384 

22 

16 

51588  31675 

23 

17 

51588  33484 

24 

18 

51508  29685 

25 

19 

51588  21332 

26 

21 

51588  32896 

27 

21 

51588  29177 

28 

22 

51588  29851 

29 

23 

51588  22353 

31 

24 

51588  33660 

31 

25 

52588  28881 

32 

26 

52580  38996 

33 

27 

52508  22884 

34 

28 

52588  41991 

35 

29 

52588  26124 

36 

3f 

52588  36988 

37 

31 

52588  26894 

38 

32 

52588  33888 

39 

33 

52588  26258 

4f 

34 

52588  33854 

41 

35 

52588  29342 

42 

36 

52588  33834 

43 

37 

53488  48618 

44 

38 

53488  24848 

45 

39 

53481  19222 

46 

41 

53488  22361 

47 

41 

53488  35624 

48 

42 

53488  34556 

49 

43 

S3488  27289 

51 

44 

53488  28289 
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be  moved  as  desired.  Friendly  aircraft  can  enter  the  system  at  a 
variety  of  intervals,  altitudes,  and  airspeeds.  WWs  hunt  for 
threats  to  attack  by  searching,  identifying,  locating,  and  then 
launching  their  weapons  at  the  threat.  WW  tactics  can  be  changed 
as  the  requirements  of  the  mission  dictate  or  at  the  desire  of  the 
WW  crew.  Self-protection  jamming  can  be  selected  by  either  WW 
or  attack  aircraft.  Enemy  threats  will  fire  at  an  aircraft  when 
the  aircraft  comes  within  the  threat’s  range  as  long  as  the  threat 
is  not  engaged  with  another  aircraft.  Early  warning  radars  account 
for  threat  radar  command  and  control  functions;  their  control  over 
the  associated  radars  can  be  changed  as  desired. 

Changing  the  WW's  altitude  from  60  meters  to  200  meters  did 
not  effect  friendly  attack  aircraft  survivability.  Leading  the 
attack  force  into  the  threat  area  by  30  seconds  as  opposed  to 
accompanying  the  attack  force  did  not  influence  attack  force 
survivability.  Further  development  of  the  model  to  include  turn¬ 
mode  capability  for  the  WW  weapon  and  a  tactic  for  pre-emptive 
weapons  launch  in  anticipation  of  threat  radar  radiation  is  recom¬ 
mended  . 
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